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“—~] copy of Interim Report No. 2 & attached papers 
: “with this copy of letter, 
UNITED STATES DEPARTMENT OF AGRICULTURE 2nd copy of Interim 
AGRICULTURAL RESEARCH SERVICE Report sent in separate 
S0iL. AND WATER CONSERVATION RESEARCH DIVISION package, 
NORTHWEST WATERSHED RESEARCH CENTER 


_ P.O. BOX: 2700 {4 x 
tags ai \/BOISE,,DAHO 83701 OF ee 
Ste A295 a oe (( : 5 Bec phar ae MELELK ‘ 
sy ie J December 14, 1971 
Fg . 
Subject: Interim Report No, 2 
ARS-BLM Cooperative Studies 
Reynolds Creek Expt. Watershed 
A ieys Glenn H, Lipscomb, Watershed Mgmt, Specialist 


Watershed Staff, P-350 - 
BLM Portland Service Center 

710 N, E, Holladay Street, P.O. Box 3861 
Portland, Oregon 97208 


Two copies of the report have been mailed under separate cover, 

The Interim Report No. 2 has been bound with the 1970 Annual Report, 
Which is considered as a supplement, The 1971 Annual Report will be 
transmitted to you in March 1972, Attached to one of the bound reports 
are copies of all the supporting reports mentioned in the Interim Report, 


- You will note that five of the attached paperg have heen completed by 


Utah State University and one paper by the University of Idaho, The 
papers by Utah State have been supported in part by the BLM funds as 
a cooperative effort in formulating analytical models to describe mois- 
ture movement in soils, The rainfall simulator, which has now been 
field tested, will be used extensively next year to obtain information of 
moisture movement into and through soils for different soils and vegeta- 
tive cover, These data will be used in the analytical models developed 
cooperatively with Utah State University, 


A cooperative endeavor with the University of Idaho, again supported 


in part by BLM funds, has resulted in the formulation of a general 


energy budget model for computing water loss from sagebrush range= 
lands, This model Wi]l be tested in the Reynolds Creek Watershed 
against soil moisture measurements obtained at the eight experimental 


- sites, 


As stated in our review conference of November 18, 1971, which was 
attended by Mr, Bill Meiners, all the items listed under the work plan 
for F¥ 1972 will be essentially completed by June 30, 1972, The next 
meeting to review our Rrogress and to formulate a work plan for 1973 
was scheduled for early April. We will be in correspondence with you 
relative to thig meeting and set a firm date in February, 
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An additional copy of the Interim Report No, 2 is being transmitted to 
Mr, Bill Meiners, District Office, along with a copy of this letter, 


This copy of the report is for review by the District and State offices 
of the BLM, e 
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Table 1. --Cover measurements and erosion condition on 


Figure 1,--Location of experimental Sites, 
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2,.-lats study site. 


3.--Nancy's Gulch study site. 


Figure 4,--Reynolds Mountain (East) study site. 
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Figure 


9.--Sketch of a typical runoff-sediment installation 
for microwatersheds of 1-3 acres, 


6,--Plot collector, tanks, and splitter for runoff 


and sediment measurement. ... 


io CrOpLOt i at,§ by ouit, «for charting vegetation 
changes, Upper Sheep Creek, 


8,--Scaler for gamma system used to determine 
soil moisture, 


9.--Scintillation probe and standard for gamma 


system... 
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Figure 10,--Brush removed treatment, Upper Sheep Creek... 


Figure ll.--Sprayed Plot (2,4,5-T), Upper Sheep Creek. . 


Pigure 12,--Grazed site, Upper Sheep Creek, 
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ARS-BLM COOPERATIVE STUDIES 


REYNOLDS CREEK WATERSHED 
INTRODUCTION 


As stated in ibe Report No, 1 (March, 19.70), the purpose 
or the ARS-BLM cooperative wee Al is to accelerate research on 
precipitation-runoff and sedimentation characteristics of watershed 
areas with specific soil-vegetation Characteristics, together with 
quantifying the effects of certain management and treatment practices 
applied thereto. These joint efforts, during the reporting period, 
have been concerned with the water balance and sediment production 
on sagebrush rangelands as specified in the Agreemet Renewal for 
BY 19°71. 

Research was conducted in the Reynolds Creek Experimental 
Watershed, as documented in the 1970 Annual Report of the Northwest 
Watershed Research Center (attached as supplement) with emphasis 
on infiltration (Res. Outline, Ida-Bo-105, 6), two-dimensional flow 
systems (Res, Outline, Ida-Bo-105, 5), evapotranspiration (Res, Cut- 
line, Ida-Bo-106.1), sedimentation (Res, Outline Ida-Bo-107. 1), and 
cover-soil conditions for different levels of management (Res, Outline, 
Ida-Bo-106,1). The latter Research Outline entitled: "Evaluation 
of Cover Productign, Herbage Yield, and Soil Conditions for Different 


Levels of Vegetation Management," was formalized during the reporting 
g & g 


2 


period, Progress in the development of study sites and a summari- 


zation of research investigations are included in this Report. 


SMOG una 

Locations of the eight study sites, developed to study the 
interrelationships of cover, herbage yield, soil Ondo ee infilira- 
tion, evapotranspiration, and sediment production under different 
levels of vegetation Mmahagement are shown in Figure l. Site infor- 
mation was included in Table 1 of Interim Report Not 1970). © Dhe 
Nettleton site was selected, fenced, and partially instrumented during 
the past year, The development of a study site on granitic soils 
Was postponed until next year. | G&G 

Feneing for livestock exclosure or control was completed at 
all active sites prior to the 1971 grazing season and close-mesh fenc- 
ing to exclude rabbits was installed at the Lakebed Flats site,  Phote- 
graphs of three of the eight study sites, Lakebed Flats, Nancy's Gulch, 
and East Reynolds Mountain, are included as Figures 2,3, and 4, 
respectively. 

INSTRUMENTATION 

Runoff and Sediment 

Instrumentation for runoff se sediment measurement on 2 to 
3-acre micro-watersheds was installed at the Lakebed Flats and 


Nancy's Gulch sparse cover sites, Figure 5 is a sketch of a typical & 
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Figure 1. Location of experimental sites, 
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FIGURE 3. Nancy's Gulch Study Site, 
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FIGURE 4, 


Reynolds Mountain (East) 
Study Site. 
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2 
runoif-sediment installation which includes: (1) a fiberglass cutoff 
wall extending to backhoe excavating depth, (2) a catchment tank for 
coarse sediments, (3) a 30-degree V-notch weir and Water-level 
recorder to measure runoff, and (4) a series of single-stage sus- 
pended sediment samplers with intakes at regular depth intervals to 
sample Eamon concentrations at varying flow rates. » Also, a i4- 
it. by 72-ft. grazed, sparse vegetative cover plot at the Nancy site 
was instrumented with borders, a collector, tanks, a water-level 
recorder ang splitter (Figure 6) to measure runoff and sediment 
production during storm events. Similar installations will be installed 
at the four other sites with sparse vegetative cover, Additional 
runoif and sedimentation facilities at Upper Sheep Creek and Reynolds 
Mountain were described in Interim Report No. 1 (1970). 

Soils Moisture and Vegetation 

Installation of access tubes for measurement of soil moisture 
using both the gamma and_neutron systems was completed at all sites 
except the sparse vegetation study site at Upper Sheep Creek, 

Round portable exclosures seven feet in diameter were con- 
structed to protect herbage yield plots from grazing while plants 


‘were maturing, Triangular exclosures similar to steel gate panels 


were also constructed for the same use on dense vegetation sites, 
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RTiGURE 62 Plot Collector, Tanks, Recorder and Splitter 
for Runoff and Sediment Measurement, 
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Instrumentation for measurement of temperature, relative 
humidity, and wind were installed at Nancy's Gulch and Upper Sheep 
Creek, Such instrumentation was already in operation at Lower 
Sheep Creek and Reynolds Mountain, 

INVESTIGATIONS 
Infiltration and Flow in Porous Wadia tl 

Studies concerned with infiltration and flow in porous media are 
discussed in detail in the 1970 Annual Report, Research Outlines 
ida-Bo-105.5 and Ida-Bo-105, 6, A summarization o? these studies 
is presented in the following paragraphs: 

ihe.accuracy and utility of computer solutions of the aint 
differential equations that describe unsteady-state, axisymmetric infil- 
tration were evaluated for use as 2 vrediction model; Laboratory 
tests were conducted on large soil-cores by using a rainfall simulator- 
gamma-probe- infiltrometer, This instrumentation was discussed in 
Interim oe Noses) (1970), 

A mathematical model describing transient flow of water from 


circular infiltrometers was formulated and analysis of preliminary 


| 2 
humerical solutions was accomplished by Jeppson (1), 2) Pertinent 


Re a ee eee eee REET CPO MINOR ome 


1/ Research on infiltration and flow in porous media was conducted 


cooperatively. with Utah State University and Oregon State Univer- 
Sity. 


2/ Numerals in parentheses refer to publications, page 23. 


features that describe the flow characteristics obtained irom 32 
solutions for varying initial conditions and for 12 soil types were 
summarized in tables and figures, Relationships between depth of 
moisture penetration and the following characteristics were developed: 
(1) lateral movement of the wetting front, (2) rate of moisture applica- 
tion, and (3) initial soil-moisture tension. 

In a companion report, Jeppson (3) described the determination 
of hydraulic conductivity-capillary Sues aes relationships from satura- 
tion-capillary pressure data for soils. The Burdine Theory was used 
to obtain relative permeability from the pressure-saturation data, 

The report describes a computer program used to evaluate Burdine . 
Lo 
Integrals and describes input to and output from the program, 

Solution to transient, vertical, moisture movement based upon 
Saturation-capillary pressure data and a modified Burdine Theory 
was described by Jeppson (4) in a later report. This formulation 
and solution method is consistent with the solution method developed 
earlierWill 2neonege three-dimensional, axisymmetric movement of 
water applied at the ground surface within a circular infiltrometer, 
This formulation was deemed necessary, so that comparisons of the 
solution results from the two different cases would indicate quanti- 
tative effects on the flow pattern of the component of radial, moisture 


movement, 
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Vertical moisture-movement solutions were obtained (4) for 
several initial moisture contents and for several application rates 
for a sample of disturbed soil taken from the summit area of the 
Reynolds Creek Experimental Watershed. Saturation-capillary pres- 
sure data for this soil were obtained in the laboratory. In addition, 
laboratory measurements of the hydraulic conductivity corresponding 
to a number of Pep ieheiouie ssarbs were obtained. Using the satura- 
tion-capillary pressure data in the Burdine Equations for evaluating 
hydraulic conductivity gave good agreement with the laboratory meas- 
urements of hydraulic conductivity. 

In a later report that is not included in the Res re Report, Wei 
and Jeppson (5) obtained solutions to the problem of steady-state 
infiltration of moisture that moves through partially saturated, homo- 
geneous soils toward a water table from a circular area of application, 
Solution techniques were formulated, and solution results presented, 
The solutions indicate that significant radial movement of moisture 
(spreading effect) occurs and causes higher infiltration rates at the 
edge of the horizontal source circle than near the center, The infi]- 
tration rate is closely related to various soil parameters that char- 
acterize the hydraulic properties of soils, Also presented are 
several distributions of the relative permeability or effective satura- 


tion .on the surface, along the axis of symmetry, and on the plane 


ize 
. «& 
including the axis of symmetry. In addition, the report indicates 
how the distributions of relative permeability or of effective satura- 
tion are related to the soil parameters, 
During the ee of 1971, saturation-capillary pressure data 
were obtained from four 32-inch diameter, undisturbed soil cores. 
These cores were taken from different locations on the Reynolds Creek 
Experimental Watershed and included different soil types. Small cores) 
2-inch in diameter, were extracted from the large cores for obtaining 
conductivity-pressure data in the laboratory. In addition, a mathe- 
matical model of transient, vertical moisture movement through 
layered soils was essentially completed. . 
In another study the adequacy of a numerical solution of the. . 
steady-state, two-dimensional flow system resulting from infiltration 
of water from melting snow was tested. Jeppson and Schreiber (6) 
described the mathematical. model, and Schreiber, et al. (7) evaluated 
the model by field data obtained from a northern-slope section of the 
Upper Sheep Creek jigs ee 
The model (6) provides for soil heterogeneity, soil anisotropy, 
and partially saturated-saturated flow regions. A representative 
solution, illustrates the cane bitin and flexibility of the model. 
Field measurements obtained by Schreiber, et al. (7) during 
1970 included soil moisture, soil capillary pressure, piezometric 


e 


fluctuations recipitation, snowpack volume, snowpack water content. 
a 3 3 Pp ? 








snowmelt, and streamflow. Soil properties were determined in the 
laboratory. 

Mathematical solutions and laboratory data indicated that saturated 
hydraulic conductivity SY Are ent, varies with distance up slope 
from the stream channel, Other parameter variations, such as 
anisotropy, could He Specified to bring solution results in agreement 
With observed conditions, 

A different solution technique, Kirchhoff Transformation, for 
the solution of transient flow of water from a circular infiltrometer 
was described in a supplemental report by Jeppson (2). it was sug- 
gested that the numerical solution might be improved for situations 
containing iste of nearly saturated flow by means Of the Kirchhotr 


Transformation, Example solution results are presented in the report. 
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Studies concerned with evapotranspiration from sagebrush range 


Sites are summarized in the 1970 Annual Report (Res, Outline, Ida- 


ae 


\ 


Bo-106.1.) Further information on the studies was reported by Belt(8). 


Estimates of evapotranspiration from sagebrush rangelands by 
the Bowen ratio method and by energy balance-combination equations 
were made, using data vet iese by ee eee instrumentation 
during the summers of 1969 and 1970. The Bowen ratio estimates 
were found to be smaller than those obtained by other methods. 
3 / Research on evapotranspiration from sagebrush rangelands was 
conducted cooperatively with the University of Idaho, 


oe ee ee SS a Se ae ae eee a pe a ae 


Aa rire mit Ee st Sh Bt ne Fak 


CS 


The apparent difficulty in using the Bowen ratio method in semi- 


arid environments lies in the assumption of equality in the transfer 


coefficients of heat and water vapor and in Obtaining reliable meas- 


uremenis of very small] vapor p 


tessure gradients. The energy 


balance-combination equation for computing evapotranspiration that 


requires humidity data in addition to surface temperature and radia- 


tion required in the basic SPT Shag 


Since the computations were fou 


y balance equation is not applicable 


na to. be redundant, 


Sufficient data have been obtained to evaluate the inequality of 


the transfer coefficients of heat-vapor and momentum, This step 


will make it possible to obtain independent estimates of evapotirans- 


resistance concept, to compute 


vegetated rangelands, 


Runoff and Sedimentation 
entation 


Cc 


Piration for testing a proposed energy balance equation, based on the 


evapotranspiration from sparsely 


During the 1969-70 water year the only important runoff produc- 


ing events occurred during the January 21-27, 1970, storm and the 


May-June snowmelt period, January precipitation Varied, fromy2.96 


inches at 3915 feet elevation to 


the January precipitation reache 


9.952 inches at 6800 feet, Although 


d near record amounts, only minor 


flooding occurred and the high elevation snow absorbed most of the 


rain, Also, sediment Production from the January runoff was low 


at the Reynolds Mountain and Upper Sheep Creek Sites, 
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The snowmelt eee was somewhat later than normal and melt 
runoff from higher elevations Meee mainly in late May and early 
June, Snowmelt runoti and sediment production were about normal. 
Further details are contained in the 1970 Annual Report (Res. Outline 
Ida-Bo-107.1). Processing and compiling of runoff and sediment 
data for the 1970-71 water year are scheduled for completion in 
December 1971 and will show above normal winter precipitation and 


snow accumulation, 


Soil Moisture and Vegetation 


>) 


Duplicate micro-plots of 3 it. by 3 ft. were established at all 
study sites for each type of vegetation treatment, Initial surveys 
of these plots were wey ENE regarding identification of species of 
grasses and forbs, and the measurement of the basal area of grasses 
present, 


These plots were also charted and photographed. Any fluctua- 


tions in either species composition or basal area of grass clumps 


caused by treatments imposed will be compared with this initial 


se S  eae yey 
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information. Figure 7 shows one of the microplots on the brush- 


ro 
as 


removed treatment at Upper Sheep Creek. Clippings were also taken 
from all sites for measurement of herbage yield. 

Soil moisture measurements were obtained at three of the study 
Sites: Lakebed Flats, Nancy's Gulch, and Upper Sheep Creek, dur- 
ing the period of this report, The neutron depth gauge was used for 
measurements below six inches. Use of the gamma density probe 
was initiated in 1971 for the purpose of following moisture change in 
the 0-6 inch zone. Figure 8 shows a scaler for gamma measurements 


being used at Upper Sheep Creek while Figure 9 shows the scintillation 


probe of the gamma system. 


tudies were initiated in 1969 to evaluate the effects of various 
brush treatments on herbage yield. This initial study was incorpora- 
ted into noe Outline, Ida-Bo-105.4. Results.from the study in 1970 
are reported in the 1970 Annual Report, page 7-7. | Figures 10 and 


ll show the sites from-which brush was removed and sprayed, res- 


—— 


pectively, with all grazing animals excluded, Figure 12 shows the 
grazed treatment. 

Cover transects were recorded for all sites and each site was 
assigned an erosion condition class derived from the BLM Soil Sur- 
face Factor rating sheet, COPS Re Agee ne and erosion condition 


are given in Table 1, 
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for Charting Vegetation Changes 
(Upper Sheep Creek) 





. y FIGURE, 8 Scaler for Gamma System Used 
\ | to Determine Soil Moisture, 
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FIGURE 9, 


FIGURE 10, 











Scintillation Probe and Standards 
for Gamma System. 
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Upper Sheep Creek, 


ata See ER ee Rae Pt OE ee = — 


ae Ae Sa am 


' 
7 7 oe oe | a. 
7 ee me 
a irs Pia _ at 
oars Any 
~abnei® bus sdord cobelitiatod 
eteye ekaumai) tot 





FIGURE 11. 


FIGURE 12, 











Sprayed Plot (2,4, 5-T) 
Upper Sheep Creek, 


Grazed Site 
Upper Sheep Creek, 
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TABLE 1: __GOVER MEASUREMENTS AND EROSION CONDITION ON STUDY SITES 
Smal Large Bare Erosion 
Site and Treatment Vegetation Litter Rock Rock Ground Condition 
th dh ts te b Class 
Lakebed Flats | 
Grazed 40.6 7 24.9 4.7 0,0 29.9 Slight 
Exclosure - Not grazed 39.4 2554 4.9 0,0 30.3 Stable 


Nancy's Gulch 


a8. 6a phe 2 le od Vy. 21.1 Slight 


Grazed os Je Ver 

Exclosure - No brush treatment — 43.4 2bek 14.3 0.7 16.4 Slight 

Exclosure - Sprayed 41.3 25k 13. 1 0. 6 Ls Slight 

Exclosure - Sage removed ; 28.4 30.9 12.9 Eee 26.6 
Lower Sheep Creek 

Grazed 48.4 19.1 24,7 0.4 1.3 Moderate 

Exclosure - Not grazed 55.6 18.4 Peo 0.6 ee Moderate 
Upper Sheep Creek (North) } | 

Grazed | 40.7 15.1 23.4 1s eat Moderate 

Exclosure - Not grazed 52.0 15.6 19,5 1.3 = i2 8 Moderate 
Reynolds Mountain (West) 

Grazed 50.6 Lies 33,0 0.0 ek Stable 

Exclosure - Not grazed 48.4 11,1 35.8 0.3 4.3 Stable 
Upper Sheep Creek (South) 

Grazed 78.9 19.8 0.0 0.0 132 Stable 

Exclosure - No brush treatment 87.6 Sept as “0.0 0,0 0.4 Stable 

Exclosure - Sprayed _ 89.4 a. 6 0.0 0.0 i.v Stable 

Exclosure - Sage removed 89.4 8.6 0.0 0.0 2a) Stable 
Reynolds Mountain (East) 

Grazed $6.3 1224 tha! 0.0 Tai Stable 

Exclosure - No brush treatment 83.7 1339 0.1 0.0 ces Stable 

Exclosure - Sprayed 83.6 14.7 0.4 0.0 Se Stable 

Exclosure - Sage removed 37.0 38.4 0.4 One 4.0 Stable 
Nettletons 

Grazed 68.1 2.0 oe ee 4.3 Stable 

Exclosure - Not grazed 74.4 4 1.4 0.3 Sect Stable 
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21 
SUMMA RY 


Data collection was initiated on eight special study sites (Figure 
1, page 3). These sites were developed to study the interrelation- 
ship of cover, herbage yield, soil condition, infiltration, evapotrans- . 
piration, and sediment production under different treatments, ‘Tabu- 
lations of cover measurements and erosion condition classes for 
separate treatments are given in Table 1, page 20. 

Instrumentation yet to be completed on the study sites fecienic 

soil moisture access tubes on Upper Sheep Creek (North) and per- 


manent runoff-sediment measurement plots on four sparsely vegetated 


_ sites, This work is scheduled for completion during FY 1972 with 


one additional study site to be installed on granitic soils, 

The cooperative Meeeioations of a two-dimensional, steady-state 
watershed flow system was Benassy completed and reported in 
publications 6 and 7, ther investigations concerned with infiltra- 
tion led to the development of a mathematical model to describe the 
three-dimensional axisymmetric flow of water from an infiltrometer 
through partially saturated soil and a model for calculating the one- 
dimensional, vertical movement of water in soil, These Uy aaet and 
solution techniques are reported in publications 1,2,0, 4, and 5, 
Laboratory and field data have been obtained to test the applicability 


of these models. 


22 
Preliminary analysis of data on evapotranspiration phtaned 

by micrometeorological methods indicates that by. considering the 
inequality of the transfer coefficients for heat-vapor and momentum 
such methods may be used to obtain independent estimates of evapo- 
transpiration, Data already obtained will be used to test a proposed 
energy-balance equation to compute evapotranspiration from sparsely 
vegetated rangelands, A summary of evapotranspiration data is 


found in publication 8, 
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23. 
PUBLICATIONS 


Jeppson, R.W., 1970. Transient flow of water from infiltrometers-- 
formulation of mathematical model and preliminary numerical solu- 
tions and analyses of results, Project report PRWG-59c-2, Utah 
Water Research Laboratory, Utah State University, Logan, 


Jeppson, R.W. 1970. Formulation and solution of transient flow 
of water from an infiltrometer using the Kirchhoff Transformation. 
Project report PRWG-59c-3. Utah Water Research Laboratory, 
Utah State University, Logan, 


Jeppson, R.W, 1970, Determination of hydraulic conductivity - 
Capillary pressure relationship from Saturation-capillary pressure 
data from soils, Project report PRWG-59c-4. tah Water Re- 
search Laboratory, Utah State University, Logan, 


Jeppson, R.W, 1970, Solution to transient vertical moisture 
movement based upon Saturation-capillary pressure data and a 
modified Burdine Theory, Project report PRWG-59c-5. Utah 
Water Research Laboratory, Utah State University, Logan, 


Wei, C-Y¥., and Jeppson, R.W. 1971. Finite difference solutions 
of axisymmetric infiltration through partially saturated porous media, 
Project report PRWG-59c-6, Utah Water Research Laboratory, Utah © 
State University, Logan, ; 


Jeppson, R.W. and Schreiber, D.L, 1972, Solution of a two- 
dimensional, steady-state watershed flow System: Part I, Des- 
cription of mathematical model. A portion of Paper No, 71-255 
Presented at the Annual Meeting of the American Society of Agri- 
cultural Engineers at Washington State University, Pullman.  Sub- 
mitted in July 1971 for possible publication in the Transactions of 
the ASAE, 


pcnreiber 1)! Ty a. Jeppson, R.W., Stephenson, G.R,, Johnson, C.W,, 
Cox, (hb, Moy and Schumaker, G.A, 1971. Solution of a two-dimen- 
Sional, steady-state watershed flow system: Part II, Evaluation 

by field data. A portion of Paper No, 71-255 presented at the 
Annual Meeting of the American Society of Agricultural Engineers 

at Washington State University, Pullman. Submitted in July 1971 

for possible publication in the Transactions of the ASAE, 


Bert. totus) Opring evapotranspiration from low sagebrush 
range in southern Idaho. Research Project Technical Report, Project 
A-014-Ida, Water Resources Research Institute, University of Idaho, 
Moscow, Idaho (October), 
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PREFACE 


STATEMENT OF L.ICSION 





Congress, in 1960, established the Northwest tiatershed 
Research Center to gain basic information on runoff character- 
istics, including water yield, from plateau and foothill 
grazing areas of the Northwest and basic information on runoff 
and sediment problems of the Northwest wheat-producing areas, 


The research mission is to gain a better understanding of the 
wien ms 


role of the land anc the influences of vegetation, climate, 
and land management on the movement of water and sediment. 





Objectives o£ the research program concerning the hydrologic 
4 te = 
performance of range and agricultural iands are: 


1. To measure water and sediment yields, evapotrans- 
piration, changes in storage of soil moisture and 
ground water, precipitation and cther climatic 
variables on selected study areas at representative 
locaticns in the Northvweste 


2e To define the influence of soil, vegetation, 
climate, geology, and land management on the 
disposition of precipitation and on sediment 
yield for the study areas. 


3. To develop a hydrologic model for the prediction 
of Water and sediment yields from range and 
agricultural lands in terms of readily obtainable 
climatological and physical data for different 
management levels. 


Research is conducted in cooperation with the Experiment: 
Stations of Idaho, Oregon, and “Washington, and with the Ovyhee 
Soil Conservation District, the Soil Conservation Service, and 
the Sureau of Land Management. esearch ¥indings are docu- 
mented in the technical literature, and hydrologic data are 
released following evaluation for use cy other groups or 
individuals. 


"ie Annual Resort for the Year 1970 


LTITRODYCTORY PATLSRIAL 


This annual resort is a combined effort of the professionel, 


noneprofessional, and clerical staff 
shed kesearch Center, cise, Idano,. 


Ae STATION OPERATION 





t “personnel: 


the Northwest vViatere- 


ALL personnel employed during the report- 


ing year are included in the following list, except for 
temporary exsloyees on short-term apocintments. 


Professional 





Mame a 








We Re Hammon Research Hydrauiic 
Engineer (RIL and 


Center Director) 


Cc. W. Johnson Research Hydraulic 


Engineer 


(a eae he bt 


Research Hydraulic 
Engineer 


Researcn Soil 
Scientist 
Starting 7-26-70 


Ge. A, Schumaker 


Ge R.e Stephenson Geolo is 
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oO 


as 
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io tM 
e 6 


j+ 
e 


nesearcn 
Responsicilities 





Hydrologic models. 
Precipitation, 
LAST rAt Ons ae 
evanotranspiration. 
Center administra-= 
tion and technical 
suoervisione é 


Runorx investiga=-* 
tions, 
Sedimentation. 
Weir construction. 


Snowmelt studies. 
Evapotranspiration 
stucies. 


SUrTaAcS Lunoli, 
ii ritration, and 
ground watere 
Hydrologic models. 
Sedimentation. 


Soil and vegeta-# 
tion studies. 
YVegetative response 
to management and 
treatment. 

Soil moisture. 


Hydrogeoclogic ine 
vestigations. 
GeonorphoLlogy « 
studieS.e 
Geochemistry of 
ground water. 
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Name 
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Mame 





Nee Ot Lot 

















sotanist 
Lich attending Colo- 
radco State University 
een AgGronomist 
TO 3-21-79 
Matnematician 
To G-6-70 
Clerical 

Title 

Clerk-Stenogranher 


ee: rales 2 - my 
Clerk-Stenograpnher 


General Schecule 








Engineering fechni-+- 
cian 


Engineering 
cian 
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L=2 
<esearch 
wesponsibilities 
Hegetation intlu-e 

ences, 
PRecien tac Lott 
chervacterretics, 
Snov desosition,. 


tative response to 
management anda 
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Secretarial, 
Jooxkeeninge 
Personnel, 


ctensgraphy 
Pitino cand travel. 
Ciprer. 
otenography 
rrocurement,. 
Recent i ond Sih 


Duties 


Data collection 
eB reduction. 


anunofs data 
collection. 
Sediment sazcling. 























Suoprofrfessional Continued) “a 
General Schedule 
Nae Title Duties 
D.C. Ropertson Engineering Techni- L, Data collection 
cian and reduction, 
2e instrumentation. 
aetnt ef melic ee) Engineering Techni- 1. Hycdrogeolsgical 
cian Gata collection 
and analysis 
2. Lata Dprocurenent 
scheduling. 
3, Drattigicg 
M. D. isurgess Blectronic Techni- Instrumentation, 
cian 
To 10-38-70 
Wage Grade 
Naiae Title Duties 
is 
Me. I. Gasser Soil Research Helper Data reduction. 
To 9-26-70 
R. F. Hermann Carpenter Leader Construction, 
Re M. Hoagland Automotive Mechanic Automotive renairs,. 
Lee Perkins Maintenance Worker 1. ®ield instrumenta- 
CLIN 


2e Maintenance, 


K. We. Trautman Engineering Equin- 1. 2quiprient main- 
ment Operator tenance, 


Ze ROaG construction, 


C. *,. Ward Carpenter Construction. 


on to the foregoing, a number of temoorary employees 

to the extent of 502 wan-days in construction, 152 

in instrumentation, 350 man-days in data acquisition 
ays clerical, or. .a total of,1294 
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g “a 22) IFAC Li ties 


The research offices are located in a leased building at 306 

N,. Sth Street, Soise, Idano, Facilities located at Gowen Field 
(Boise Airport) consist of a shop housed in a quonset building 
and a soils, geology, and instrument laboratory housed in a 
prefabricated steel building. The Reynolds Creek Experimental 
Watershed field headquarters is located in the watershed at 
Reynolds, idaho, 50 miles scuthwest of mbOLSe, and consists of 

a quonset buiiding and storage sheds on 3 acres of leased land, 


3e iInstrumentation— 1/ 


A tacular summary of instrumentation is found in Table 1, 
(pe 1-5). Location maps are shovm in Figures 1 and 2 for 


Reynolds Creek ‘iatershed and Raboit Creek ‘latershed, respec- 
tively. 
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Instrumentation additi 
(1) installation of 
Reynolds Mountain Via 
Gulch study area and Wiice 
Gevelopment and ia steaias 
Reynolds Mountain West 


}4e 
QO 
Bey 


the renorting year include: 

gamma prove access tubes at 

heen Creek Jiatersneds, Nancy 

studysareagys(2)..soring — 

water Measuring devices at 

exrshec anc near the Summit Water 

: o shed; (3) installation a Universal Surface Precipitation 

$ Rs Gage, snow EGE and precipitation gages in the drainage 

4 basin of the test Fork of veynolds Creek; (4) an improved 
solarimeter at ihe Lower» Sneep Creek tieather Station; (5) 
gravity suspended sediment samplers at Moscow, fdaho (Thomoson ~ 
pce Shas) and the Summit “iatershed trib utary to Reynolds 
Creek; (6) a runoff-sediment plot and tanks at the Nancy Guich 
study area; and (7) an infiltrometer rainfall simulator test 
facility for undisturbed soil cores at the eynolds Field 
Headquarters, Aliso, several study sites were fenced, 
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LA Locations are found by use of the federal system of land 
divisions where each section in the Neynolds Creek Water~ 
shed numbered as shown in Figure 1 *rom the upner left 
(northwest) corner in rowe of 10, In the Rabbit.Creek 
watershed, Figure 2, the reguiar.secticon NuMmberS Are uSeGe 
Bach section is subdivided into 100 squares, as indicated 
in the subdivided section in the unper right corner of 
Fig gure 1, Within a numbered section, a location is found 

doy moving dowm any section to the row, and then across to 
the proner colur & precipi ation gage is designated 

by a 6—digit pmipert The first three digits renresent 

the section anc the last two digits represent the location 
of the smallest subdivision within the sections 
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Cooperative 
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Ground-Water 

Study | Basin (86) 
Nancy - ‘Galen 
Study Area (98) 
Sheep Creek (117) 
Watersheds (138) 
Climatological 
Stations C127) 
Reynolds Mt. 

Study Basins (166 & 176) 
Climatological 
Stations (176 
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(1) The cooperative work in the Keynolds Creek Water- 
shed is for comprehensive studies on the effects 
of vegetation uncn the hydrologic performance 
of watersheds. The 3LM furnishes supplemental 
funds and the required esperimental sites on 

| public lands as assistance in conducting these 

studies, The funds have been used primarily in 
support of investigations under esearch Cutlines 
Ida-30-105.4, 105.5, 105.6, and 106.1. 

(2) Data on the Rabpit Creek siatershed (Frail Lands 
area) are collected, processed, and reported : 
with onerating funds suppiied by the SLM DLSTIACt 
Office, Boise, idaho. A map of the Rabbit Creek 
Watershed, showing the data network, is included 
aS Figure 2, F.7.< E 


be University of Idaho 


(1) Studies of the ground-water flow systen and 
ae water balance for an irrigated £iecld (Research 
§ Cay Outline Ida-i0-103.3) has progressed to the final 
analysis stage. A Doctor's thesis by Mr. David 
Allman and papers covering this work will be 
prepared in 1971 and the Outline terminated. 


(2) Field data on evapotranspiration obtained by 
use o£ micrometeorological instrumentation has 
been processed by Mr. George Selt and prelimi- 
nary analyses made (Zesearch Outline Ida-50- 
106.1). At least one paper will be prepared 
during 1971 covering the analysis or data 
obtained during 1969, fdditional field data 
Will be obtained in 1971. 


() 


(3) Collection of climatological, scil moisture, 
runoff and sediment data from the Thompson 
Watershed near Moscow, Idaho continued through 
1970 under direction of Dr. Myron Molnau,. Also, 
3 tivo gravity samplers for suspended sediment were 
designed and fabricated at the University for 
use on the Thompson and Reynolds Creek Watershedse 


Ce University of Nevada 


Data collection and analysis of the ionic con- 


| 3 ms centrations in a ground-water flow system is 
A She approaching completion. (Research Outline © 


Ida-B0-103.4.) The chemical analysis has been 
; performed by Dr. John Snarp of the Desert Ree 
search Institutee 
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de Oregon State University 


Worx has continued on development of techniques 
and instrumentation ror determining hydraulic 
properties of soils in the field by Dr. xoyal 
Brooks and Mr. Neil Biggs. (This work is re-= 
ported under Kesearch Outline Ida-o-105.6.) 
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e,. Utah State | 


Investigations on the flow of water from an 
infiltrometer and two-dimensional flow systems 
have been conducted by Dr. Roland Jepnson and 
Mx. Chi-Yuan Wei. Xevorts have been prepared 
covering analytical solutions and leboratory 
data obtained for verification. (Work reported 
under Research Cutlines Ida-5o0-105.5 and 105.6.) 
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Hydrologic data c: rom the instrumentation as 
listed in Table 1, p. 1-5, and as located in Figures 
1 and 2 on Pages 1-6 and 1-7, respectively, have been 
tabulated or entered on punch cards for orocessing 
and analysis. These data for the water year of 1970 
are included as an Appendix (attached to reserve- 
working copies of this Report) containing information 
on runoff, precivitation, snow water content, soil 
moisture, sedimentation, water-table depths, wind, 
solar radiation, and humidity. 
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REPORTS BY INDIVIDUAL EXPERIMENTS 
CRIS Work Unit No. SWC-Ol-f£o-1 Code No. Ida-Bo-100.1 
Ditle: Snow accumulation and melt on study areas. 


Location: Northwest \¥/atershed Research Center, 
Boise, Idaho. 


Personnel: Lloyd NM. Cox and Freeman IM. Smith. 
Date of Initiation: August 1, 1961. 


Expected Termination: Originally planned, November 1968. 
Recommendation, November 1971. 


Objectives 


1. To determine the physical and meteorological factors con 
tributing to nonuniformity of snow accumulation in SR ces 
covered study basins on mountainous terrane, 


i 


ie) 
s 


To determine the influence of the controlling physical and 
meteorological factors on snowmelt from the above areas. 


( Need for study: 


A substantial proportion of the runoff from the sagebrush rangelands 

of the Northwest has its origin in the melting of snow. Any proposals 
| to improve the quantity or timing — of flow from snow-fed streams by 

| manipulation of vegetation or by other practices require that the be- 

| havior of snow be well understood. There has been little research on 
\ the behavior of snow in shrub areas anywhere, and almost none in the 


Sagebrush areas of the Northwest. 


Destructive late winter and spring floods in the Northwest frequently 
originate from rapid melting of snow at low elevations characteristic 

of the sagebrush zone. Although there is little likelihood of modifying 
‘Snowmelt rates enough to alleviate this threat, knowledge about the 
behavior of snow in the sagebrush zone will be helpful in devising better 
warning and forecasting techniques that may reduce the danger to life 


and property from snowmelt floods. 





Design of the Experiment and Procedure to be Followed: 


Because of the exploratory nature of this experiment, no ordinary sta- 
tisticol design is applicable to the study. However, the general plan 
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is *to choose two unit source areas--both shrub-covered--and at ran- c 
domly located points, to determine the following variables: snow depth 
and water content, slope, aspect, vegetation height, total and net solar 
radintion, runoff, temperature, wind, and humidity. Multiple regres- 
sion analysis will be used to identify the characteristics contributing to 
heavy snow accumulation, | 


The areas of heavy and light snow ceccumulation will be delineated by 
the use of photogrammetric techniques, The photogrammetric data 
will be used to construct cross sections to show snow depth in relation 
to ground-surface features, 


Experimentcl Data and Observations: 


Destructive flooding from melting snow is usually associated with a 
rein-on-snow event, The rain by itself adds very little heat to the 

snow to produce melting. v/hat usually happens is that the condensa- 
tion heat and also convective heat exchanges are also acting at the 

time of rainfall, thus contributing to snowmelt. If the snow is quite 

cold and the crystals quite fine, then the snow can be termed as dry 

and has the capacity to store water. A rain-on-snow event under these 
conditions may not produce any runoff until the storage capacity of the 
snow has been satisfied. Much depends, of course, on the existing — 
meteorological conditions and the depth of the snowpack, 


Figure i is a graph depicting - rain-on- snow event with a snow depth 
of ll inches and water equivalent of 1.28 inches. The curve at the top 
of the ficure represents the changes of air tempernture with time during 
this event. The next curve shows the changes in water equivalent with 
time as measured by 2 universal gage and the melt from this gage is 
represented by the second curve from the bottom, ‘The bottom curve 
is a plot of the precipitation 2s collected with an é-inch unshielded rain 
gage located adjacent to the site. 


Precipitction was observed as rain from 0800 November 23, 1970, until 
0100 on November 25, 1970. ven though the rainfall was not constant, 
water continued to flow out from the bottom of the snowpack at a fairly 
constant rote until 1600 on November 24, 1970. After this time, bursts 
of precipitation produced surges in the water being collected at the bot- 
tom of the snowpack, indicating thet the storage capacity of the snow had 
been reached, With the cessation of rainfall at 0100 on November 25, 
1970, water continued to drain from the snowpack until 1000 on November 
25, 1970. Snow started to fall at 0700 on November 25, 1970, thus 
producing a2 subsequent increase in the water equivalent measurement, € 
f 


l1/ Figure 1 follows page 2-4, 





an 
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The total precipitation collectec by the unshielded rain gage for this 
event starting at 1000 on November 22, 1970, through 1200 on November 
25, 1870, was 2,17 inches, while the shielded gage caught 2.4 inches, 
The water equivalent measurement showed a net incrense of 1,10 inches, 
using o tentative calibration. The melt collected curing this time 
period wns 1.44 inches of water. | 


Comments, Interpretation and Future Plans: 
Pe es ee ne eis 


Data collected during the rain-on-snow event showed that the net change 
in water equivalent for the event resulted in an incrense of 1,10 inches 
of water. Total melt collected for this same period was 1.44 inches of 
water. Cf the 2,17 inches of precipitation collected by the unshielded 
rain goge, 1.82 inches was observed as rain, the remaining was snow. 
This indicates that at least 0.38 inches of rain (1.82 -1.44 = 0. 38 
inches) was stored within the snowpack. 


Data collected from another snow study site during this same rain-on- ~ 
snow event did produce similar results in the total nmount of melt 
water collected. Unfortunately the rain gage at this other site was 

not in operation and no precipitation was collected, 


Previous studies (Research Cutline Ida-Bo-102.6) have indicated that 
the precipitation caught in an unshielded rain gage is not representative 
of the actual catch. If the calibration coefficient, 2s determined from 
these studies, is used to calculate the actual precipitation using the 
unshielded total of 2,17 inches and the shielded totel of 2,40 inches then 
a value of 2,67 inches is the totel precipitation for this event, This 
would mean that about .13 inches more precipitation fell than could be 
accounted for by the increases of water equivalent ond melt water col- 
lected (1.10 + 1.44 = 2.54 inches) ty the universal gnge. This is a 
very close agreement, considering that a tentative cnlibration was used 
for converting the pressure recording on the universal gage to water 
equivalent and that the calibration coefficient for the rain gages may be 
slightly different than a value of 2,0, 


One other source of heat that enn produce appreciable snowmelt comes 
from the long wave radiation that is radiated from the low warm cloud 
cover during the spring melt senson. An effort was made to collect 
net radiation data during the previous spring snowmelt season, but this 
effort was unsuccessful because of instrument failures. Duplication of 
criticel cormponents in instrument recording systems is 2 must if con- 
tinuous deta are to be collected xt remote snow research sites. 
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Sufficient replication of meteorological measurements and snowmelt 
collection measurements is necessary for adequate evaluation of heat 
energy flexes during special snowmelt events. Quality measurements 
of net radiation, solar radiation, precipitation, air temperature, rela- 
tive humidity, windspeed, snow profile temperature, and snow density 
measurements along with simultaneous measurements of water leaving 
the bottom of the snowpack are all absolutely necessary for developing 
and testing a meaningful snowmelt model. 


C) 
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CRIS Work Unit No. SWC-O11-fBo-1 Code No. Ida-Bo-102.1 


Title: Precipitation characteristics of a northern, 
mountainous, semiarid watershed. 


Location: Northwest Watershed Research Center, Boise, Idaho. 
Personnel: Panis Smith! and W. R. Hamon 
Date of Initiation: June 1961 


Expected Termination Date: Originally Planned: December, 1969 
Recommended: December, 1971 


Objectives: 


1. To develop methods for evaluating precipitation rates and amounts for 
watersheds of different sizes. 


2. To determine seasonal distribution of precipitation with respect to 
amounts, character, and areal extent. 


3. To develop depth-duration-frequency and depth-area-duration relationships 
for the Reynolds Creek Experimental Watershed. 


4. To establish general precipitation-elevation-aspect-slope relationships 
from precipitation data obtained in the Reynolds Creek Experimental 
Watershed. 


Need for Study: 


No dense recording rain gage networks exist in the Northwest. Such a network 

is necessary to delineate thunderstorms andsiwurmvariability. Furthermore, 
Weather Bureau data collection stations are generally located in or near the main 
cities. These cities are generally along the main stems of major streams cr 

in their valleys. Thus a sample of the range watershed areas is not available 
from their records. In addition to this, there is a dearth of rain gages 

capable of recording intensities or even individual storm data. 


The Soil Conservation Service is constantly called on to make estimates 


of precipitation away from the large population centers. Very little data 
is available to them on which to base their estimates. 





1/ 
~ On leave without pay during reporting period. 
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Design of Experiment and Procedure to be Followed: 


Data Procurement: 


1, A recording rain-gage network with a basic density of one gage 
per square mile will’ furnish the basic data. At selected locations 
additional gages will be added to furnish data to obtain 
refinement of the spatial variability of precipitation. 


2. Data on elevation, aspect, and slope will be obtained from 
topographic maps. 


3. Comparable precipitation data will be obtained from U. S. Weather 
Bureau gtations in the general region. 


Procedure: 


Precipitation-elevation-slope-aspect relationships will be 
established and basic data reduced to a common datum plane for 
storm analysis to obtain objectives. 


Experimental Data and Observations: 


Precipitation data collection was continued from the Reynolds Creek and Cc, 
Rabbit Creck Wetwork, Figures 1 and 2, pages 1-6 and 1-7, respectively. Data* 
on the precipitation gage networks and other instrumentation is noted in 

Table 1, page 1-5. 


Each precipitation site is instrumented with dual gages--one unshielded and 
one shielded gage. A mathematical model that uses the shielded and unshielded 
gage catches to computed "actual" precipitation, independent of wind, 

has been presented in previous Annual Report. 


An isohyetal map of Reynolds Creek Expermental Watershed for precipitation 
caught es unshielded gages, during the 1969-70 water year is presented as 
Figure 12/ Isohyets of the shielded gage catch in excess of the unshielded 
gage catch is shown in Figure 2. Computed precipitation totals, obtained 
from totals of shielded and unshielded gage catches, for the 1969-70 water 
year have been used to prepare an isohyctal map,-: Figure 3. Values obtained 
for this latter isohyetal map of computed total precipitation are up to 

4 percent too small in the higher precipitation zones since yearly rather 
than storm or monthly totals were used in the calculations. 


; 


2 
Zi Figures follow page 3-4, 


€, 


eee 
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Water year totals of unshielded and computed precipitation, along with 


runoff and evapotranspiration estimates, are tabulated for 4&4 watersheds in 
Table 1, 


All previously obtained precipitation data are to be converted to better 
estimates of actual precipitation by applying procedures developed under 
Research Outling Ida-Bo-102.6. 


TABLE 1. --Watershed precipitation, runoff, and evapotranspiration. 
eee 
Sine ape erer peng spn ap oe eects ae SDI aii els aaah anaemia 


/ 
Watershed Precipitation RunOTS nT, P.E.T.2! U/A 
and (U) (A) (R) (A-R) 
PAaCVAC LON hance, PFC, fee —) = eet tee Sate ew ot Inches- - ----+---+ - (%) 
ee 
Summit 9.25 10, 00 0 10.00 ah: 93 
(4180-4800) 

Lower Sheep eM Matic Be 14.88 0.03 EE: oe Mahe 78 
(5200-5428) 

Salmon wee =~ 18.85 CRE iene Be A 86 
(3675-6200) 

Tollgate Pe ele Py 30705 10, 14 20,51 29 78 


(4600-7300) 


dy: A computed 
The symbol U represents unshielded precipitation gage catch and d 
Sageressch. Oe ci yiration. 

2/ 


P.E.T. represents potential evapotranspiration. Value extracted from 1968 
Annual Report. 





Comments, Interpretations, and Future Plans: 


The computed water year precipitation, Figure 3, for the Reynolds Creek 
Watershed exceeds the precipitation caught by unshielded gages by about 1 
inch where the total precipitation caught by the unshielded gage was only 

10 inches, but exceeds the unshielded gage catch by nearly 14 inches at 
higher elevations where the total unshiclded gage catch was about 31 inches. 
Expressed as percentages, the observed unshiclded precipitation gage catches 
were about 89 percent of the computed values where totals were about 10 
inches, but only 69 percent at the highest clevations where the total 
unshielded catches were near 31 inches. The greater percentage of 
precipitation occurring as snow and the higher wind speed at the higher 
elevation contribute to the resulting smaller percentage of actual precipitation 
being caught by the unshielded gage. 
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; Water balance studies and runoff models require reliable input data if - 
results arc to have any applicability to areas outside the study area. In 
the case of Reynolds Creek Watershed, inadequate precipitation data would 
make it impossible to develop and test predictive hydrologic models. The 
importance of realistic precipitation data is demonstrated by the data on 
watershed precipitation, runoff, and evapotranspiration prescntcdsin Table 1. 
The evapotranspiration estimates in this table are obtained by use of the 
water balance equation with storage changes neglected. 


As demonstrated in Table 1 for 4 watersheds within the Reynolds Creek basin, 

the unshicld precipitation gage catches (U) are completely unrealistic. 

Use of these values would indicate an evapotranspiration lost for the water 

year of only 13.63 inches for the Tollgatc watershed as compared to an estimated 
potential cvapotranspiration of 29 inches. In comparison, the computed 
precipitation value for the watershed is G7 percent of the potential 
evapotranspiration. On the basis of computed precipitation and runoff data, 
Table 1, essentially no water yield occurred in the arca recciving less 

than 15 inches of precipitation for the water year or for elevations less 

than 5200 feet. 


Refined procedures for computed actual precipitation, as developed under 
Research Outline Ida-Bo-102.6, will be utilized to obtain realistic values 
of precipitation to be used in the required analyses to reach the objectives 
of this Outline. 


A paper 3/ has been prepared on the Reynolds Creek precipitation gage we ee 
It contains information on the dual gage procedures and on tentative 
precipitation-elevation relationships. 


———————— 


= Hamon, W. R. The Reynolds Creeks Precipitation Gage Network in 
Southwestern Idaho. Approved by Division for publication in, ARS 
Precipitation Facilities and Related Studies, ARS 41-176 (Edited by 
D. M. Heorshficid). 
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Precipitation 
Computed Totals 
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Figure 3° 
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CRIS Work Unit No. SWC-Oll-fBo-l Code No, Ida-Bo-102. 6 
Titles Evaluation of precipitation gage performance, 
Location: Northwest Watershed Research Center, 

Boise, Idaho 
Personnel: Bey My ae! and W.R. Hammon. 2 
Date of Initiation: December 1964, 


Expected Termination 
Date: Originally Planned: July 1, 1969. 





Recommended: December 1971. 
Objectives: 
1. To determine which of several commercially available pre- 


cipitation-measuring devices most accurately measure the 
actual amount of precipitation, 


f\D 
e 


To develop criteria and specifications for installing precipita- 
tion gages to provide the most accurate measurements of 
precipitation, 


3. To develop methods and procedures for adjusting records 
from precipitation gages with ''standard'' installations to more 
accurately reflect the point precipitation. 


Need for Study: 


Man has always been interested in precipitation because of his depend- 
ence upon it to sustain his way of life. First efforts to systematically 
measure rain are lost in antiquity, but they probably coincided with the 
early development of the art of making clay bowls and vessels. Writ- 
ings from as early as 400 B.C. describe containers with standard 
dimensions being used to collect rainfall for forecasting agricultural 
conditions, 


The art of measuring precipitation has changed little in the past 

2500 years. Modern-day gages are still simple containers with well- 
defined dimensions in which the catch is measured by reference to a 
standard volumetric container, use of a dip stick, or by weighing. 
Accepted practice is to place gages, when possible, at locations that are 


We ne Ean 


1/ On leave without pay during reporting period. 
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protected from the direct sweep of the wind and are no closer to 
buildings, trees, fences, etc. than three times the height of these 7 
objects. 


In intensive hydrologic investigations, it is important to know as pre- 
cisely as possible the input of precipitation to small watersheds or study 
areas, This requires two things: (1) a gage network of sufficient den- 
Sity to adequately sample the aerial variability of the precipitation, and 
(2) devices that will accuarely measure the precipitation at a point, The 
first of these, network density, can be determined by standard statistical 
Sampling procedures, For the second, however, it is questionable whether 
or not modern-day precipitation gages will provide accurate, quantitative 
point measurements. To resolve this question, an evaluation of the per- 
formance of precipitation gages under varying conditions of topography, 
exposure, and wind is required. It is only after such an evaluation is 
made that data from either individual gages or networks can be used with 
confidence. 


Design of Experiment and Procedures to be Followed: 


Evidence to date indicates questionable accuracy of precipitation data 
from gages at windy exposures, Experiments have demonstrated gross 
inaccuracies in snow catch under windy environments, In this study an e 
attempt will be made to determine: (1) Whether or not there is a ciffer- 
ence between the amount of precipitation caught in the various gages and 
the assumed "true" catch of a pit gage; (2) if there is a difference, what 
factors of wind, installation technique, topography, orientation, etc. 

cause the difference; (3) the criteria for assessing the site selection, type 
of gage, and installation techniques under various physiographic conditions 
to provide more accurate measurement of the amount of precipitation that 
reaches the ground surface; and (4) procedures and techniques based on 
site characteristics to correct precipitation records. 


Sites will be selected on the Reynolds Creek Experimental Watershed that 
are representative of the major topographic-elevation complexes. At each 
of these study sites precipitation measurements will be made by using an 
unshielded gage and shielded gage with orifices at the same height and 
exposed to horizontal wind movement. The same installation will be made 
at specicl evaluation sites with the following additional installations: (1) 
unshielded gages in a vertical profile, (2) snow pillows, (3) sunken pre- 
cipitation gage, (4) vectopluviometer, and (5) anemometers., 


Data obtained from profile measurements of precipitation and wind will 
be used to compute the actual vertical flux of precipitation. These com- 
puted values for windy, snowy conditions and data from snow pillows and ¢ 
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a vectopluviometer will be used to calibrate the relative catches of 
shielded and unshielded gages to obtain estimates of actual precipitation. 
These generated data will be used in a relationship between catch, wind, 
and temperature for correcting the catch of a single gage to account for 
the wind effect. 


Experimental Data and Observations: 


Data collection continued at the precipitation gage evaluation site on 
Reynolds Mountain (176 X 14), at the snow pillow site (176 X 07) about 
one-fourth mile to the east, and at the dual-precipitation gage site 
(166 X 94) about one-fourth mile to the north. The types of data col- 
lected at the 3 sites are indicated in Table l. 


The tipping-bucket, heated, gages were added to the experimental site 
during the fall of 1969, but little useful data were obtained during the 
snow season of 1969-70 because of heater blowout or failures, 
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TABLE 1,--Data collected at sites on Reynolds Mountain. for evaluatiogy 
of precipitation gaze performance, — ; 








Types of ake and att OV 
Height of Sensor, Ft, sites— 
. 176 X 14 176-507 166 X 94 
LS ere pa) x 
U (5) x 
Ueki x 
U (10) x bg eS 
U (14) zz 
S (5) ze 
3 (10) “a Ze xX 
a (5) oe ¢ 
lame ee) PS 
vf (10) Zn v4 
UL, (10) x 
oy (10) 4 
Sve x 
sc x 


1/ Instrumentation: U is recording unshielded gage; S is recording 
"  ghielded gage; U, is tipping-bucket, heated, unshielded gage; 
S, is tipping-bucket, heated, shielded gage; W is recording 
cnemometer with wind vane; Sp is snow pillow, and SC is snow 
course, 


2/ Instrument sites can be found in Figure 1, page 1-6. (Also see 
Yootnote 1/ on page 1-4.) ¢ 











Comments, Interpretations, and Future Plans: 


Sufficient deta will have been obtained at the clos of the 1970-71 snow 
season to complete the evaluation of precipitation gage performance 
and to establish firm values for the coefficients in the mathematical 
model to be used in calculating “actual” precipitation. This model, as 
developed in the 1968 Annual Report, is as follows: 





~ rare 
5 =e (1) 
= an QP (2) 
U _ _-(b-a)W 
ze (3) 
presi FBP ioe (4) 
= A eC we) 
B= a (4a) 
; U 
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where fi is the computed "actucl'’ precipitation; U is the unshielded 
gage catch; 5S is the shielded gage catch, a and b are coefficients, and 


o” 


Band C are termed calibration coefficients, 


The model requires data from the shielded and unshielded gages with 
orifices at the same height (high enough to escape blowing snow) and a 
rigid shield (constrained leaves at an angle of 30 degrees from the ver- 
tical) on the shielded gage. 


Profile dato on precipitation catch and profile windspeed, snow pillow 
and snow course data, and data from the unshielded and shielded gages 
(10 fect) at the 8 sites (Table 1) have been used to establish a tentative 
value of the calibration coefficient, B, or 


(6). 
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This calibration coefficient will be more accurately established when 
sufficient data have been obtained which should be achieved during the 
snow season of 1970-71, It has, however, been reasonably established 
that the calibration coefficient is not a function of type of precipitation. 


With an established value of the calibration coefficient, B, in Equation 4, 
the "actual’' precipitation con be computed by use of measured precipi- 
tation by on unshielded and shielded gage. (Equation 5 is equally applic- 
able since C is also known.) 


Precipitation data, however, are universally obtained by either an un- 
shielded or shielded gage except for the dual gage network that was 
established in the Reynolds Creek Experimental Watershed in late 1968, 
The earlier data in the watershed and data elsewhere need to be ad- 
jJusted for more adequate estimates of actual precipitation, particularly 
where snow is the precipitation type. 


Data for the period, November 1987 to December 1970, including data 

on windspeed, temperature, and catches by unshielded and shielded gages, 
have been analyzed to establish procedures for adjusting catches by either 
2 shielded gage or unshielded gage of the type used in the study. 


Data for 122 storms in excess of 0.1 inch were grouped according to the ¢ 
following 4 temperature ranges: 


with plotting made of values for U/S and for, wind ("7) on semilogarithmic 
paper, according to Equation 5, in Figures—'1 and 2, €urves drawn through 
the origin and mean are shown in the figures and grouped in the bottom 
groph of Figure 2, From these curves, the value of the coefficient, 

a, was established for each temperature range as shown in Table 2, 
With o reasonable estimate of the calibration coefficient, B, known 
(Equation 4), values of the coefficient b are established, Equations, 

1 and 2 can now be used to compute actual precipitation, given measured 
precipitation by either an unshielded or a shielded gage and data on wind 
and temperature, It is significant to note that with B = 2.0, (b—w)= 
1,0a, and the ratio U/S is replaceable by S/A in Figures 1 and 2. 


¢ 





2/ Figures follow Page 4-8, 
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TABLE 2,--Coefficients to correct for the effect of wind and temper- 
ature on precipitation caught by shielded and unshielded 
precipitation gages. 











b 
For B= =a = 2,0 
Coefficients en nor 
Lemperature, “i, 
pes i Se CeO eee eo) LS eS 
2 0,070 UL 0.030 0,053 
b 0, 014 0, 054 0, 060 0.106 


Note: bawz- (2,0 X -) 





The derived coefficients are usec in Equations 1 and 2 to obtain values 
of the ratios of shielded and unshielded gage catches to actual precipita- 
s tion for two selected windspeedsand the results noted in Table 3, 


The research task remaining is to establish a reliable value of the 
calibration coefficient and to test its variability with different types of 
dual gage installations. Further tests are to be made on the use of 
heated gages to eliminate the capping and riming of gages that invalidates 
the use of the collected data in ony fashion. 
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TABLE 3,--Values of the ratio of unshielded, U, and shielded, S, 
gage catch to computed "actual" precipitation, A, for 
ranges in temperature and for two windspeeds, 


See EEE—EE———EE_EEEE eee 
Ratio and 
Windspeed Temperature, °F, 


timsss® 82° =< t < 35°. 32° eater t S250 
WOACOGt Colh, 


OE I ae Pe Na ee ee ar ee ee at ee ey vee ioe oes ee 


W * 30 m.p.h, 
U/A 0.66 0. 36 0.17 


S/A 0, 81 0,60 0, 41 0.20 


Z. 
A Ee ee been prepared outlining the aspects of the dual gage and 
profile technigues for calculating actual precipitation, 


47 Haimon, W.R. Dual gage and profile techniques for calculating 
actual Precipitation, Working Group on Measurement of ces oa 
tion, Commission for Instruments and Methods of Observations, 

world Meteorological Organization, Gensra, Switzerland, October 1970, 
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Figure 1. Ratio of unshielded (U) to shielded (S) gage catches 
in relation to wind and temperature. 


TT pan ES OE IT Pa eee ST PE I Ee 








1.0 
0.9 € 


0.8 9°< T < 23°F, | 


0.6 | : 


U/S 


0.4 : 


os SITE 


(x) 176X007 
(6) 176X114 


ae WIND, M.P.H. 


S20 TP Ss5ten 


U/S 


Soa 1 > 23°F. 


0.3 


0,2 9°<T < 23°F, 





5 10 1S 20 25 30 
WIND,M.P.H. f € | 


Figure 2. Ratio of unshiclded (U) to shielded (S) gage catches 
in relation to wind and temperature. 
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CRIS Work Unit No. SWC-Ol-fBo-l Code No, Ida-Bo-102.7 


Title: Evaluation of pressure pillows, and hydraulic 
weighing and catchment devices for snow and 
precipitation measurements, 


Location: Northwest //atershed Research Center, Boise, 
Idaho. 
Personnel: L. M. Cox and W. R, Hamon; M. W. Nelson 


and J. A. Wilson, SCS. 


Date of Initiation: peptember 1967 


Expected Termina- 


tion: Originally planned, August 1969, 





Recommendation, July 1971. 


Objective: To determine accuracy and utility of steel 


diaphragms, butyl or plastic pressure pillows, 
and a combination hydraulic-weighing and :catch- 
ment device, for use in 5CS snow surveys and 
ARS hydrologic research studies, 


Need for Study: 


The Soil Conservation Service is currently using many types of pressure- 
sensing devices to measure snow-water content. Measurements obtained 
by 12-foot-diameter pressure pillows have been superior to those ob- 
tained by snow-tube sampling, These pillows made of butyl or plastic, 
however, are subject to leakage and physical rupture. The initial 

cost is excessive, 


Field tests are necessary to determine the accuracy and utility of new 
types of snow-water-measuring cevices. Data obtained by such devices 
need to be correlated with snow-tube measurements to make past rec- 
ords fully useful in the preparation of water supply forecasts, 


Besides the measurement of snow-water content, a universal device is 
needed that will also measure the total precipitation and snowmelt. The 
availability of snow-water content and total precipitation data each day 
is essential for the prediction of peak flows, low flows, and monthly 
increments of flow. Data obtained by such devices are required for 
hydrologic-water balance studies, 


at cn. 
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Design of Experiment and Procedure to be followed: 
eee OG tOMlOWOCC: 


1, Commercially available 12-foot diameter butyl or neoprene Pillows 
will be used as standards for comparisons, Steel diaphragm 
pillows will be constructed to SCS Standards, A universal gage 
to measure snow-water content, total precipitation, and snowmelt 
will be an ARS-designed instrument, Installations will be as noted 
in Table 1 below: 


TABLE 1.-- Snow pillow installations 
dentate sie siroonn selectins med oe 


| hee ee ec 


Locatians 
Reynolds Bogus Atlanta Trinity 
Min, Basin Summit Mtn. 
ee ee 
flev., in feet (m.s.1.) 6, 800 6,120 7,500 1, 180 
Av. Snow depth in feet 3 6 8 10 


(Sst. ) ° j 


Butyl/plastie pillow i/ 2/ 
(1Z-ft. diameter) x = xX 


Steel pillow 
io awOv Ott) Po 


xt 
. 


Steel pillow 
(4-4 ft. by 5 ft.) 
(srouped) * * 


% 


Universal gage 


So 
4 


(5-ft. diameter) x 


oh 
7 


i/ New installation, 
2/ “Stablished, 


2, All instruments will be celibrated with the 12-foot pressure pillow 
and independent calibrations obtained by gravimetric measurement 
of snow-water content, Samples will be obtained by snow tube 
and by other feasible methods. o 








’ : 
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Experimental Data and Observations: 


Snowfall for the 1969-70 snow season was quite heavy with accumula- 
tions at the higher elevations approaching those of the previous year. 
Travel to the Atlanta Summit and Trinity Mountain sites was very 
difficult. come snowmelt data were collected at Trinity Mountain but 
hazardous travel conditions prevented regular servicing of the record- 
ing equipment. 


The Soil Conservation Service is in the process of installing a new 
telemetering system at the Trinity Mountain study site. This new 
system has a greater range plus the capability of handling more inputs 
from different sensors. The SCS plans to install new equipment to 
telemeter the water equivalent measurements frorm two 8 by 10-foot 
steel pillows and the universal gage. Once the system is working, a 
transducer will be added to telemeter the data from the melt collector 
of the universal gage, The 1z-foot butyl pillow at this site developed 
an unrepairable leak, and it was replaced with another 8 by 10-foot 
steel pillow, The 12-foot butyl pillow at the Bogus Basin site has also 
been replaced with an 8 by 10-foot steel pillow. 


Correlation coefficients for pressure-pillow response and snow -course 
measurement are given in Table 2. At both the Reynolds Mountain 
and Bogus Basin study sites, all measurements were equally corre- 
lated for the total snow season. In the past it was found (last year's 
Annual Report) that all measurements usually correlate very well dur- 
ing the accumulation period, but show a distinct lower correlation dur- 
ing the melt period. The 12-foot butyl pillow performs somewhat 
differently during the melt period as compared to the universal. gage, 
This difference in performance appears to be closely related to the 
rate of snowmelt with low snowmelt rates givins good correlations 
ana hich snowmelt rates reduced correlations. | 


The longer the snowmelt is delayed in the spring by cool, cloudy 
weather, the greater is the potential for higher snowmelt rates, As 
the spring season progresses, more energy becomes potentially avail- 
able to melt the snow because of the occurrence of higher air temper- 
atures and subsequently greater vapor pressure eradients thus increas- 
ing convection and condensation melt; higher sun angles for greater 
radiation melt; and a greater opportunity for the acvection of heat 
energy from bare soil surfaces to the existing snowpack. The-per- 
formance of snow pillows.during 2 period of late season melt is 
demonstrated in Figure 12! This graph shows the changes in water 
equivalent of the snowpack, with time as measured with a 12-foot butyl 
snow pillow (top curve) and a universal gage (middle curve). The 





1/ Figure 1 follows page p. 5-6, 
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TASLE 2, --Correlation coefficients for pressure-pillow wind Wt e 
and snow-course measurements, 
nonce ication ct Rinne, he SIRO RS NS ENR 


eno ereeainlieince nest leapakinteeshiaananin tiie aeietneiainabex~ sinh uel Perms Le MO Uae wt, obs) ee ae 
—— 
Units Correlated Reynolds [fountain and 


Bogus Basin Sites Combined 
—_—$S $i ee ites ombined 


aes 


Univ. Gage by 12-foot butyl eed 


Snow Course 
Measurement by 12-foot butyl eee, 


Snow Course 
Measurement by Univ. Gage -o 


h ios 


eee — 


snowmelt rate as collected and mensured with a universal gage is 
shown oat the bottom of the graph. These data were collected at the 
Bogus Basin snow research site with about 51 inches of snow and 
23.0 inches of water equivalent. 


The changes in water equivalent of the snowpack as measured by the € y 
12-foot butyl pillow shows (with one exception) 2 continuous linear 
decrease with time for the four-day period. This exception was 
brought about by a rainstorm that occurred at 1400 to 1600 on May 19, 
1970. No rain gage is located at the site so no estimate of the total 
precipitation during this time period was possible, The pillow meas- 
urement shows a characteristics increase when the rain first fell on 
the snow and a sharp decrease os the water started coming out of the 
bottom of the pack as indicated by the melt collected, (A further dis- 
cussion of rain-on-snow can be found in Research Outline Ida-Bo-100.1). 
The lz-foot butyl pillow measurement does not show the decrease in 
melt per any as indicated by the daily decrense in snowmelt collected, 
Very little evidence of the diurnal melt pattern is reflected in the 
lé-foot butyl pillow measurement. 


The water equivalent as measured with the universal gage does demon- 
strate three distinct slopes for the water equivalent changes during the 
four-day period. Each particular slope corresponds to a specific melt 
period with decreasing slope being associated with decreasing melt, 

The water equivalent as measured with the universal gage follows the 
diurnal melt cycle much better than it follows the 12-foot butyl pillow. 
The heavy rainstorm that caused the disruption in the 12-foot butyl g 
pillow measured on May 19, 1970, was not as noticeable in the univer- 
sal gage measurement, 
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> From 1600 on May 17, 1970, to 0900 on May 20, 1970, the decrease 
| in water equivalent of the snowpack as measured with the universal 
gage was 2,75 inches as compared to 2.10 inches for the butyl pillow. 
The total amount of water collected from the bottom of the snowpack 
was 4.60 inches during this same period. This would indicate that 
at least (4.80 -2,75 = 2,15) 2,15 inches of water leaving the snowpack 
came from the rainstorms and condensate, Boise valley air temper- 
ature ranged from a low of 49°F. to a high of 65° ir, during the 
period. Up to 0,50 inch of rain was recorded in the valley, but 
heavier storms are usually the rule for this research area which is 
about 4400 feet higher in elevation. Measurements of precipitation, 
air temperature and relative hurnidity are very urgently needed at this 
Site for a more complete evaluation of the snowmelt process, 


Comments, Interpretations, and Future Plans: 


The feasibility of maintaining Atlanta Summit and Trinity Mountain 
as research sites was recently questioned because of the increased 
travel time required to service those sites. Further discussion 
revealed that these are the only deep snow sites (8-10 feet) we have 
is and that Trinity Mountain is also one of the most important snow course 
9 = measurement sites, Data from this site are used to predict stream- 
flow intc the Boise River irrigation and flood control complex, Be- 
cause of its great importance for streamflow forecasting, the decision 
was made to maintain the Trinity Mountain site. 
With the installation of the new SCS telemetering system, the potential 
for better data collection from this site will be greatly increased, This 
new system will make it possible to telemeter all types of meteorol- 
ogical data from both Trinity Mountain and Atlanta Summit. Measure- 
ments of air temperature, relative humidity, windspeed, soil temper- 
ature, snow temperature and radiation will be made as soon as econ- 
omicclly possible. 


The 12-foot butyl pillows are beins gradually phased out. On the 
average, the life of these pillows is only about five years. They are 
being replaced by the cheaper and more versatile 8 by 10-foot steel 
pillows. 


The difference in the slopes of the water equivalent measurement 

(Figure 1) as made by the universal gage and the 12-foot butyl pillow 

clearly demonstrate why the two units have a lower correlation dur- 

| ing the melt period, During lone periods of high melt, water tem- 

) he porarily mounds on the surface of the 12-foot butyl pillows, thus giv- 
ing a false impression as to the time melt rate. In an old ripe 


ars 
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Snowpack, the universal gage does not exhibit this effect, because 
as water is added to the snow by rain or melting, it runs out through 
the collector at the bottom. Very little water is held in the old snow 
to ternporarily increase the water equivalent measurement. New, 
cold, low density snow would act to the contrary and under certain 
conditions, “much of the rain added to the snowpack may not run out 
the bottom. (See Research Outline Ida-Bo-100.1, ) 


Downstrear runoff predictec from water equivelent measurements from 
a 12-foot butyl pillow would obviously be in error during periods of 
high continuous melt, Many forecasters in the “/estern United States 
are currently using telemetered water equivalent data from 12-foot 
butyl pillows. The true measure of what is happening is indicated by 
the melt collected from the bottom of the snowpack. This method also 
gives additional information on condensation water and rainwater, 


This study has brought further questions to mind concerning the pro- 
cedures tQ collect snowmelt. This research outline will be termia 
nated following the 1970-71 snow season. A new outline will be ini- 
tiated and will involve a a study concerning whet size of melt col- 

lected is required for accurate measurements of snowmelt. Also ¢ 
more study is needed on the rate at which melt water can be ve tioue cee 
from 2 snowpack and still mcintain a representative sample of current 
melt conditions. Under high melt conditions, if the water is removed 
too fast, then a greater area may be drained than the area of the 
eollecton; : 
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Figure J. Changes in water-equivalent with time as measured with a 12-ff, butyl snow pillow and a universal 
‘gage during late season snowmelt. Indicated snowmelt was collected from the 57-inch snowpack 


Using the melt catehmand cnt nas HF th. ty 








24 


IIZWMONS 





ap en Sop ee 


we Oy 


9 


macnt ee thin 
Aenrefs ape) ads bist ~~ 7 ; 


ans ae i e; - oT 


XV 3 Wwather ie: prone to the snow by: Ad OF 2 ane ‘ing, ¥ + rune 
oO the soltupian int av the splaagions Very Uitie is & F- tat 
vay 




















K 





= <n. j aiouk os 
e bo “BOF pore rily ae ise the water aguay eh cy Divan 7 ee 
ae Ya : A r j : . 
= 2 Se) sol} “Low Mons! ty snow would act) to the: contrary cunts 
Sit i. \s ; find! 
5 4 =) copitions, mauch of the rain addéd¢ to the snownies ancy met 
ws @ ; a me.* ire ae fda. 2 an ote * ne, 
~ Eg - ! “Se bottern, (oe Hesearcn utline Ide-Bo-100.L,? ' 
. a | ©) i PP =z 
> 8 ” — Ze © 
— f 
4 a . > Tene 
“4 ron ; Downstr: . wf redicted [6 
“ a +. — 
pF r + _ £ 


a 
water équivy aie : ee 
ert 
: 


E = . : Lt \ ! o t . pai N ok, wc’ } di¢ Lat is) ‘ } ein erro aided i ,p oe ae 
“ae Nee . . ray Sa oe v 
hone _ \), (. ’ ) 77, ow aay Gi x ae a . eto rs in. he o> ay 1 % ed 7 


| i 7 
3¢ mi. | ager nt) ing (Alemetere@ water equivalent py Gog dees 
“ar: 4 : bom f 7 a 
> ‘yl pillows e true ae ¥ happenin ra ins ote tk 


ca ; Aa J ‘ ‘ ‘ ; | ‘ 
= ~~ yb & . ; i+ ¢ Slice ‘ Ys from, : ‘ } Bor: ou ier SnnwaniGk, The ‘phal cag md 4 
a sd | ] 4 ‘ as i swe i J 
o- = » 7 : a” ii | ? 4 : , a « 
=o : - + 
: Viti, i. 
7 i. 
i ‘ 
{> 
i iT 
: rc . 
Pe - 
' a a 
2 y= me 
j 
od 
ms : 
é oll . 
=. a ) 


4 
= 
, 7 
7 
7 


. elie 
\ ; a 
, 
¢ ‘ \ q rs 
~~ - re = ;. 7 
- é ‘ 
> kl . - 
“ q 
| 7 
a ~ 
i , a j 
i - a 
q 
”~ 2 \. a 
; 7 
r . ; 
4 ‘ 
- = 
P Le fa *. 
¢: , 
® a 
; 7 < 
a one y 
a zs 
»’ e 
. : ~- 4 
. a2 et aa” 
4 : - 
. ¢ ke» f 
> <*s 5 
3 : & oD ‘ 
« = > | } 
al 
- $e ; = 
« - d 
~ & he : * 
~~ & * 
edt v ) 
- oo» = f 
= wa a] j 
~4  . “~ 
ys 3 , 


4 


ex 
rs 
re 
ow 
3 
— 
— 









rivA ti" fi 


a & +4 
| seed onan : 
rs 





. 


© : 7 


' 





6-1 


CRIS Work Unit No, SWC-Oll-fBo-1 Code No, Ida-Bo-104,1 
Title: Resistance coefficients for steep-rough channels, 
Location: Northwest «/atershed Research Center, 


Boise, Idaho. 


a 


Personnel: C. W. Johnson, and W. R. Hamon, 
DBD, E, Overton (USDA Hydrograph Laboratory); and 
Dr. HH. =, Judd (cooperating through Utah State 
University) 


Date of Initiation: March 1967 
Expected Termina- 
tion Date: ‘Originally Planned: December 1970 
Recommended Termination: July 19% 
Objectives: 
1, To test the applicability of classical flow resistance equations 


to flow data from natural channels in the Reynolds Creek Ex- 
perimental Watershed, 


2, To characterize channel geometry and bed roughness 
elements for use in investigating resistance coefficients, 
Ss. To determine the relationship of the friction factor or 


roughness coefficient to channel characteristics. 


Need for Study: 


The prediction of the resistance to flow of natural streams is vital to 
the ultimate goals of mathematical simulation of the hydrologic per- 
formance of upland watersheds. The successful routing of runoff 
through a complex watershed stream-system is directly dependent upon 
an adequate prediction of channel resistance. Field research is needed 
to develop quantitative procedures for the representation of resistance 
coefficients; i, e., friction factors or roughness coefficients from easily- 
obtained field data on channel dimensions, and on size and spacing of 
roughness elements. Field engineers need an objective procedure for 
predicting depth of flow associated with the design discharge in the 
design of water-conveyance channels and streambank protection works. 


‘ 


Design of Experiment and Procedure to be Followed: 
eee et BE PP oLowed 


Cooperation of the USDA Hydrograph Laboratory, Utah State University, 
and the Northwest Watershed Research Center is continuing until analysis 
of available data from laboratory flumes and natural channels is com- 
plete, A thorough understanding of flow resistance in laboratory flumes 
is a necessary prerequisite to analysis of roughness in natural channels, 
A study of the effect of laboratory channel roughness parameters on 

flow resistance is progressing through the use of computers and results 
should be helpful in analyzing the more complex problem of roughness 
measurement in natural channels. 


Data on channel geometry, stream discharge, size and spacing of rough- 
ness elements and slope of the water surface or channel are used to 
determine channel resistance, The Chezy eguation 


and Larcy-Weisbach equation 


Q=A V8gRS 


are used to determine resistance or roughness coefficients C and ¢ 

where © is discharge, A' is cross-sertional area of the channel, R 
is the hydraulic radius (both A’ and Rdenote inclusion of a roughness 
parameter), S is the energy slope (bedslope for uniform flow), and g 


is the acceleration of gravity. 


t 


Channel reaches on Reynolds Creek Experimental Watershed streams 
are instrumented to provide data from natural streams and existing 
structures are utilized where possible, 


Gary 


ixperimental Data and Observations: 
een ee ee Eee ee tet: 


Analysis of a great bulk of data from laboratory flumes and natural 
channels was completed during the reporting year and no additional 

field data were collected. The determination of relationships between 
the resistance coefficient and parameters such as effective depth, rough- 
ness intensity, channel width and mean size of bed material Was not 
equally applicable to laboratory flume studies and natural stream chan- 
nels with large bed elements. 
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Comments, Interpretations, and Future Plans: 


A procedure was developed for locating the plane of zero velocity near 
the channel bottom of Inboratory flumes and natural streams with large 
bed elements. Also, an equivalent spacing for randomly-spaced bed 
elements was developed by use of regularly-spaced patterns of hemis- 
pherical bed elements in a laborctory flume. The random bed elements 
produced greater roughness coefficients than regular patterns of equal 
roughness intensity and increased the energy losses from 34-61 percent 
of that from regular patterns tested. Further, a study of wall ‘effect 
in leboratory flumes of various widths showed that narrow channels 
produce greater resistance to flow than wide channels of equal rough- 
ness intensity and cross-sectional area. Thus, a greater understand- 
ing of the influence of random bed elements on flow resistance in 
natural channels has been achieved, 


A draft of a report, ‘Constant Resistance Coefficients for Large Bed 
Stream Elements, '' is in the process of review for publication, Termi- 
nation of this outline will follow approval of the report for publication, 
Further work will require a new outline and objectives, | 

Work covering the initial investigations under this outline has been 
presented (ASAK National Water Resources Engineering Meeting, Mem- 
phis, Tenn,, January 26-30, 1970) and an approved manuscript is pend- 
ing publication, 


Manuscripts: 

Prepared for Review: 
Cverton, 0,E., Judd, Hari &., and Johnson, ©C.W. Constant 
resistance coefficients for large bed element streams. Proposed 
for publication through the Utah Water Research Laboratory, 
ijoenn, Utah, (1971, ) 

Pending rublication: 
Overton, D.E., and Hamon, W.R. A new resistance standard for 


rigid open channels. Under revision to meet acceptance for pub- 
lication in ASCE, Jour. of the Hydraulics Sivision, (1971. ) 
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CRIS Work Unit No.: SWC-Oll-fBo-1 Code No. Ida-Bo-105,4 

Title: Evaluation of cover production, herbage yield, 
and soil conditions for different levels of man- 
agement, 

Location: Northwest v/atershed Research Center, Boise, 
Idaho, 

Cooperation: The Bureau of Land Management will provide 


partial funds and furnish public lands on which 
to conduct most of the investigations. The 
Soil Conservation Service will assist in obtain- 
ing soils data and cooperate on other aspects 
of the study relative to their mission, 


Personnel: G.A. Schumaker and J. F. Zuzel. 
Date of Initiation: February 1971 


Expected Termina- 


tions ace: Field work: December 1974 
Interpretations and Summary: December 1975, 
Objectives: 


l. To determine the effect of grazing management and treatments on 
yield of herbage, cover production, soil moisture regimes, and 
soil surface conditions at selected sites, 


2. To study changes in plant density and plant composition as a 
result of grazing management and treatments. 


Need for Study, and Literature Review: 


Quantitative data on herbage yield from rangelands under different 
levels of management are needed to guide managers of public lands 

in optimizing multiple use of the range, Such data are of equal im- 
portance to the Soil Conservation Service, other public bodies, and 
landowners. These needs require more discerning information on how 
vegetation and soils respond to imposed treatments, including con- 
trolled grazing. Grassland management and production factors are 
intimately related to other facets of the Reynolds Creek Experimental 
Watershed research program; particularly water use, runoff, and sedi- 
ment production, Information is also needed with regard to methods 
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of increasing cover and to the rate of recovery of native range 
following intensive grazing practices. This information is needed 
to gain insights on the total range utilization process; that is, a 
stable cover that gives optimum grass production in relation to 
water yield with a minimum of sediment loss. - 


The $0-square-mile Reynolds Creek Experimental watershed, rang- 
ing in elevation from 3600 to 7000 feet, m.s.l., offers a variety 
of soils and climatic zones for study. Vegetative cover is exceed- 
ingly variable in the Watershed with cover-area percentages as 
follows: 


Percent Cover rercent of Area 
0-25 64 
26-50 10 
51-75 12 
76-100 14 


dh rhesus oh 


Moderate sediment yields occur from these sparsely vegetated areas, 


Such sites offer an opportunity for optimizing sediment reduction 
through management control. In the more densely vegetated areas, 
in the higher precipitation zones, herbage yields can be optimized 
while controlling erosion. 


The need to learn more about vegetation influences on watersheds 
and to develop ,watershed-oriented management has been pointed out 
by Dunford (1)— Cover and cover management play an important 


role in watershed hydrology and the effects of cover have been inves- 


tigated by various research workers. 


The effects of grazing and subsequent cover on the hydrology of 
salt-desert type rangeland was studied by Lusby (3) at the Badger 
Wash near Grand Junction, Colorado. Cover measurements taken 
ten years apart showed an incrense in bare soil and exposed rock 
and a decrease in ground cover where grazing continued during the 


study period, lLDifferences in runoff between the grazed and ungrazed 


treatments were related directly to percentage of bare soil on the 


watersheds under study. The greatest improvement occurred ‘within 


the first three years after cattle had been excluded. 


The results of Meeuwig (4) obtained from rainfall infiltrometer, data 


emphasize the importance of vegetation and cover in raaintaining in- 


filtration capacity and soil stability on high elevation herbland on th 


Wasatch Sront in northern Utah. Cther work reported by Mee uwipnts 


shows that certain soil properties can be used as predictors of 


1/ Numerals in parentheses refer to corresponding items in the 
Literature Cited on p. 7-10. 
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infiltration although infiltration was not highly correlated with any 
one soil property. He was able to identify those site factors which 
are important in precipitation detention and soil erodibility. In ad- 
dition, a study by Dunin (2) shows oa greater infiltration and higher 
moisture storage from the more productive catchments near Bacchas 
Marsh, Victoria, 


Studies concerning rangeland vegetation in the Reynolds Creek Experi- 
mental «/atershed are essential to determine the effects of vegetation 
upon the hydrologic performance of watersheds. Mcet of these studies 
are to be made on watersheds where data will also be collected on 
infiltration rates, runoff and sediment yield. The watersheds sample 
the range of climate, geology, and soil on the 90-square-mile water- 
shed. The watersheds range in size from 30-500 acres. Small to- 
talizing watersheds in operation include the Northeast Summit, Sheep 
Creek, Reynolds Mountain, and Nancy Gulch, 


Design of =xperiment and Proceduresto Be Followed: 


A. Variables under Investigation 
l. wcological Zones 


a. Species composition and cover 
Db. Elevation 

c. Precipitation 

d. Physical characteristics of site 


(1) Soil type 
(2) Slope 
(3) Depth of topsoil or erosion 


Seasonal Rainfall 
Cover Management 
Soil Moisture 
snow Cover 


cn Won 


B, Variables not under Investigation 


i. rlant Disease 
Trenatments: 
A. Y©oor cover sites 


l. Exclosure of one acre or more 
2. Grazed 
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Good cover sites 3 ~ 


1, &xclosure of one acre or more with three types of 
brush treatment 


™ “Brush removed 
A Brush burned 
- Brush remaining 


Qe & 


ie} 


a Grazed 


Development of Experimental Sites 


site Selection 


1, Six to eight study sites will be selected on shallow soils 
where poor cover exists. Sites selected will be of vary- 
ing soil types and within varying rainfall zones: therefore, 
sites will be different ecologically, 


a, ‘two sites will be selected on moderately deep soils where 
© good grass cover exisis. 


o. Two sites will be selected on deep soils where a dense 
cover of sagebrush exists in the higher precipitation and 
snow accumulation zones. | 


Installation of Exclosures 


At each study site an aren of one ncre or more will be fenced 
for the purpose of excluding grazing animals. fn adjacent 
study aren will be selected which will be comprrable to the 
exclosure with respect to soil, slope, aspect, and cover, This 
latter site will be heavily grazed to maintain - poor cover site, 


Instellation of Soil Moicture Gites 


Soil moisture access tubes will be installed -t study sites for 
nonitoring soil moisture and to determine infiltration rates by 
use of a rainfall simulator. 


Vegetative Site Installations 
2 
1, At each site open to grnzing six areas of 12 ft,” will be 


protected from grazing with cages. New sites will be 
established before each grazing season, 
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ca) 


ee vWithin each treatment -t ench site two nine ft. areas 
will be marked for noting changes in species composition 
and cover, 


Description of Range Condition of Sites 
See 5 ain war seh fe St et ie si kn een he sd calico 
ae ae 3razing treatments are imposed on study sites, inventory 
data will be collected relative to erosion condition, ground cover, 


etc, 


Dat Collection Procedures 


A, LUry Matter Production 

1, Ory matter production will be measured twice each 
season, Clipping will be ai Eee when early season 
grasses have headed and again after late senson grasses 
have produced heads. 

4. All annuals will be removed, dried and weighed, That 
portion of perennials which represents annuel growth will 
be removed by clipping and weighed. 

3. © roduction will be Shera on individual species when 


feasible. 
B. Species Composition Sites 


1, lant species and size of grass clumps will be measured 
twice each season on the permanent reference sites denoted 
for measuring species composition. 


¢. 2och reference site will be photographed oat the time of 
data collection in order to more carefully observe changes 
over the period of. study. 


C,. ‘Soil moisture will be measured at weekly intervols during the 
growing season, Additioncel determinations will be made follow - 
ing sizeable storms, 


Dats to be Cbtained 


1, Yield of dry matter, 


4, «© eriodic inventory of plant species and cover, 


ieee 3 SD RR TEE rte nt een (Te a ecg ee 
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3.  90il1 moisture content. . te 


es 


. é/atershed inventory of physical characteristics for each 
study site as developed under BLM guidelines; both at 
initiation and conclusion of study, 


Experimental Data and Observations: 
nnn 


Six shallow soil sites were developed during the 1976 season, 

Fencing the exclosures was completed nt each of these snallow-soil 
study sites with the exception of Upper Sheep Creek and the Nettleton 
private holdinz, All instrument sites have been selected and marked. 
Installation of soil moisture access tubes was started in 1970 and was 
more than 50 percent completed. Instrumentation for soil moisture 
at the shallow soil site on Upper Sheep Creek and the Nettleton pri- 
vate holding was incompleted. Akout one-half of the Nancy Gulch soil 
moisture sites was completed, 


A deep soil site was selected on Reynolds Mountain 2t the 7000 feet - 
elevation in addition to the north slope, deep soil Upper Sheep Creek 
site where much of the exploratory work prior to preparation of the 
research outline was conducted. As outlined under procedures for 
conducting this research, the deep soil sites will receive two types 

of brush trectment along with the exclusion of grazing animals. Brush 
removal was completed at the Reynclds Mountain site and the other 
brush treatment will be imposed before growth begins in 1971, 


Additional soil moisture access tubes were installed at the Upper 
Sheep Creek site. These sites will permit soil moisture mexsure- 
ment to a depth of six feet. Permanent reference sites were also 
marked for the measurement of changes in species composition under 
the different treatments, 


Current data consist of that from exploratory work. The effects of 
the trentments initiated in 1969 on the -Upper Sheep Creek study area 
were very much in evidence during the 1970 season. On treatments 
where sagebrush had been removed send where the sagebrush had been 
sprayed, the native grasses displnyed vigor and were becoming estab- 
lished where brush cover had previously been in existence. Grass on 
the ungrazed treatment where sagebrush remained clso displayed some 
vigor owing to the absence of grazing animals. Flerbage yields 
taken in August show the production from the four different treatments, 
Table 1. 
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TABLE 1.-- List of species clipped from Upper Sheep Creek, August 1971, and weight of dry matter removed 


from cach treatment. 





=) = 
ee re — Brush Removed | No Treatment No Trdatmenee 
No Grazing ‘ No Grazing No Grazing Grazed 
Species Species Species Species 
and Total Yield and Total Yield and Total Yield and Total Yield 
Stipa sp. Stipa spp. Bromus marginatus Sitanion hystarex 
Carex sp. Sitanion hystrix Stipa’ spp. Poa spp. o 
Sitanion hystrix Poa spp. Sitanion hystrix Fescus idahoensus 
Poa sp» Bromus marginatus Agypuron trachycaulum Stipa 
Bromus marginatus Sedge Fescus idahoensus Foxrbs 
Annuals & Forbs Annuals & Forbs sage Herbaceous sage 


Rabbit brush Rabbit brushie 
Annuals & Forbs 


wee em meee eee ele 


-~-Weight of Dry Matter Removed-- 


773 lb./acre 793 1b./acre 267 1b./acre 





Lie : : : 
=" Big sagebrush and rabbit brush covered 25-50 percent of the area of sites measured. 


Species were identified, however, animals had grazed clipping sites prior to measurement, therefore, 
the grazing influence was not measurable. 
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Data were collected from the perraanent reference sites within ¢ 
each of the treatments, Plant species were identified and the basal 
cover of the grasses was measured. Where shrubs were present, 
measurements of aerial c over were also made. Other similar 
measurements have been taken from these sites annually during the 
study period, the progressive change in species composition and 

basal area can be observed. 


Soil moisture measurements were t-ken periodically with the neutron 
probe from the different regimes of management during the growing 
season, Table 2, 


The aspect of slope at the Upper Sheep Creex study site is an im- 
portant factor in snow accumulation, olso the type of cover can affect 
the depth of snow accumulation. Snow depth measurements were 
taken ecrly in 1970 at random points on each of the trhee treatments, 
Table 5S. 


Comments, Interpretations and Future Plans: 
Dit acetate Ie anal eae Aiel salpoak othe ede 5S ot Ik ee es Se 
A, Comments and Interpretations 


The excellent herbage yields obtnined from either of the 
brush treatments indicates that native range grasses have an 
ability to establish themselves once sage competition has been 
removed, 


Current soil moisture mensurements taken to 7 cepth of three 
feet indicate a slightly grecter cepletion of store moisture for 
the cepth of soil from the brush removed site, Further studies 
will indicate whether this trend continues. 


Snow depth measurements indicate that less snow accumulated 
where sagebrush had been removed, It is very doubtful that 
any differences were great enough to alter the weter balance on 
the brush removed treatment, 


B. Future Slans 


The development of other study sites in 1970 reached a stage 
where data collection can be initiated during the growing season 
Ob ISH: 
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TABLE 2.--Soil moisture measurements within exclosure at Upper Sheep 
Creek during 1970 growing seasons. 


A. Treatment Within Exclosure 


1. Brush Removed 


Date 4-23 5-1 5-6 5-15 5-18 6-16 7-22 8-5 8-19 

Depth 

Measured 
i So 3.4 4.6 aS u. 3 1.6 25m 1.8 1.8 
2 6.83 Pent, se pe oye Gok Ly eye 2.0 9 
3 Thad See fod 2 bye cree eae Lay 
a Fes 


Sie Cee 8 ee a ES eT a 
Date 4.-23 5-1 5-6 5-15 5-18 6-16 1222. ,0°o 8-19 
Depth 
Measured 
£ 6.3 5.0 Sy 2 Lad 4.8 2.0 B43 i a jena 
2 Oui 437 Sy) Sr 2 532 dade ar, 2.9 ere 
3 6.5 Ae 4.6 &.5 2.0 ae 3.0 2 
4 ei 3.5 


3, No Treatment 


Date 5-1 5-6 5-15 5-138 6-16 7-22 o-)5 8-19 

Depth 

Measured 
it 4.8 4.9 4.9 (2 Seats) fabes € B22 away bye 
2 4.9 Loma bs 4.6 2.0 3.9 Sok Zo 
3 3.9 fo 1.6 3H 20g Leo 
4 
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TABLE 3--Average depth of snow on April 16, 1970, at Upper © € 
oheep Creek exclosure., 


Treatment 

Position on Olope Brush Removec Srush Sprayed fo Treatment 
Near Creat 8.0 TO tie 
Midslope 8.4 10,4 10.8 
Near Base 6.4 i120 | Oso 
Average 7.6 10. 8 10,7 

Literature Cited: 

1, Dunford, =. G. 1967. Techniques in grassland watershed ¢ 


research, XIV, UFRO-KONGHRiZS5 Section ll, pp. 444-462, 


Oe unis os, 1970, Changes in water balance components 
with pasture management in South-Eastern Australia, Jour. 
of Hydrology 10:90-102, 

8. tisusby avis. Cot-019 70! <eHydrologie@ and Soi0 te eiiceuer on eo racine 
versus nongrazing near Grand Junction, Colorado, Geologic 
Survey Research, B232-B236, 


4, Meeuwig, R.O,. Infiltration and soil erosion as influenced ‘by 


vegetrtion and soil in northern Utrn, Journal of Aange Manage- 
ment 25 (3): 185-188. 


3, Meeuwig, R.O, 1969. Infiltration and soil erosion on Coolwater 
Ridge, Idaho, USDA Forest service Research Note INT-103, 
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CRIS Work Unit No. SWC-011-fBo-1 Code No. Ida-Bo-105.5 
Title: Field testing and evaluating mathematical models of 
a two-dimensional, infiltration flow system resulting 
from snowmelt. 
Location: Northwest Watershed Research Center, Boise, Idaho. 
Personnel; D. Ly, sochreiper, G. R.sStephenson, C. W. Johnson, 
L. M. Cox, G. A. Schumaker; R. W. Jeppson (Utah 
State University). 
Date of Initiation: July 1970 


Expected Termination 


Date: 


Originally planned - December 1972 
Recommendation - December 1972 


Objectives: 


1. 


3 


To field test a numerical steady-state two-dimensional flow model and to 
develop and verify a transient model by use of data from a highly 
instrumented watershed slope. 


To use the verified mathematical models and a sensitivity analysis to 
determine the influence of the various parameters (saturated conductivity, 
capillary pressure, bubbling pressure, and soil property coefficients) on 
the resulting flow patterns. 


To use the verified mathematical models to obtain solutions for other 
slopes on the Reynolds Creek Experimental Watershed. 


Need for Study: 


halt 


mati 
The 
resp 


Any 
not 


ecent years intensive studies of the various subsystems of the hydrologic 


_--cyele have been conducted. Research has been designed to establish mathe- 


cal descriptions of the physical flow processes within each subsystem. 
ultimate goal is a physically based, digitally simulated, hydrologic 
onse model to describe the entire hydrologic cycle over any watershed. 


comprehensive study of watershed precipitation-runoff relationships must 
neglect the flow system that occurs in the soil mantle. Input to such a 


flow system is generally provided by infiltration of snowmelt or rainfall 


thro 


ugh the soil surface. Recent studies have demonstrated that the 


infiltration characteristics of the soils and the partially saturated flow 


syst 
soil 


em in the soil mantle are the result of several interacting fundamental 
properties. The occurrence and magnitude of overland flow, quantities 


of interflow and baseflow, and recharge to the groundwater aquifer are all 


depe 
desc 


ndent upon the soil properties. If soil mantle flow systems can be 
ribed adequately, progress towards the goal of simulating the hydrologic 


cycle will be achieved. 
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This study should provide an important link for any comprehensive runoff ¢ 
prediction model that includes snowmelt on a watershed. Snowmelt is an 
important contributor to spring and early summer runoff from many watersheds 

in the Northwestern United States. Ultimately, this study will be used to 
develop more accurate methodology for predicting surface and subsurface flow 
from a watershed. Information concerning the flow system resulting from 
snowmelt infiltration is essential if Northwestern lands are to be managed 

for optimum multiple use. 


Design of Experiment and Procedure to be Followed: 


Steady-state mathematical solutions for a two-dimensional flow system 
resulting from infiltration on a watershed slope were analyzed initially in a 
1969 exploratory study and in project report PRWG-59c-1, “Numerical Solution 
of the Steady-State Two-Dimensional Flow System Resulting from Infiltration 
on a Watershed,’ by R. W. Jeppson, Utah Water Research Laboratory, Utah State 
University. It was found that the partial differential equation boundary 
value problems based on basic soil properties were within reach of solution 
capabilities. 


Since a steady-state model of a watershed cannot provide complete predictive 
information for the naturally unsteady watershed flow system, effort will also 
be devoted toward the development of a comparable transient mathematical model 
for the partially saturated flow system of a watershed slope. This 

transient model would include a source of moisture such as snowmelt or rain, 
evapotranspiration and deep percolation all incorporated as boundary condition 
in the partial differential equation initial-boundary value problem. 


In developing the model for the transient flow system and in testing its 
solution capabilities, the solution from the steady-state model will act as 
a "bench mark’? toward which the transient solution should approach in an 
asymptotic manner at longer time intervals. Consequently, the steady state 
solution results will also act as a check on the accuracy of transient 
solutions from more sophisticated mathematical models. 


In order to evaluate the effects of soil parameters, to establish the actual 
flow pattern, and to demonstrate the usefulness of the mathematical solution 
results, a watershed profile was selected and heavily instrumented. The 
northeast slope of the Upper Sheep Creek Watershed W-17, a subbasin of the 
Reynolds Creek Ex Sn Dy oe Watershed, was selected. The instrumentation sites 
are shown in Figure 1l~, a contour map of the Upper Sheep Creek Watersheds, and 
in Migsure 2 fsa neaaete. cross section of the slope. 


The instrumentation from which data were obtained during the water year 1969-70 
includes: (1) seven single probe and two dual probe sites for measuring soil 
moisture with a neutron probe, (2) two single piezometers at sites P-1 and P-2 
for measuring the potential surface of the saturated region, (3) three sites 
containing batteries of tensiometers at various depths for measuring soil 





1/ 


— All figures pertaining to this Research Outline are located in order G 
following page 8-10. 
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i ) ai) moisture tension beneath the snowpack, and (4) two identical drop-box weirs 

with cutoff walls for measuring stream-channel flow upstream and downstream 
from the instrumented slope. At the beginning of water year 1970-71 nine 
more piezometers were installed. Two of the additional piezometers were 
installed at site P-1 to provide a battery of three piezometers at various 
depths. The single piezometer that was located at site P-2 was discontinued 
in place of a battery of four new piezometers at various depths. A battery 
of three piczometers was also installed at a new site, P-3, in the middle of 
the slope. 


A seismic survey was conducted late in the water year 1969-70 on the 
instrumented slope, This survey was necessary to complement and complete the 
profile cross-section information obtained from driller's logs of the. soil 
moisture and piezometer sites. 


Experimental Data and Observations: 


Moisture content data for the soil profiles at the soil moisture and dual 
probe sites were obtained on a weekly basis throughout the year with a 
neutron probe. During the snowmelt season soil moisture data were obtained 
twice weekly. Gravimetric samples were taken occasionally to spot check the 
neutron-probe soil-moisture data and to determine bulk density, porosity, 
and saturated conductivity. Representative data for soil moisture, bulk 
a density, porosity, and saturated conductivity are listed in Table l. 
; ey) Representative soil moisture content profiles obtained from the neutron probe 
lp a data at sites SM-1, DP-2, and SM-5 are illustrated in Figures 3, 4, and 5, 
respectively. 


The two single piezometers from which data were collected during water year 
196S-70 were located at site P-1, close to the stream channel, and at site 
P-2, on the upper portion of the slope just below the usual snowpack area 
(see Figures 1 and 2). These piezgometcrs were used to determine roughly the 
potential surface of the saturated region. The hydrograph obtained during 
the snowmelt season from the lower piezometer P-1 is illustrated in Figure 6. 


Three multipie-cell tensiometers were located in the soil under the snowpack 
to measure capillary pressures at various depths in the soil. This data is 
listed in Table 2. 


Periodic surveys were made during the snowmelt season to obtain the location, 
volume, and water equivalent of the snowpack. The average density of the 
snowpack remained essentially constant (0.49-0.55) during the snowmelt season. 
Figure 7 illustrates the changes in a typical cross-sectional area of the 
snowpack during the winter and spring of 1970. 





The snowmelt data were analyzed by two methods. One method was used to 
indicate only the total volume of water that melted from the snowpack. each 
day. The second method was used to derive a soil surface intake rate by 
accounting for the decreasing ground surface area beneath the snowpack. 

hace A, Daily melt rates from the cooperative ARS-SCS snow pillow site at Reynolds 
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TABLE 1.--Soils data from the north slope of Upper Sheep Creek Watershed W-17. 
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Water Content, % volume 














Neutron Probe Gravimetic donip eset Average Average Average 
Site SM4 Site SM5 Bulk Poro- Saturated 
Danie goer /049292705 = 9-72-70 9959-70 9-2-9 3-10-70 Density; sity,- Conduc- 
feet Actual Average Actual Average em/cm a eres 
em/min 
0.0-0.5 2a 6 17.8 94 64.6 0229 
97 a5 0 1 aw L6n2 .85 66.9 . 04.4.1 
Lies 13.06 12,4 14.9 Pee Load 
1.0-1.5 15.5 Lie a0) 66.2 . C396 
1. 34220 16.0 21.6 eed. 66.0 . 0540 
13756 pete 14.4 15.5 i720 22° 6 
2.052 .5 ay al 23.5 . 84 638.3 . 0463 
25-370 1D a2 25.7 . 89 66.4 . 0405 
1502 LGm'l 16.5 19.0 19,2 24.6 
3.0-3.5 pee oat G9 66.4 .O0G19 
3.5-4.0 21, 92/ 27. 92! at ea A 
17.9 18.4 24.1 fame pe 3) 
1/ ' ; - 
~ Samples were talien from a pit located approximately 25 feet west of piezometer battery P2. Dry- 
weight moisture contents were determined and converted to volume basis. 
2/ 
Samples of constant volume (1306.6 cm?) were taken from the pit mentioned in footnote 1/ with a 
Soiltest core sampler. 
3/ 


Values were computed using a value of 2.65 for real specific gravity. 
4/ ‘Values were determined ‘using -25-gram samples -in constant-head. pexmeameters. 


a/ Value was left on dry-weight basis, since bulk density value was missing. 
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TABLE 2.--Soil moisture tension data from the north slope of Upper Sheep 


Creek Watershed W-17. 
nn REE EEREEEEEEEEEEEEEmEEEeeeeneeeee ee 
V 
Tension in Inches oz Water below Tensiometecr Depth 





Mili- Site T-1 2/ Site T-3 
1970 tary Tensiometer Depth Tensiometer Depth 
Date Time 6.00 18.25 2.75 7.00 15.25 19,25 
ae POO Se Rie oa oe ee 6.77 0.145 Ht HY Te” Gk 
L05Sei, ee f Sh ere ees. 9.96 +12.35" + 9.48 16.97 
1250 eaak bis jaar Pres. 25 + 9,23 16.02 
B-25io0O3014, ea lank . yovigadg .2 PUL 79N2 NE OF 20.52 22.67 
11G.sige: eat work andy 201G POs 250 72 21.08 O17 
1620 1926 1496499 21.90 31.08 23.06 24.83 
5-26 1015 12.39 19914.94 27.15 - 42.90 rd Bee 20.07 
1100 + 0.24 + 9.8 272 297k2 43752 MT ony Mae anid chs Be 
1205 + 6.24 + 8.49 27.65 44.52 cle he Wi ie Be, 
1300 +24.24 + 7.61 27083. 2045197 es Bebe es Be 
14.00 +4224 416.43 20-02 “46752 28.02 25.77 
1500 442.24 + 5.49 28.40 47.40 28.27 a ika? 
1600 442.24 + 4.70 261402 Sagto7 28.52 27.65 
1700 +4252G Cer Ait 26252°° 49. 15 53077 +10.48 
5-27 1445 439.99 +20.35 FINO es 2002 31.90 +19.73 
5-28. 1310 +40. 49 luenbye, 3 ei 20d TB, RT 34.77 £2323 
1445 440.37 454.24 34.33 63.08 se, FF +21.29 
5-30 1650 +3970 53124 40.579 +°9.48 +10.48 +20.98 
6-4 1230 +38.76 +50.17 47.93 12.24 42.99 +16. 69 
6-8 1345 3-35.28 449.45 SR. 25 easly 0.51 +10.45 
6-9 1200 +24.74 445.49 SEP SLION 16.69 + 3.86 + 7.33 
6-10 1200 +33.84 441.99 SAG Poae 1:94 2 + 3,86 + 7.33 





1/ ; : : 
— Tensiometers at Site T-2 did not function properly. 


& - 2 . 7 & 
— Depth is in inches below ground surface. 


3] 


= Plus (+) sign indicates positive pressure head in inches above 
tensiometer position. 
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Mountain, 4.0 miles southwest of che snowpack on the north slope of Upper Fé 
Sheep Creek Watershed W-17, were used as a guide in determining daily melt 
rates at the latter site. The Reynolds Mountain snow pillow site is 
comparable to the snowdrift area on the north slope of Upper Sheep Creek. 
Exposure is approximately the same, and the elevation, 6780 feet, is only a 
difference of 500 feet. Snow melt and rainfall data from the Reynolds 
Mountain snow pillow site and rainfall data from the Lower Sheep Creek 
Weather Station (1.5 miles northwest of the north slope of Upper Sheep Creek 
and at elevation 5410) and from the Black Mountain rain gage (0.7 miles 
northeast of the north slope of Upper Sheep Creek and at elevation 6610) are 
listed in Table 3. Results of the two methods of analyzing the snowmelt data 
are illustrated in Figure 8. 


Stream-channel flow is measured upstream and downstream from the instrumented 
slope by two identical drop-box weirs with cutoff walls. No overland flow 

or surface runoff from the slope was observed during the snowmelt period. 

hannel flow is attributed to subsurface flow from snowmelt infiltration with 
subsequent seepage into the stream channel. The difference between the 
hydrographs from the two weirs is illustrated in Figure 6 and gives an 

indication of the snowmelt hydrograph after it is attenuated by flowing downslope 
through the soil mantle system. 


Comments, Interpretations, and Future Plans: 


A comparison of Figures 6 and & indicates that equilibrium or steady-state 
conditions have not been fully achieved. The snowmelt hydrograph peak precedes\ 
the lower piezometer hydrograph peak by two weeks,.and the latter precedes the 
channel-runorfi peak by five days. 


The adequacy of the steady-state mathematical model developed by Jeppson for 
this study in defining variations in hydraulic conductivity or other soil 
parameters from an unsteady field situation needs further study. In brief, 
more field measurements are needed. However, the model does indicate what 
types of additional cbservations should be obtained. 


Ac the time the computer problems were formulated, not att of the data needed 
to completely define the field flow situation were available. For example, 
the two initial batteries of tensiometers at sites T-1 and T-2 were located 
during the autumn of 1969 just below the expected extent of the 1969- 70 
snowpack accumulation. However, the snowpack did not extend downs lope. as far 
as had been anticipated. An additionai battery of tensiometers were installed 
at site T-3 beneath the snowpack during the spring of 1970, but prior to the 
heavy snowmelt period. Furthermore, measurements to determine the lateral 

and vertical variations of saturated hydraulic conductivity were not complete. 
In fact, the work plan called for using the steady-state model to indicate 
what field measurements would be most useful. Consequently, judgment was used 
in estimating some of the parameter values needed to define the computer 
problems. The other parameter values were varied in a systematic fashion to 
bring computer solutions to closer agreement with certain observed features of 
the actual flow system. é 














TABLE 3.--Comparable snowmelt and rainfall data. 





Snowme lt Rainfall Rainfall Rainfall 
at Reynolds at Reynolds at Lower at Black 
Mountain Mountain Sheep Creek Mountain 
May Snow Pillow Snow ae ete MAD ERED Rain 
1970 Site‘ Site; Station> Gage=4 
Date Inches Inches Inches Inches 
1 0.10 is 
2 0435 : 
3 0. 84 : : 
4, 0.90 
5 0.174 Cz 
G 0.69 0.07 0, 04. 0.03 
7 0.66 neers iors te ae 
3 0.92 0.43 0. 09 Os 13 
9 0439 die 0.09 0.06 
10 Oal3 a29 Oul2 0,24 
ih. 0.02 0.05 a shchs evens 
12 0.01 Qoiatds 0.02 0.03 
i3 0. 00 0.09 0.03 0.02 
14 0..03 sit ns 2 0. 04 
15 0.389 
16 Lae 
Ly 1.27 
18 ry mee oe aes is & pears 
rate, Oe) gets a ~~ 0.04 
20 0.52 Qa 1D ialel 0.08 
Zs Late Pec eA 0. 02 
22 0.83 0.03 Oude oe 
ea 1.49 0.35 0. 34 Pa 
24. Zee l ire Piccks 0.01 
ya) 2.60 sap ee T Piaget 
26 appe he) 0, Oy. D032 0.06 
27 Lege 9 5 de ass cae 0.02 
C O¢s5 5 dha fo 
29 0.01 oar 


L/ 


Gage 176407 located at elevation 6780 and 4.0 miles southwest of the 
north slope of Upper Sheep Creek Watershed W-17. 


Gage 127407 located at elevation 5410 and 1.5 miles northwest of the 
north slope of Upper Sheep Creek Watershed W-1/7. 


Gage 128487 located at elevation 6610 and 0.7 miles northeast of the 
north slope of Upper Sheen Creek Watershed W-1/7. 
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On the basis of the soil tension data (Table 2) and the soil moisture 
measurements (Figures 3, 4, and 5), the capillary tensions were specified 
in the computer problems as either 0.5 or 1.0 feet of water. Drillers' 
logs indicated that soil material obtained during the installation of the 
piezometers and soil moisture cubes contained an increasing amount of fines 
with depth. Consequently, a function allowing decreasing values of hydraulic 
conductivity with depth below the ground surface was used in practically all 
o£ the computer solutions. 


The values for the boundary coordinates in the computer solutions were 
selected upon examination of the field data. Piezometer data were not 
adequate to define the position of the water table throughout the watershed 
slope for a steady-state flow system; however, the saturation data obtained 
from soil moisture measurements suggested that the water table created by 
snowmelt would be a relatively small distance below the ground surface. 
Therefore, boundary coordinates from the base of the snowpack to the stream 
channel were selected to give a water table at a depth of one foot or less. 
Boundary coordinates for the impervious layer were selected to give depths 
to the restrictive layer as suggested by the drillers' logs and as shown 

in Figure 2. 
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t is not practical to include all of the computer solution results in this 
report. Flow nets tor four solutions are presented in Figures 9, 10, 11, 
and 12, to facilitate the explanation of factors that appear to influence 
the downslope flow within a watershed profile. Every other streamline and 
equipotential line that was computed is plotted om these flow nets. The 
dashed line separates the saturated and the partially saturated flow regions 
(flow above the dashed line is partially saturated, whereas, the flow below 
the dashed iine is saturated). 


The first two flow nets, Figure 9 and 10, present the solution to problems 
with all parameter specifications identical, excent for the variation in 

the saturated hydraulic conductivity, K_. Furthermore, in both problem 
specifications the change in saturated Rydr caulic Pekduceseuey with respect 
to the stream function is identical such that K_ at the ground surface 
equals 1.75 plus K_ just above the impervious layer, and such that Ky varies 
linearly between the two surfaces. 


In the problem illustrated in Figure 9, no variation in the hydraulic 
conductivity was specified with respect to the potential function. In essence, 
the soil at the bottom of the slope was assumed to have similar hydraulic 
properties to the soil at the top of the slope. Since the shape of the 
resulting ground surface profile illustrated in Figure 9 is steeper near the 
channel and flatter near the top than that of the actual watershed as shown 

in Figure 2, some parameter specification or combination of specifications 
needs to be modified, 


While a number of the parameters used to specify the problem may be changed 

to allow closer agreement between the computer solution and the actual profil 
Shape, the most logical change seems to be to vary the saturated hydraulic G 
conductivity, Kj, such that its value is greater at the bottom than near the 
top of the slope. 
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Solution to several problems (not illustrated in this report) were obtained 
in which Ko was specified to vary linearly with the potential function and 
to be approximately twice as great at the bottom as it was at the top of 
the slope. With only linear variations in K, and without variations in the 
anisotropy parameter or changes in the boundary coordinates, it soon became 
apparent that complete agreement between computed and actual profiles could 


not be achieved. 


Regions of nonconformity, except near the stream channel, were eliminated 
systematically by superimposing positive and negative sine curves on a 

decreasing linear function of K, with potential function. A flow net for 

such a formulation is illustrated in Figure 10. These results indicate 

that the actual saturated hydraulic conductivity varies with distance up 

the slope from the channel bottom. More precisely, K decreases from the 
streambed up the slope for about one-third the distan Ee, then just below ¢ 
the dual probe site DP-2, it increases sharply. Ne aes decreases again 

until near the top of the slope, where it increases again. ; 


The sharp increase in saturated hydraulic conductivity required by the results - 
of the solution in the vicinity of the dual probe sice DP-2 may result from 
the actual decrease in surface soil depth at this location. Furthermore, 
saturated hydraulic conductivity data obtained in the laboratory support 
this hypothesis. Values of K, were obtained from Sora tny bas soil samples 
taken from three pits spaced Eon the slope. Samples taken from a pit 
near site DP-2 yielded K_ valucs that were _ peer eie twice as large as 
K, values obtained for che soils in pits located near the piezometer sites 
Pe land P= 2, 





It is noteworthy that other parameter specifications and variations could 
have been used to bring solution results in agreement with the actual case. 
Agreement couid probably have been PRIPneY 29, by varying the anisotropy 
parameter in the upslope direction, or by specifying different boundary 
coordinates. 


Different boundary coordinates were specified for other solutions, in which 
the water table was maintained within a depth of one foot from the ground 
surface, and in which the capillary tension head in the soil beneath the 
snowpack was maintained at one foot oi water. The flow net that is 

illustrated in Figure 11 resulted from such a solution. In this and other 
such solutions not included here, only linear decreases of K, with potential 
function were defined. Of these solutions the profile sepa in Figure Ll 
corresponds best to the actuai profile. However, this srofile does rise 
above the actual profile just upslope from the dual probe site DP-2, at the 
position where the restrictive layer reduces the depth of the ran, soil. 
Again with these coordinates, the actual profile cannot be duplicated 
unless K, is sharply increased in the vicinity of DP-2 while the anisotropy 
parameter is maintained constant, or vice versa, or a combination of the 
two parameters is taken. 


The solution, illustrated in Figure 12 used still other boundary coordinates 
in its specifications. Here it was assumed that a slightly larger increase 
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in the saturated hydraulic conductivity with stream function exists than 
in previous solutions, and that K_ does not vary with potential function. 
Saturated hydraulic conductivity in the horizontal direction was assumed 
to be twice as great as in the vertical direction. Therefore, this 
sokution illustrates the shape that the watershed vrofile would have under 
steady-state flow conditions with the depths of soil shown, with the entire 
water table close to the ground surface, and with saturated hydraulic 
conductivity characteristics the same at the top as at the bottom of the 
slope. The computed profile is much steeper in the areas of shallower 
soil than is the actual watershed. In these shallow soil regions either 
the water must emerge from the surface or K, must be greater than in the 
remaining profile. Should the water emerge from the surface, it would 
tend to make this area more susceptible to erosion, and thus provide a 
natural mechanism for further reduction in the depth of soil capable of 
transmitting moisture. : 


Future plans include the development and verification (with field data) 

of a numerical, transient, two-dimensional flow model. R. W. Jeppson, Utah 
Water Research Laboratory, will develop the model, and the Northwest 
Watershed Research Center staff will obtain the necessary data for 
verification from the 1970-71 water year. To meet these needs, the new 
piezometer batteries mentioned in an earlier section of this report were 
installed. Tensiometer batterics will be installed at the base of the 
snowpack prior to the snowmelt season. Since it was demonstrated last year 
that good tension data are difficult to attain, an effort will be made to 
modify the tensiometers. An ideal situation would include continuous 
recorders and no hydraulic lines on the tensiometers,. The recording rain 
gage 130031 that is located very close to the study area (Figure 1) will be 
modified so that it will provide reliable data (no data were available for 
water year 1969-79). Further geologic and hydrologic investigations will 
be conducted to determine if all the water from the large snowpack on Upper 
Sheep Creek Watershed W-16 flows through the drop-box weir W-16 before 
entering the stream channel of Upper Sheep Creek Watershed W-17. 


A manuscript in preparation since last year's report is the following: 


JeppsonzOR. 2W295 Schreibery eb.’ ii. Stephenson, G. Rh. Jonnsou, ue ta. 

Cox;°L. "Mi, and Schumaker, G?-A.°°1971.-" Solution of a watershed 

flow system resulting from snowmelt with verification by field data. 

Abstract accepted for presentation at the 1971 Annual National 

Summer Meeting of ASAE, June 27-30, Pullman, Washington. Manuscript 
1 


i 
to be submitted to SWC for anproval to publish in the Transactions of 
the ASAE. 
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Figure 1, Instrumentation and topography map of Upper Sheep 
Creek Watersheds W-16 and W-l7, 
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Soil moisture content profiles obtained with neutron probe at Site SM-l, North 
Slope Upper Sheep Creek Watershed W-17, 
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Figure 6, Hydrographs for piezometer P-l and stream channel, 
North Slope Upper Sheep Creek Watershed W-17, 
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Changes in snowpack cross-sectional area at station 5 + 00, North atone Upper Sheep 
Creek Watershed W-17, 
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Figure 10, Computed flow net resulting from a linear variation of saturated hydraulic 
conductivity with stream function and a nonlinear variation (positive and 
negative sine curves superimposed on a decreasing linear function) with 

- potential function, waive | 
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Figure 12. Computed flow net resulting from different linear variations of saturated 
hydraulic conductivity with stream function and with potential function. 
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CRIS Work Unit No. SWC-Ol1-fBo-l Code No. Ida-38o0-105. 6 
Title: Developing, testing, and evaluating an analytical 


infiltration model, 


Location: Northwest //atershed Research Center, Boise, 
Idaho, 
Cooperation: The Utah State University Water Research 


Laboratory will cooperate in the analytical. phases 
of this study. The Oregon State University Agri- 
cultural fngineering Department will cooperate in 
assessing soil properties and in laboratory and in 
field testing. The Bureau of Land Management, 
USDI, will contribute support and participate in 
pursuit of the study. 


Personnel: Ljutimpcareluer, ai. Bh. Hammon, and ts. A. 


Schumaker (ARS-C:/C); R. W. Jeppson (Utah State 
University); and R. H. Brooks and =, N. Biggs 
(Oregon State University). 


Date of Initiation: February 1971 


Expected Termination 


Date: Field Work: November 1974, 


oo 


Laboratory: | september 1972, 
Interpretation and Summary: December 1974, 


Objectives: 


To develop: or adapt mathemetical models in the form of partial 
differential equations to descrise steady-state and transient 
one-dimensional and three-dimensional axisymmetric flow through 
partially saturated soils. 


To determine quantitative means by use of the mathematical 
models to adjust for lateral ''spreading effect'' of moisture 
movement from a circular rainfall simulator under various soil 
types and conditions. 


To test and refine the mathematical models by comparing results 
with laboratory and field determinations of infiltration and to 
establish the relative influence of the several interacting physical 
processes on infiltration, 
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4, To determine the parameters in the saturation-pressure re- « 
lationship (moisture characteristic) and saturated conductivity 
by the use of parameter optimization in conjunction with the 
mathematical models and infiltration data. 


9. To determine the potential quantity of water retainable by 
various soil-vegetation complexes that is independent of in- 
filtration for different initial soil moisture levels. 


Need for Study: 


Infiltration of water into soil profiles is the critical hydrologic com- 
ponent in watershed management, overland flow prediction, sediment 
generation, natural and artificial ground-water recharge, and irrigation, 
With increased public interest in the management of land and water 
resources, better analyses and predictive methods for describing infil- 
tration and flow in porous media are needed. 


This study should provide an important link for the development of a 
comprehensive hydrologic response model. More accurate information 
concerning the flow system resulting from infiltration is essential if 
agricultural lands are to be managed for optimum multiple use. ~ Ps 


Much data have been collected around the country in past years using 
small infiltrometers on many different soil types. This study should 
provide a quantitative means to adjust the data for the lateral spreading 
effect of the water beyond the boundary of the ring. Such a correction 
for lateral spreading would not only increase the value of the reams of 
infiltrometer data that are available but make it possible to adjust field 
data to meaningful quantitative values. 


Several studies that have been’ conducted within the past 20 years have 
fully demonstrated that infiltration characteristics are dependent upon 
fundamental hydraulic properties of soil in relation to water movement. 
These studies have demonstrated that, by defining such scil properties 
adequately and describing their effects upon the flow system by means of 
partial differential equations, solutions can be obtained which show close 
agreement with laboratory and field oiservations. Solutions to well- 
formulated boundary and initial value problems permit various components 
and features of the flow system to be isolated and studied, thus provid- 
ing knowledge and insight into extremely complex and otherwise unde- 
cipherable cause-effect relationships, 


The partial differential equations that describe partially saturated flow in «& 
a porous media are nonlinear. Consequently, the majority of available 
solutions have been obtained for the assumption of one-dimensional flow. 



































The more general boundary and initial conditions have required inves- 
tigators to obtain solutions by numerical methods, Solution of the 
partial differential equations requires a relationship of saturation and 
permeability to capillary pressure. Considerable research has been 
devoted to defining relationships between these variables in a partially 
saturated soil-water-air flow system. 


The problem which will be studied under this research outline is the 
one-dimensional and the three-dimensional axisymmetric Partially satu- 
rated flow regimes that result from moisture applied at the surface of 
a soil core and on plots beneath a rainfall simulator. This study will 
provide needec analytical methods for describing the infiltration com- 
ponents as part of the total effort by the Northwest Watershed Research 
Center in developing a physically-based, computer-simulated response 
model of the entire hydrologic cycle of a watershed. 


If the hydraulic properties of soils are defined adequately, resulting 
effects upon flow systems in soils can be described by solutions: to the 
defining partial differential equations. Various components of a flow 
system can be isolated and studied by formulating proper boundary- 
and initial-value problems. Knowledge and insight can be provided into 
extremely complex, undecipherable relationships. 


The theory of porous media flow used in this study is also applied to 
describe the two-dimensional flow system resulting from melting snow 
over the upper portion of a watershed slope (Research Outline Ida-Bo- 
105.5). The equipment and some of the procedures developed for this 
study will be utilized in determining the hydraulic properties of soils, 
particularly the moisture characteristic, as a cooperative study with 
Oregon State University. 


Design of #xperiment and Procedure to be Followed: 


l. Variables-- 


Hydraulic conductivity of the various soil layers. 
Moisture content of the various soil layers. 

soil moisture tension. 

Soil profile and physical characteristics. 

. Rate of water application. 

. Free water detention of soil-cover complexes. 


oad ww wv 


rh 


2. Data to be Obtained-- Data will be obtained for the variables listed 
above. Hydraulic conductivities, roisture contents, and moisture 
tensions will be obtained in 3 different ways: (1) under carefully 
controlled laboratory conditions at Cregon State University, (2) under 
field-simulated conditions at the Reynolds Creek Rainfall Simulator 
Laboratory, and (3) in situ; The soil profile and physical charac- 
teristics such as porosity and bulk density will be determined by 
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digging open pits and taking gravimetric Samples. The rate é 


of moisture application can be controlled and will be varied 
systematically to determine its effect upon the resulting flow 
system, 


Procedure--Iliathematical models of steady and transient one- 


dimensional flow and three-dimensional axisymmetric flow through 


partially saturated soils are currently being perfected by Dr, 

R. «/, Jeppson at the Utah State University Water Research 
Laboratory under a cooperative study. These models will satis- 
fy objective No. l, 


Objective No, 2 is being pursued concurrently with objective No, 
1, <«©reviously published soils data, not necessarily from the 
Reynolds Creek Experimental :’atershed, will be used initially to 
investigate the importance of lateral spreading. As more data 
are obtained from the Reynolds Creek Watershed soils, they will 
be used to verify the model and also to satisfy objective No. 3. 


The laboratory phase of objective No, 3 will be acooperative 

effort with Oregon State University under a cooperstive study. 

The laboratory investigations will be conducted at Oregon State 
University and at the Reynolds Creek Rainfall Simulator Facility. 
Disturbed and undisturbed soil cores from Reynolds Creek Experi- 
mental ./atershed will be tested in the laboratories. Results will 
be compared to those obtainec from the mmthematical models, 


Following completion of the lnboratory phase of objective No, 3, 
field testing phase will be initiated. The rainfall simulator with 
its gamma probe infiltrometer is portable, so it can be trans- 
ported to various sites on the Reynolds Creek ixperimental 
Watershed. Data similar to that obtained in the rainfall simu- 
lator laporatory will be obtained in situ for several different 
soil-vcgetation complexes, Clerieh 


Since soil moisture tension and conductivity data are difficult to 
obtain, especially in situ, parameter optimization will be utilized 
to achieve objective No. 4, <A suitable analytic or parametric 
mocel for the saturation-capillary pressure relationship is to be 
obtained in a cooperative study with Oregon State University. 
Results from the laboratory phase of objective No, 3 will be 
utilized in a parameter optimization technigue, using the mathe- 
maticol model and a modified Surdine Theory as on alternative 
to using soil moisture tension and saturation concuctivity data. 
If parameter optimization can successfully be usec to replace 
soil moisture tension and saturated conductivity date, then ensu- 
ing studies would not have to include such data, 





« 
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In order to separate out the potential quantity of water retain- 
able by soil-vegetation complexes and the potential guantity that 
is independent of infiltration, the rainfall simulator will be used 
to apply a given quantity of water at the same initial soil mois- 
ture conditions for a range of time intervals, Data obtained in 
this manner will be used to extrapolate to an instantancous eten~ 
tion of water for the various soil-vegetation complexes, 


Experimental Data and Observations: 


1. Instrumentation--The conceptual anolytical infiltration model was 
conceived and forraulated on the premise of available or developable 
mathematical models, access to large clectronic computers, and of 
acquirable new instrumentation for field mensurement of the necessary 
hydraulic and physical soil properties. 

A rainfall simulator-gamma probe-infiltromecter has been built as a 
cooperative effort with the University of Idcho and put into operation 
during 1970. The 6 by 6-foot unit is composed of 2 by 3-foot double , 
compartment modules supported at a height of 8 feet by 2 collapsible 
frame, Capillary needles on 3 by 3-inch centers extend frorn the water 
compartment through the air compartment and an orifice in the lower 
plave from which air flows for controlling drop size. ., Controls sare 
available for regulating the water and air pressure. | 


Two sets of modules have been built with different sizes of capillary 
needles, The eae needles, 0.027 inch (I.D,), will deliver water 
at a rate of 0.5 to 4,0 inches/hour, and the smaller needles, 0,016 
inch (T7bs): ae deliver water at a rate of 0.15 to 2.0 inches/hour, 

with drop sizes controlled by varying the sir pressure, 


The two-probe (gamma) density gage is operated by a pulley and drive 
mechanism from the top of the rainfall simulator. Access tubes extend 
through the modulcs and connect to access tubes in the soil, The con- 
trol system makes it possible to obtain density changes (water content 
changes) simultaneously with the movement of water into the soil pro- 
file, Soil water content, therefore, can be monitored at tensiometer 
locations for devclopment of capillary presaure-saturation relationships. 


The $roxler Modcl SC-10, two-probe (gamma) density gauge, first used, 
was found to be too temperature sensitive to meet field measurement 
requirement. That gauge has been replaccd by a Troxler Model 2376, 
two-probe density gauge, which utilizes a trickirig differenticl, > Pulse 
Height Discriminator eee system, The PHD operates with a scintil- 
lation probe consisting of a photomultiplicr tube which is Baticrtisy and 
mechanically coupled to a Thallium activated Sodium Iodide crystal 
doped with Aimcricium 241, The function of the PHD is twofold: 
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(a,) It selects pulses of a predetermined height and € 
Within a pecborlehines window to allow density 
measurements with small composition error. : 


(b) It compensates for changes in the gain of the 

photomultiplier tube and the characteristics of 

the Sodium Iodide crystal due to temperature, 
Two tests were made with the new -H=O system to determine its 
temperature stability. First, a series of density readings were 
made in the laboratory at 3 temperature levels using the magnesium 
standard, Second, an actual infiltration test was run in the labora- 
tory with the probe remaining stationary at 4 inches below the soil 
Snes 


Data obtained in the first test are tabulated in Table 1. The means 
and standard deviations for the 3 samples and for the combined data 
are tabulated at the bottom of the table, 


The second test of the PHD system was made in connection with a 
laboratory infiltration test on a large disturbed soil core to observe 
capillary pressure and saturation changes with passage of the wetting 
front, The counts as a function of time for the range of saturation, « 
99 to 85 percent, are shown in Figure 3,—' The corresponding change 
in saturation anc capillary pressure is plotted in Figure 7. The gamma 
probe was stationary during this test at 4 inches below the soil sur- 
face, 


Capillary press Bue or soil water tension measurements, such as those 
in Figures 4 anc 7, were obtained with ceramic cups and pressure 
transducers, This equipment will be used for tension measurements 
up to 0,8 bar, and psychrometers will be used for tension in excess 
of 1-2 bars. ‘Scanivalves are available to alow multiple readings of 
tensiometers, but experience has demonstrated that a sufficient time 
must be allowed for equilibrium in the system before reliable readings 
are obtained. 


A small laboratory facility has been built at the Reynolds Creek Experi- 
mental Watershed Headquarters to accommodate large soil cores with 

a 5-foot deep pit the size of the reinfall simulator or 6 by 6 feet. 

The rainfall simulator-gamma probe-infiltrometer is portaple: a and will 
be used in both the laboratory and field studies. 
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10179 10227 10189 
10064 10203 10144. 
10139 10224 10167 
10154 10239 10264 
10194 10367 10069 
10169 1024.7 10276 
10096 10164 10266 
10254 
10279 10202 10402 
10213 10348 10122 
10084 10071 10092 
10105 1034.6 10267 
10181 10274, 10164 
10128 10165 10139 
10280 10123 10288 
10283 10144 10227 
10262 10371 10122 
10345 10193 10202 
10173 10230 10115 
10264 10164 1041 

10248 10301 

Z, = 10193 Hy = 10228 Zz = 10212 
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2. Laboratory Uata--To date, moisture content data, capillary pressure q 
(moisture tension) data, and hydraulic conductivity data have been obtained 
by two different methods using disturbed soil samples extracted from a 
pit on the Summit -/atershed W-12, a subsasin of the Reynolds Creek Ex- 
perimental v/atershed, These data have been obtained (a) by using a 2- 
inch diameter core under carefully controlled laboratory conditions at 
Oregon State University and (b) by using 14-inch and 60-inch diameter 
cores at the Reynolds Creek Rainfall Gimulator Laboratory. 


The capillary pressure-saturation data obtained for the 2-inch cores are 
illustrated in Figure 4 and are listed in Table 2. The data in Table 2 
have been smoothed to provide for a continuous polynomial curve-fitting 
between each adjacent three data values, The intrinsic permeability- 
capillary pressure data obtained for the Z-inch cores are illustrated in 
Figure 5, 


The 14-inch core was packed with the disturbed Summit soil sample to 

see if the results that were obtained fromm the 2-inch core under carefully 
controlled laboratory conditions could be duplicated under a field-simulated 
situation with the rainfall simulator and its associated gamma probe, The 
14-inch diameter was necessary to allow for the span of one foot between 
“the gamma probe and its detector. \/ater was applied by the rainfall & 
simulator over the entire 14 inches of surface diameter to provide a ver- 
tical, one-dimensionel moisture movement, 


Following the test with the 14-inch core, a 60-inch diameter core was 
packed with the disturbed Summit soil sample to obtain data to verify 
Jeppson's mathematical model of transient, three-dimensionel exisymmetric 
flow through particlly| saturated soils, since water was only applied to an 
area having a 14-inch diameter, data were obtained to quantify the lateral 
"spreading effect'' of moisture moverment from a circular application area. 


Three tests using the 60-inch sample were run, During the first test 
tensiometers were located near the center of the sample and the water ap- 
plication area at cepths of 3 and 12 inches below the soil surface. Addi- 
tional tensiometers were located at a depth of 3 inches and at a radius of 
9 inches (2 inches beyond the area receiving water), 12 inches, and 15 ... 
inches from the center of the sample. The effect of lateral. spreading 

at the 3-inch depth is illustrated in the saturation-time curves in Figure 
6, 


A second test was run with the 60-inch sample four weeks after the first 

_test. The second test was deemed necessary to perfect the techniques and 
/instrumentation used to collect the data. In the previous test quite a lot « 
of scatter was noticed in the capillary pressure data. This scatter was 
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attributed to excessive switching of the Sanivalve used to connect the 
hydraulic lines leading from the tensiometers to a single pressure 
transducer. To eliminate the excessive switching of the Scanivalve 
required in the first test, only two tensiometers were used during the 
second test. These tensiometers were located at a 3-inch radius from 
the center and at depths from the surface of 3 inches and 9 inches, 


respectively. cata scatter was essentially eliminated in this test, 


A third test was ceemed necessary, since saturation-time curves ob- 
tained for the second test with the gamma probe did not correspond closely 
to the saturation-time curves derived from the capillary pressure data 

(for the second test) and the capillary pressure-saturation curve (Figure 4) 
developed from the 2-inch core. The discrepancies noted in the results 
of the second test were attributed to the possiffility that the gamma probe 
and tensiometers were not at corresponding elevations in the sample when 
readings were taken, The third test was run 6 weeks after the second 
test, Two tensiometers were located near the center of the sample and 
at a depth of 4 inches below the soil surface, Extreme care was used in 
placing the gamma probe exactly at the same elevation as the two ten- 

_ siometers., —uring the test the gamma probe was not moved from that 

/ position, A better agreement was achieved between satursation-time curves 
obtained with the gamma probe and derived from the capillary pressure data, 
The capillary pressure-time curve obtained with the tensiometers and the 
saturation time curve obtained with the gamma probe are illustrated in 
Figure 7, 


3, Mathematicel Ifodel Data--The solution to the partial differential equa- 
tions that describe transient, one-dimensiona] vertical flow. and transient, 
three-dimensional axisymmetric flow through partially saturated soils depends 
upon obtaining values for (a) the change in saturation with capillary pres- 
sure, (b) the relative hydraulic conductivity (ratio of partially saturated 
hydraulic corductivity to gaturated.hydraulic conductivity), and (c) the change 
of relative hydraulic conductivity with capillary pressure. Of these needed 
quantities, it is more practical to obtain capillary pressure-saturation data 

in the laboratory or field with tensiometers and either Neutron meters or 
gamma probes. #ycraulic conductivity-capillary pressure data are much 
more difficult to obtain. Consequently, it is desirable to obtain values 

for relative hycraulic conductivity and the change of relative conductivity 

with capillary pressure from saturation-pressure datas or from some func- 
tional relationship involving pressure. ‘Such functional relationships have 
been proposed in the literature; i.e., the Erooks-Corey Equation. This 
relationship dic not give close agreement to the data obtained from the . 
soil taken from the Summit on the Reynolds Creek Experimental Watershed, é 
~ However, reasonable results are-obtained by using a modification of the 
Surdine Theory with capillary pressure-saturation data. . 
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pressure data (for the Summit soil) with that given by the Brooks-Corey 
Equation, A comparison of observed hydraulic conductivity-pressure data 
with that obtained py the Brooks-Corey “quation and with that obtained 

1 | from saturation-pressure data through the modified Burdine Theory is 
illustrated-in igure 9, Since better agreement to laboratory data is ob- 
ri tained by using the modified Burdine Theory, it has been incorporated in 
if the mathematical model, rather than in the Brooks-Corey “quation, 


, Illustrated in Figure 8 is a comparison of the observed saturation = 
i 
| 


\ 


In order to evaluate quantitative effects on the flow pattern from water 
applied at the curface of a circular infiltrometer, comparisons are needed 
between results obtained from one-dimensional vertical flow and three- 
dimensional axisymmetric flow models. The formulation cand solution 
method for both models are consistent and compatible. For the data ob- 
tained from the Summit soil, results obtcined from the two models 
indicate that boundary effects on circulor infiltrometers significantly alter 
) the flow pattern, even reducing the snturction at the surface centerline 

} oppreciably over thet which would exist for the same application rate over 
an infinite area, 

















O Figures 10 and -11 illustrate the capability of the transient, three-dimen- 

» —’' sional axisymmetric flow model, Isosaturation lines are plotted at two 

different time steps, r = K t (K_ is saturated hydraulic conductivity and 6 
t is time). For the problém illistrated in Figures 10 and ll, the- initial 

1 capillary pressure head’ inthe soil, h’, -is*'-9;0°feet of water, and the 

water is applied over oa circle of one-toot radius at a rate less than the 

intake capacity of the soil, Consequently, nowhere does the soil become 
completely saturated, 





In addition to Figures 10 and ll, Figure 12 illustrates in another manner 
the effect of latercl spreading from «2 circular infiltrometer, The satura- 
tion-time curves cre for the 3-inch depth end for radii which are con- 

4 sistent with the curves illustrated in Fisure 6 for the laboratory data. 


re) 


Comments, Interpretations, and Future /lnns: 








l. Instrumentation--The essential instrumentation components have been 
designed or procured. The major unit, the rainfall simulotor-gomma 
probe-infiltromecter is in complete operational condition. This unit, with 

2 combination of needle sizes, and cir and water pressures, is capable 

of duplicating the median drop sizes of netural rainfall up to intensities of 
over 4 inches/hour. In Figure 1, the drop diameter of simulated rainfall 
ot intensities of 0.5 and 4.0 inches/hour are compared with those of natural 
‘rainfall, The control of drop sizes and intensities also make it possible 
to produce a desired rainfall energy even though terminal velocities are 

not reached, - 
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Each of the data sets obtained on density by the PHD system at three 
different temperatures, Table 1, were found to be normally distributed 
when plotted on arithmetic probability paper. In addition, an analysis 
of variance was made to test significance between the average of the 
three data sets by use of the F-test. The ratio of the variance between 
samples, 6315, to the variance within samples, 7432, is 1.18. For 2 
degrees of freedom for the smaller variance and 57 degrees for. the 
larger variance, the value of F for a probability of 95 percent is 19. 47, 
Since the variance ratio, 1.18, is not significantly larger than 1,. the 
difference between the means of the samples is not significant. 


The differences in variance between the 3 samples is, ligewisg, not 
Significant since the variance ratio is small; iie., F = g4/oa, = 1.56 
(Table 1). Using 19 degrees of freedom for each sample, the F value 
for a probability of 95 percent is 2.14. It is concluded, therefore, that 
the observed counts (density) obtained by the density vase, utilizing the 
PHD systerm is not significantly effected over the tempersture range of 
S7letemsOe 2. 


As a result of the analysis of variance, the 3 samples of data (Table 1) 
were combined for determining the random errors of measurements for 
l-minute counts. An arithmetic probability plotting of the combined data 4 
is shown in figure 2, The standard error, S.E, (standard deviation o), 

is 85 with SD = 0.83 percent. The probable error, ©.%. (50 percent 
error), obtained from the Si’ E, is 54 with P.E, = 0.52 percent, Using 
the actual values, the average deviation, A.D., is 70. 


For the 80-percent increase in saturation, the counts decreased from 
8675 to 7120 (mean of 7738), with an average deviation, £.D., of 41 
from the curve of equal deviations (Sigure 3). A change of 1-percent 
saturation represents 51.8 counts. The A.D, of the counts is 0.5 , 
percent, corresponding to the PP, 4. of 0.52 percent for the first test. 
The probable error iw measured saturation for the saturation range 
from 55 to €5 percent is 0.63 percent or (51.8 »¥ 0,005)/41 = 0,.0063. 


Data obtained on saturation by the rainfall simulator-gamma prone- 
infiltrometer and by the tensiometer-transducer system indicates that 
procurement of such data can be reliably obtained in the field, Satura- 
tion and tension data will be obtained in the laboratory from large, un 
disturbed soil cores (5 feet in diameter and 5 feet deep) before initi- 
ating field data procurement this comings summer, 
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2. Laboratory Cata--Results obtained in the initial laboratory tests 

at Oregon State University (Figures 4 and 5 and Table 2) on the 2-inch 
diameter, disturbed, Summit soil samples demonstrate (Figures 8 and 9) 
that the Brooxs-Corey Equations did not describe very well the soil's 
partially saturated behavior. This fact lec to the modification of the 
Burdine Theory for application to the imbibition phase, 


The test on the 14-inch diameter disturbec soil core was conducted for 
three purposes. The first objective, as mentioned earlier, was. to try 
to duplicate the results that were obtainec with the 2-inch diameter 
disturbed soil core, since a different laboratory and different equipment 
were to be utilized. The second objective wes to provide cata to test 
Jeppson's transient, one-dimensional vertical moisture-movement model, 
The third objective was to develop testing procedures and to determine 
whether or not modifications to the equipment were necessary to obtain 
adequate laboratory and field data. The latter objective was necessary 


. since the Reynolds Creek Rainfall Simulator Laboratory was a new facil- 


ity, and some of the electronic instrumentation was new, 


The tests that were conducted on the 60-inch diameter disturbed soil 
sample were necessary to provide data to test the transient, three- 
dimensional axisymmetric flow model and to further perfect testing pro- 
cedures, Two major problems associated with the data obtained from 
the first two tests appear to be solved. Excessive scattering of the 
capillary pressure data was apparently caused by excessive switching of 
the Scanivalve. Poor agreement between saturation-time curves for the 
gamma probe and tensiometers was apparently caused by a discrepancy 
between the elevations of the probe and the tensiometer, 


In the interest of space, only the more pertinent of the laboratory data 
has been presented and discussed here. <Unta not presented here include: 
(a) data obtained from the 14-inch, one-cimensional test, (b) vertical 
spreading data from the first axisymmetric test, and (c) data obtained 
from the second axisymmetric test. | 





3, Mathematical Model Data--As mentioned earlier, the Erooxs-Corey 
Equations did not describe very well the behavior of the saturation- 
pressure and hydraulic conductivity-pressure trends observed for the 
Summit soil, A xiodification to the Purdine Theory provides a function- 
al relationship for relative hydraulic conductivity in terms cf pressure- 
saturation data that gives reasonably close agreement to the laboratory 
data, The results from the Burdine Theory are given by the following 
equation: 
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where K_ is relative conductivity, p is copies pressure, and - 

is the effective saturation defined by St 5 (G5 pia S), in which® 

S. is the residual saturation, The la Stter quantity ig physically the 
saturation at which moisture movement stops, but it is generally taken 
to be the value that gives as good a fit as possible to a functional re- 


oc 


lationship between oo and p. 


To date, the Burdine integrals, Equation (1), have been applied to the 
case of desaturation (drainage). For problems dealing with imbibition 
(wetting), Equation (1) must be modified, since capillary pressure p 
becomes zero for values of S slightly less than unity. This would re- 
sult in a division by zero. Two modifications have been introduced 

to the Burdine integrals. First, a constant pressure p has been 
added to the capillary pressure p; and second, the upper limit of the 
integral in the denominator has been changed to 8 , a value that:has a * 
magnitude slightly less than unity and at which Pathe capillary pres- 
sure p becomes zero. With these modifications the Burdine integrals 
of Equation (1) become 


ee as 
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2: 
The introduction of Sand Po results in additional parameters that are 


needed to describe the hydratilic properties of a soil. A value of 5 
can be obtained from saturation-pressure data when p = o, Further’ 
studies are needed to relate the value of s_ to the physicecl properties 
or to some measurable hydraulic characteristic of the soil, In the 
absence of such studies, the value of Po must be based on judgment 
that is guidec by values obtained by trial for soil for which hydraulic 


conductivity-capillary pressure data are also available. 
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, sf From analysis of the limited amount of data from the Summit soil, it 
It appears that volues for the relative hydraulic conductivity, K , are 

not highly sensitive to small changes in the value of P, particularly 

in the region in which K_ is not too much less than unity. It is be- 
lieved that a reasonable estimate of p_ will be adequate for many 
applications of Zquation (2) for two reasons. First, little flux move- 
ment exists in regions in which K.is very small, and second ,lack of 
sensitivity exists in regions where K.. approaches unity. Consequently, 
hydraulic conductivity-capillary pressure dota will not be necessary, 

but values of saturated hydraulic conductivity will be required. 


The saturation-pressure data (Table 2 and Figure 8) for the Summit 
soil were uscd in a numerical evaluation of Equation (2). Several 
values of p.| were used, and the results were compared with the con- 
ductivity-pressure data for the same soil. Laboratory data indicated 
that a value of a = 0.939 should be used for the Summit soil, but 
several values were tested in the mathematical model. Results of the 
evaluation indicated that values for K.. increased with larger values of 
p over the entire range of capillary “pressures. Increases in K were 
more pronounced in the regions of large ncgative values of capillary 
ms pressure, p; whereas, decreases in the velue of S increased K more 
) es markedly in the portion of the curve whcre p approaches zero. ~ Rela- 
tively good agrccment exists between the values of K_ determined by 
Equation (2) (with S_ = 0.939 and p = 1.0 foot) and the experimental 
data (Figure 9). erefore, Equation (2) has been implemented in the 
mathematical models of partially saturated,transient, one-dimensional 
vertical flow, and three-dimensional axisymmetric flow in porous media, 


























The mathematical problem of one-dimensional vertical moisture move- 
ment through soils has been solved for this study by finite differences 
using the Crank-Nicholson method. In this method the differences at 
the advanced time step are weighted equally with those at the current 

time step, ‘This technique leads to an implicit method that is stable 

for all incremental time steps, and, as such, requires the solution of 
a tridiagonal, cocfficient matrix to advance each time step. 


An interesting result of the one-dimensional study on the Summit soil 
sample is that the hydraulic gradient nenr the soil surface is much 
greater than the unit gradient that exists for saturated flow in a ver- 
tical column. The magnitude of the gradients above unityin the par- 
tially saturated flow are the result of capillary forces. 


A number of one-dimensional solutions was obtained by specifying sev- 
1) fa cral application rates and initial values of hydraulic head, h_, to 

correspond to specifications used in obtaining solutions to similarly 
formulated, threc-dimensional axisymmetric problems of flow from 
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circular infiltromcters. The saturation at the surface centerline of a 
the three dimensional axisymmetric problems is greater than anywhere 
else at the soil surface, as illustrated in Figures 10 and ll. In fact, 
saturation approaches initial saturation at some distance from the infil- 
trometer ring, but it also decreases within the ring. 


For all problems at all time steps, the magnitude of the saturation for 
for the three-dimensional axisymmetric case is less than that for the 
equivalent one-dimensional problem, The decrease in saturation for 

the three-dimensional case can be attributed to the radial component of 
velocity removing some water, even at the centerline. The radial com- 
ponent of velocity appears less significant for lower rates of applica- 
tion, particularly if a low rate of application is specified in conjunction 
with a relatively large value (small in absolute value) of initial hydraulic 
head i 7 


The influence of the radial component of velocity (or the spreading ef- 
fect) is different, cepending upon the point being considered. The satu- 
ration at the surfacc near the infiltrometer ring will be less than at the’ 
centerline. Table 3 gives a comparison of saturation at the soil sur- 
face from one-dimensional solutions with those at the surface centerline 
and at He infiltvemcter. ring from equivalently specified three-dimen- 
sional axisymmetric problems, 

Also of interest is the manner in which the saturation increases at a 
point within the soil profile as a function of the application rate and of 
the initial condition, The point selected for illustrative purposes is 

the soil surface for the one-dimensional problem. Values for the ratio 
of saturation at the surface divided by the initial saturation, also on the 
surface, have been plotted (not shown here) against the application rate 
for several time steps. On these figures different values of initial 
hydraulic head define separate curves. These figures illustrate (a) 

how the saturation on the surface increases with smaller values (large 
in absolute value) of initial hydraulic head and (b) how the surface satu- 
ration increases with larger rates of application. 


Future plans concerning the laboratory work include extracting a 14- 
indi diameter core from the 6@sanch diameter disturbed Summit soil 
sample that was last tested. The 14-inch core will be tested under 
one-dimensional vertical moisture movement conditions to obtain hydrau- 
lic conductivity-capillary pressure data. These data will be used to : 
test further the modified Burdine Theory that was mentioned earlier. 

At the time the 14-inch core is extracted from the larger sample, 
gravimetric samples will be taken to determine the exact bulk densi- 
¥#es and porosities at various depths. 
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TABLE 3.-- Comparison of saturation at the soil surface from one-dimensional solutions with those 


at the surface centerline and at the infjitrometer rin specified 


three-dimensional axisymmetric problems, 
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Initial Applic. Time Ro = One-Dim, R, = One-Dim. sat./ Percent 
head rate Parameter sat./axisym. sat. axisym. sat. at diff. 
ho a] 1K, 5 at centerline infiltrometer ring (r=L1") R1-Ro 
x 100 
O 


a a 


-0.859 0. 0954. .99 1.0081 1.0966 9.5 
~2,859 0.0636 1.0 1.0275 1. 1379 10.6 
~2.859 0.0477 1.0 1.0433 1.1567 10.7 
-3.859 0. 0636 RO 1.0285 1.1396 10.3 
-1.859 0.0636 1.0 1.0256 1. 1304 10.2 
+0. 141 0.0636 1.0 1.0195 1.0792 5.9 
+0.141 0.0954 1.0 1 €150 1.0898 7.4 
+0. 141 0.159 £75 1. U446 1.0900 4.3 
-0.859 0.159 5475 1.0506 1.1118 5.8 
-3.859 0.159 5 1.0293 1.1299 9.8 
+0, 141 ah 2226 3 1.0756 1.1461 6.6 
~1.859 0.159 585 1.0415 1.1157 Ti 
-1.859 0. 2226 470 1.0729 1.2041 12.1 


qf Taken from Jeppson, R. W. 1970. Solution to transient vertical moisture movement based upon 
saturation-capillary pressure data and a modified Burdine Theory. Progress report PRWG-59C-5. 
Utah Water Research Laboratory, Utah State University, Logan. 
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Following the test for conductivity-pressure data on the 14-inch core, 
tests will be conducted on a series of 32-inch diameter, undistributed 
soil cores taken from several different locations on the Reynolds Creek 
Experimental v/atershed that contain different soil types. Four of 
these cores have already been extracted. Plans call for obtaining two 
more of the large, 32-inch, undisturbed soil cores. Saturation-pres- 
sure data will be obtained from each core in the Reynolds Creek Rain- 
fall Simulator Laboratory by using the rainfall simulator-gamma probe- 
infiltrometer and the tensiometer-transducer system. Since the 
modified Burdine Theory appears to give adequate values of relative 
hydraulic conductivity, present plans do not call for obtaining conduc- 
tivity-pressure data from the undisturbed soil cores, Data on satura- 
ted hydraulic conductivity will be required, 


After data have been obtained from each undisturbed soil core, the 
Sample will be broken up, gravimetric samples will be taken, and a 
disturbed sample will be formed. Saturation-pressure data will then 

be obtained for each disturbed sample. Comparisons between the data 
from the undisturbed and disturbed samples can 2 made, in addition 

to further analyses of the hydraulic properties .: -.ae different soil 

types. An attempt will be made to determine p...ametric equations that 
adequately describe the hydraulic properties, especially the saturation- 
pressure characteristics, of the Reynolds Creek soils. r 4 


Following the laboratory studies described above, further use of the 
mathematical models will be made. Parameter optimization will be 
incorporated into the models to determine: (1) the parameters in the 
saturation-pressure relationships that evolve from analyses of the lab- 
oratory data, and (2) the saturated hydraulic conductivity, The labora- 
tory data. will provide a check on the parameter-optimization results. 

If parameter optimization is successful, then less field data will be 
necessary in future studies. 


Future plans concerning field work include in situ tests using the 
rainfall simulator and its associated gamma probe at various locations 
that include different soil types on the Reynolds Creek Experimental 
Watershed. The ficld data will be used to satisfy three objectives, 
First, a field verification of the transient, three-dimensional axisym- 
metric, partially saturated flow model is needed. Second, tensiometers 
will be installed at some of the test sites to provide field data to fur- 
ther verify the use of the modified Burdine Theory and parameter opti- 
mization. Third, for different initial soil moisture levels, ficld tests — 
will be conducted to determine the potential quantity of water retain- 
able by various soil-vegetation complexes, and the potential guantity 
that is independent of infiltration, 
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Publications completed and in preparation since last year's report 
are the following: 


Jeppson, R.wW. 1970. Transient flow of water from infiltrometers-- 
formulation of mathematical model and preliminary numerical 
solutions and analyses of results, Progress report PRWG-59c-2, 
Utah Water Research Laboratory, Utah State University, Logan, 


Jeppson, R.v/. 1970. Formulation and solution of transient flow 
of water from an infiltrometer using the Kirchhoff Transformation. 
Progress report PRWG-59c-3, Utah ‘ater Research Laboratory, 
Utah State University, Logan, 


Jeppson, R.vW. 1970, Determination of hydraulic conductivity-- 
capillary pressure relationship from saturation - capillary pressure 
data from soils, Progress report PRv/G-59c-4, Utah v/ater 
Research Laboratory, Utah State University, Logan. 


Jeppson, R.v/. 1970. Solution to transient vertical moisture movement 
based upon saturation-capillary pressure data and a modified Burdine 
Theory. Progress report PRWG-59c-5. Utah Water Research 
Laboratory, Utah State University, Logan, 


Wei, C-Y., and Jeppson, RoW. In Press. Finite difference solu- 
tions of axisyrmmectric infiltration through partially saturated “Se 
porous media. Progress report PR./G-59c-6. Utah Water Res 
search Laboratory, Utah State University, Logan. 


Penton; V.&. , andjHamon, W.R. In preparationf~l97l. A rainfall 
Simulator, gamma probe, infiltrometer, To be submitted to 
Water Resources Research, 


Jeppson, R.wW,, Schreiber, ®.L., and Biggs, E, N. In preparation, 


1971. Verification of a mathematical model of transient water flow 
from infiltrometers, Water Resources Research, 
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Figure 3. Average gamma probe counts per minute (foure 
minute readings) during passage of wetting front 
at four inches below soil surfacee 
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Figure 4. Laboratory saturation-capillary pressure curves for Summit soil. 
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Figure 6. Observed effect of lateral spreading from a‘circular 
infiltrometer at a depth of 3 inches in a disturbed 
Summit soil sample. 
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Observed time variations of capillary pressure and saturation at 
a depth of 4 inches in a disturbed Summit soil sample. 
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Figure 8. Comparison of observed saturation-pressure data 
for Summit soil with that predicted by the 
Brooks-Corey equation. 


a a nr pie 
as 7 


RELATIVE HYDRAULIC CONDUCTIVITY, K, 


003 
/ 


[tem BROOKS - COREY EQ. 


/ 


/ | 
/ 
Fata BURDINE INTEGRAL | = 
6 


/ 






0.1! 02 0.30.4 06 0810 20 3040 608010.0 20.0 300 
RATIO OF BUBBLING PRESSURE TO CAPILLARY PRESSURE, (p/p) 


Figure 9. Comparison of observed hydraulic conductivity- 
capillary pressure data with that obtained by 
the Brooks-Corey equation and that obtained a 
from saturation-capillary pressure data through 
a modified Burdine integral. 
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Figure ll. Predicted iso-saturation lines at - step f= 2.085 for flow from a circular infiltre ‘er. 
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CRIS Work Unit No.: SWC-O11-fBo-l1 Code No. Ida-Bo-106,.1 


Title: Natural evaporation from sagebrush range- 
lands, alfalfa, and stock ponds in a semiarid 
environmente 


Location: Northwest Watershed Research Center, 
Boise, Idaho. 


Cooperation: The University of Idaho will furnish a port- 
able micrometeorological instrument trailer 
for initial acquisition of data at study 
sites in the Reynolds Creek Experimental 
Watershed. 


Personnel: L. M. Cox and W. R. Hamon;: G. H. Belt, 
University of Idaho, School of Forestry 


Date of Initiation: Novemer,1968. 


Expected Termination: November, 1973 — 


Objectives: 


1. To determine the evaporative loss of water from sage= 
brush rangelands, irrigated alfalfa, and stock ponds 
while observing pertinent meteorological parameters and 
the soil moisture statuSe 


2. For predictive purposes to develop relationships for 
associating the evaporative loss with meteorological param: 
eters, type and degree of surface cover, soil moisture, 
and potential evaporative demande 


Need for Study: 


The water and attendant sediment yielded from a watershed or 
land area represents a balance between that entering the area 
as precipitation, that stored, and that lost through the proc-= 
esses of evaporation and transpiration, For semiarid, sage= 
brush rangelands, as represented by 50 million acres in the 
Northwest, the percentage of available water lost to evapotran= 
spiration is around 95 percent where the yearly precipitation 
is near 10 inches. At higher elevations with annual precipi- 
tation of 30-35 inches, and where the major portion of precipi-~ 
tation occurs as snow, as little as 50 percent of the total ; 
available water is lost to evapotranspiratione 


In the extremely low precipitation areas there is need to con- 
serve all the moisture and provide as dense a cover as possible 
to prevent excessive erosion when convective rainstorms OCCULre 
The August storms of 1968, produced 2 to 3 inches of precipita- 
tion and sediment yields of 2 to 5 tons per acre on such lands. 
In higher precipitation areas, a potential exists for improving 
Management and increasing forage yields for a greater site 
utilization of available watere 


LO-2 - 


A predictive capability for the runoff and sediment that ‘can 

be expected from rangelands is required by the Soil Conservation 
Service, the Bureau of Land iianagement, and land ovmers to Cary] 
out multiple-purpose management of these areas, and to opt cimize\ 
forage production, water yields, sediment reduction, wildlife 
protection, ana recreational potential. A complete understand-~ 
ing of the evaporative process is essential for developing pre- 
dictive relationships for the evapotranspiration component for 
definable soil-vegetation complexes under a particular level of 
management » 


A capability of measuring and predicting the water lost by evapo- 
transpiration is required! to“apply an infiltration model in the 
prediction of the surface runoff component and resulting sediment 
yield from any land area. Improved management, based on a know- 
ledge of available water and the physical factors controlling 

the amount of water lost to evapotrans pir agtoas can reduce the 
sediment loss by possibly 50 percent in the S- to 12-inch pre- 
cipitation zone and ey increased forage yields by as much 
as 80 percent in the 20 to 30-inch precipitation zone. Improved 
efficiency in the design of water retention structures will. 
result from improved estimates of runofxs and sediment since the i 
estimates of soil-water storage will be improved with better in- 
formation on evapotranspiration losses, 


Design o Experiment and Procedures to be Followed: 





The evaporation and evapotranspirat ion studies wil 
in the outdoor hydrology laporatory of tne Noexrthie 
Research Center-—Reynolds Creek Experimental Va 
dures for measuring and predicting evapotranspi 
with sparse vegetation cover, and unsaturated s 


De cone aa 
*t Watershed 's, 
sae Rela 
LOi 


a 


The primary measurements to be optained in the field twill be 
those concernin ng the various parameters contained in combination 
formulas for estimating evapotranspiration, (See 1968 Annual 
HepOrt.) In jaccition, proliles ox wind and temperature along 
with humidity at tio epee for Bovien ratio calculations will be 
obtained, Specifically, the air temperature will be sensed at 
elevation z, where the logarithmic wind profile extrapolates to 
V = oe . 


e 


soil~heat’ storage; ‘soil heat | flux; “soil moisture; type, height, 

and cover of vegetation; soil oa and percentage of rock pave- 
mente The vegetation and soil data will be used for designating 
definable soil-vegetation sane The evaporative >serformance , 
of these units as related to exsisting meteorological conditions - 
will be used for snecifying a predictive model, 


The required meteorolocical data will be complemented by data on 
1 


The energy balance and Sovwen ratio procedure using AT and Ae at 

two heights willbe useditov obtain a e measurements of 

independent meas- | 
;O11 moisture » 





evapotranspiration where applicaple. Also, 
urements of evapotranspiration will >be made ¥rom 
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determinations and =ron lysineters Soil moisture data till be 
odptained by neutron and gemma Baneoes from access tube networks ! 
on definable soil-vege tati ion complesres, Lysineters of the hy- 
draulic-vieighing type have been installed at several study sites, 
These lysimeters consist of an inner tank (5 feet in diameter 
and 6 Teer ceed ) oOo Suntastiroecatsorlithatiresis-omsa coil of i2- 
inch eutyl tuping 2ilied: jack Liquid end oe sed in en outer tank, 
The putyl tubing is connected to br sure trans= 
ducex for monitoring. pressure change 


pres 
eaten 770m the loss or: 


Ade ta on or “water's 


rrelininary data on teroerature, erat: cet anc tiie profiles and; 
on energy com>sonents have been obtained as a cooperative endeavor 
with the Department of Forestry, University of Idaho. This co- 
operation will be continued to test measurement apparatus and 
Sercors, entero Oprain Geta 2or evaluating transter coefficients 
for dry, sparsely vegetated’ areas’. 


Mroecranentary Data end -Opservat2 ons : 


Lysinmeter data were collecte 
Site, and they tilere of margine 
COLrdi ng or Oo lens. 


ad only fron the Sheep Creek Study 
Pua lrcy pecause “Ol pressure tTe- 


fhe “yseneters Gn the Grrioated “altalta: field vere ‘completely 
enguizteu Sy a high vwater table*"because or a Teak ‘in the’ outer 
instrument tank, 
The lysimeters at Reynolds Mountain required servicing rollow- 
ing the winter, but late season snovdrifts made it impossisle 
<0 move a crane to the site until late sumer, Consequently, 
the high evanorative losses that occur dering the cpring season, 
when soil moisture is available, vere not measured. 

Son Ol ti rere te 
Shi morsture cate colrecti or, 6 tone T2ith brecipitation and) run- 
off, wac continued for the networks in the Sumut, Sheep Creek, 
and the 2eynolds Mountain study basins, Soil moisture data 
collection was concluded in Cummings Field (Ground-iater Study 
Basin). Gnforsation on these networks is: contained in, Taobe’ 2, 
p. 1-5.) Additisnal > “seil moisture’ sites vere installed curing 
1970 for use in studies under tesearch Outlines Ida-50-105,4 


and 105.66 


All soil moisture data have been placed on cards for cata 
hancling and analysise 


Energy Dudget Measurements 


TLih Ge OG, Lk Sas eee Gre ee 1+ of the University 


During. che. la SD 3e 
ed energy budget and meteorological data at three 


te 
of Idaho, optaine 
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Sites in the Reynolds Creek tlatershed as part of a cooperative. 
study. Comments on these shea were made in the 1969 Annual A 
Report. the data were used during the past year in an analysis 


of evapotransdiration from sagebrush rangeland. 





ee 


Comments, Interpretations ANG «BIDtUTe Ui Lanes 
Lysimeters 

Temperature induced pressure fluctuations are still a 

1 


in the pressure recording systems used fox the soil S 
eekly data is about the best these systems are good for, and 


Le 3 
then only curing periods when evaporative losses are highe Con- 
tinuous pressure readings may be imsossible WLthOuUT AeC.ePOTICrs 


LS currently underuay 2®or testing 
using a d.c. pressure ae poe curre 
research people e Consequently, the re 
changed to weekly manometric measurements 
completed. 


iD 

ed by snow 
ms will be 
Chas OTC 1s 


SOlt, fOr Se cime 


Soil moisture cata, collected. since1964on.+ 
Creek, and Reynolds liountain study-pbasinsewi pe analyzed durin 
1971 in connection with studies under “esearch Outline 105,4, 
the soil moisture data obtained as part ox the study uncer . 
Research Cutline Ida-5o0-103.3 have seen analyzed in assessing 
the water valance of an irrigated ficld.. (Reter to this latter 
Cutline for comment). 


e Summit, Sheep 
. : 
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watershed, Evapotransniration 


rrecipitation data obtained *rom cual precipitation gage sites 
(shiclded and unshielded gages) have seen used in a computational 
model to optain realistic values of atmosnheric precizitation. 
(See Research Cutline Ida-30-102.1). Svapotranspiration, EeTe, 
estimates for the water year of 1669-70 have been obtained for 
4 sub-watersheds in the Reynolds Creek jjatershed by use of the 
Water salance equation, with storage changes, *S, neglected, as 
MOTCC 1713 eI) 1 Gs elt 

WEOM, Tacle oly the Sat O. Of redo, oO. COLMr Ged precipitation, exe 
pressed as a percentage,,. ranged =rom 190% Ln she. Sumit, basin 
to 67% in the higher elevation area o£ the Tollgate “Jatershed, 
Also, Sor the sume areas, the evapotranspiration, e:mressed as 
a percentage of "potential" evapotranspiration, ranged from 27 
tO S7%e 

As the above indicates, water balance computation can be made 
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Table 1. “Jatershed evanotranspiration for the trater 
year 1969-70. 
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Summit 10.00 0 O 10.00 af 
(4180-4800 ) 

Loviexr Sheep 14.35 0.03 O 14,85 32 
(5200-54258) 

Saimon 13035 3.45 O 15.40 
(3675-62C0) 

Tollgate 30.65 10.14 O 20.51 29 
(4600-7300) é 
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Evapotransoirat tvonwkcTimates 


The Sicld data obtained by micro-meteorological instrumentation, 
asta cooperative etlor’ vith the University of Idaho, in the 
late spring of 1969 have becn analyzed to obtain radiation 
balance estimates, and to commute evapotranspiration by the 
Bowen ratio method and by energy balance equations. 


At the low sagebrush Sheep Creek site, under clear 


skies, the 
surface elbedo was observed to be relatively constant at 15.0 


t 
+ 1 percent. Net radia ation,, Rev vie Sine GUI MO Deal linea 
function (rt = 0.99) of incoming solar radiation, Giit, and can 
be expressed as: 


de Erie Wed at Ss ee 6 er 2 ba a (1) 


where R, and Siil are Giver 1 Lye/iiia. f Gm cloudy days a similar 
relations hii exists, but there is more scatter in the cata. 

Even so, R, can be estimated with reasonable accuracy (+ 0.05 
Lyse) by uSe of equation 1. A close correlation was also found 
between the soil heat #lux: and, net radiation. (Sample data are 
contained in a Research Report—), 
Ne 


1/ P.E.T. represents potential evapotranspiration. Values 
extracted from 1968 Annual Xeport. 


of SoeLe fe GS eae uid Bed Sati aspiration from Low Sagebrush 

a Range in Southern Ideho, Research Project Technical Comple- 
TL tcl te. .-cer rae ee Research Institute, University 
of Ideho, Moscow, Idaho (October, 1970) 


1D 


The relationship of evanotranspiration, rag to net Tediationg 
4? {Or the low sagebrush site where water is limiting ene idee 
the regression equation 


5 


LE = 0.342 R, + 0.14 (2) 


with an r value crfs0Os//wand stancarc Gtroreo" estimate or 22 
Ly./Day or aoprouzimately a 20-30 percent error, 


RATCS OL VGVAPOLTANSULSCLOM, dist esscOIm 2m aver age daily dbasis, 
at the Sheep Creek site during June, were estimated using the 


energy balance-Bowen ratio procedure. £Verage daily Dele rates 


ranged from O, OF 4 to O.12 inch under differing conditions of soil 
moisture and radiant SUPE GY supplye aaa energy xlixc accounted 
for 28-46 eres ox net radiation during daylight Hours. SOLe 


heat flux was relatively eee aera 11-15 oercent of net radia- 
tion. On all but one of the 6 days of observations, more energy 
was partitioned into sensisle heat than latent heat indicating 
that soil moisture availability, not energy availability, was 
the: jlrmmriauggiactor.ductmg this period, Suregvac true even 
though significant amounts of precipitation and overcast con-= 
ditions occurred during the period of measurement. 


Hourly variations in computed £.T. were the result o£ changes in & 
the Bowen. ratio., The magnitude .of the ratio is determined “by 
both vapor pressure end temperature gradients, the latter being, 
the more variable. Variation in temperature gradients are co 
lated with mean tind velocity during overcast conditions, but dre 
primarily determined by the net radiation flux during clear skieS., 





Comparison of Methods for Estimating 3/ 
Evanotranspiration on Semi-arid Lands in Idaho= 


Micrometeorological data obtained rom several sites in 
the Reynolds Creek E:qperimental “jatershed have been used to 
evaluate various ‘evapotranspiration (E.T.) formulas for 
applicability in the dry-land situation where non-saturated 





3/ , Computations. and analysis oy Do, George elt.) collece so, 
Forestry, University ot Idaho, 
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ce Soils dominate. Bowen ratio, eneray balance-aerodynanic, 

Ma? and combination equations as presented by Tanner and Fuchs4/ 
and Fuchs, ete al.2>/ have been used to obtain estimates of 
EeT.e These equations are given below: : 


Bowen ratio 


E — (Kay -&) 
i } +x he) z 


( ey - &) 


Eeg = —(Rn- 4) +2 cp h (7, - Tz) 
© Es = —[Sksee)]} (2-6) 
+(p cps) h [Ke - €2) -(ef- ¢,)| 


POT e Lelie 


E,= -[s/cs+n][ ext) +(ptp/s)h(e-e)] (4) 


Gene. Combination 


E, =-[(r+5)/s] E, + (pep/s) h Cs eae 





Tanner, Ce. £&. and Fuchs, M. 1968. [Evaporation from un-= 
Saturated surfaces: A generalized combination method, 
Jour. of Geophys. Res. Vol. 73, No. 4. (February 15). 


ie 5/ PucCHS, 7.7 -Tanner,” Co 8.4” Thurtell,- G. We, and Glock, T. A. 
Ly 1969. Evaporation from drying surfaces by the combina- 
tion method, Agronomy Journal (Jan.-Feb.) 
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{ 
More detailed eszsplanations of these formulas tith list of (6 
symbols can be found in the referenced papers. Table 1 summar- 
izes the data requirement of the See equations (equations 
1, 42;.3, @nd 5) atid--he-equeticn Dor poten: al evapotranspira- 
tion, equation 4, 


Tale 2 
liethod equation Data Required 
Val Va[*o|*a]2s|Poltaltel® | Bale | 





bowen ravio | (0 Sen7a05505% 
| (2 ) aif, 

S-l4 Cortoination | (3) ee es . 

Potentialie.2% (9 Se ed cca 

eno . Combing tz a | (2) | | Pt fae fe | fae tae fae Toe | 


The exchange coefficient, h, is je:;ined as 


h=k* v (p+ In (2 + ope “g 


Energy balance 


\A Wee 
ia) Wee 








where D=£(R;) and R; = £(T), Tz: Vas ve (5a) 


Primary differences in the da requirements [or the estimating 


at 
equations in Teaple 2.can >be summarized as foliovis: 


pressure 
1. The Bowen ratio methcd requires measurement o£ vapor/at a 
second ievel, but MmOv ai evaluation On i. 


2e The energy balance = aerodyi 
<9 


2 
evaluation ox h, and surzace <— 
ment of vapor pressure. 


ic method requires Numer CAL 
eir-erature Dut no measure= 


£velyation of h, 
3e The Slatyer-Mcilroy combination method requiresf®surface 
eS coy ana surface vapor pressure. 


~~ 


4, the generalized combination method sequires a numerical 


evaluation of h end surface temoerature. 


The energy balance - aerodynamic and combination formulas re- 
quire. evaluations ox, bye by, usclolmthesTeves fumedion~ Thais 
uires evaluation o£ teimerature and velocity 
gréedients, and calculation o¢ 
study, T, was obtained by statist 
erature proxvile by assuming simil 





he Richardson Numper. In this 
cally fitting the air temp- 
KLiviwith the wine prorne. 


ere nee ee : Te ry en repens 


aay 


sae ey Pinscriiin nsasasale Roatan, Derek ans © eR T oer seen eee ian 


» Was ootained at height z 
Ss to zero = ngt ground 
erveanduPuchs—/ i. 


By this method, the temperature, ZG 
where the wind velocity extrapolate 
surzace temperature as used by Tann 


Oo 


ae 6 : on . - 
inpeigures S/ 1 2; and) /ph.by eStimatesoptained.py athe enercy 
palance and comoinaticn equations from an irrigated alfalfa field, 


a ‘big. sageprush isite, and .a llovmsagebisush wite,are compared with 
those obtained wy the Dowen ratio mG iO wa era plis tor che 

irrigated Blratva =1elc contain .caga,onrained.oniy under 1 
conditions. the fever data points from the irrigated site 


—— =— 


re 


result from _using,.data free of advective heat, 


Differences between estimates appear to be the resuit of 


iffer ae C 
YVALLALLONS.in data.required.oy.the several .cormulas aren with 
sperimental error. A Significant »oint concerning the liriga~ 


ted alfalia data. is the. good agreement of the several Benin are 
methods, the principal source of experimental error was in the 


measurement o£ the vapor pressure gradients. In the case of 
the irrigated .tileld, .vanor PEessuse gradients were substantially 
greater (approx. O.2mb. ) than atesehe.other quo sites., It is : 
tentatively concluded that the more accurate measurement or 

I 


%G 
vapor pressure resulted in the good agreements 
At Doth sagebrush sites, Buchices 2e4anG nas the -fouen ~raAti ro 
cedure yielded estimates of E.T. smallex in magnitude than the 
other methods, Greater scatter of noints at these dry sites i 
also obvious. £8 stated above, this is proba>sly due to the 
lesser accuracy in measuring vapor pressure gradientSe 





2yYOLg .1i. toe .,;xcbange coer sicient, hh, are a potential source 
oz -bias .explaining..departure .or ZT estimates .<zom.the 1:1 
line in Figures 2 and 3. Comparison of exchange coefficients 
can be made by use of the dimensionless coefxicient K,* 
defined as: 
Te eK Kp 
eS te ~ ( 6) 
Au 
——n fo = 
crn ) 
. ¥ T = j o~e 
where Kp = He(Ze Z,) 





ot erate lies oz Ky, 
h the ced ae 
imibear coefficients 


and, H -is -cadcudkated toy the, .orren r 
were used to compare data from thi ¥ 
form o£ the MEYPS exchance coefsicient and 





6/ Figures follow p. 10-11 
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obtained by Pruitt and Aston from the Davis lysimeter, These 
data plotted in “igure 4, although exhibiting considerable 
scatter, sugges hat the exchange coefficients oitained from 
thee LEYPS fanct and used in equetions 2, 3,° and > are or 
reasonable magn « The lack of any systematic bias suggests 
that the coefficients are not responsible for the systematic 
differences between Sowen ratio &.T. estimates and estimates 
by equations using the transfer coefficient for the sparsely 
vegetated rangeland sites, 
The estimating equations were evaluated assuming SQua laity, or 
exchange coefficients, ~Equation assumes Ig; = Ki and 2, 3, 
and 5 assume Kj; = Kijze While Reedy aaa that othéx assumptions 
may be appropriate, particularly at the more unstable Richardson 
Number, Ry, nei cher the emperimental results obtained nor the 
literature provide a Satisfactory alternative assumption. Since 
the “lc l-relationshtprtor the “irrroatred ssire aé significantly 
meena than those of the dry rangeland sites and since the 
solute magnitude of the exepatice corr Peient appears appro-= 
eoteS abies appears that the magnitude of the vapor pressure 
gradient may have been Gide raat: pa een Under-estination of 
the vapor pressure gradient and/or inequality or exchange : 
coefficients (Ky = Kp E) appear at this writing to be the mor 
ion iS 


propable e:splanat ESE the dizferences in E£.T. Be ee 
obtained. eé 
The computed values of E.T., as obtained by the energy balance 
and combination methods, are essentially identical as demonstra- 
Ted™ in Figures/ 192 "and S.*Othistsiani aoc or Sagan in the 
data-since the consination equations are simply disguised 
versions of the basic enercy balance as represented by equation 


2, Obtained by use of identities contained vanor pressure terms, 
On examination of Table 1, it is readily apparent that the 
combination equations, equations 3 and 5, récutresdata on Vapor 
pressure in addition to the data needed for the energy balance, 
equation 2, Therefore, the only purpose for the Slatyer-lMcilroy 
combination equation is to e:press the energy balance 3n terms 
of surface and air vapor pressure ceficits so that an equation 
for potential evapotranspiration, eQuatLon 2, Cam” pe cert vece 

It follows oy a simple algebraic exercise that the only inde- 
pendent equation for computing evapotransziration is the energy 
balance formula, equation 2, 


i) 


) 





TH LB uate Pare Of) hdlet eos lig Je Atmospheric and surface 
factors affected evapotranspiration,. chaise 7 SSE 69), 
inrooks ah Mae Sere PRSLOSSy investigation of energy 
and mass transfers near the ground including the ine 
fluences of soil-plant- atmosphere system. University of 





California, Davis, Caligvornia, (Final Reports Task: 3A 
99-27-005-08. AD 410-263, 
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Future studies will se directed toward developing a useable 
energy balance equation, based on the res istance concept, 

aif she sparsely veceteted rangelands, Additional field data 
will pe obtained to evaluate the air transport resistance and 
the internal resistance to flow of water cy the soil and 
vegetation. An ch eee will be made to relate the latter re- 
Sistance to vegetation tynes and available soil moisturee 
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irrigated Alfalfa 
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Figure 1. Evapotranspiration computed by different 


methods from selected 30-minute mean 
profile data, 
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Figure 2. Evapotranspiration computed by different 
methods from sclected 30-minute mean 
profile data, 
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Evapotranspiration computed by different 
methods from selected 30-minute mean 


profile data, 


Figure 3, 
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(+) Sheep Creek (Low Sagebrush) 
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(a) Cummings Field (Alfalfa) 
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Figure 4, Dimensionless exchange coefficient, K_, in relation to the Reynolds 
Number, R., as determined from evapotrans diration as obtained from 
’ ’ 3 
lysimeters, the Bowen Ratio Method, and ky use of the KEYPS Function, 
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CRIS Work Unit No. SWC-012-f£Bo-2 Code No. Ida-Bo-103.3 
T101C: Ground-water flow system under an irrigated 

field. 
Location: Northwest Watershed Research Center, Boise, Idaho. 
Personnel: G. R. Stephenson, and J. F. Zuzel, ARS-SWC; 

R. E. Williams, and D. W. Allman, University 

of Idaho. 
Date of Initiation: ‘(Te ee Tey! 
Expected Termination: Originally planned, 1969. F 


Recommendation, September 1971. 


Objectives: 


The ultimate goal is the evaluation of the contribution of ground 
water to streamflow in an area where ground-water flow is influenced 

by irrigation practices. The field study area will consist of an 
irrigated alfalfa field in the Reynolds Creek drainage basin. Specific 
objectives are as follows: 


1. To evaluate the hydraulic properties of the materials unger ly tag 
the irrigated field. . 


2. To delineate the hydrogeologic boundaries of the flow system 
affecting the irrigated area 


3. To evaluate the inflow and outflow of ground water across the 
boundaries of an irrigated field. 


4. To conduct an empirical field analysis of che effect of different, 
specified irrigation practices on ground-water flow. 


5. To construct a resistance-capacitance analog model to depict the 
variations in the ground-water flow system when evapotranspiration, 
irrigation practices, and boundary conditions for the flow system 
are varied. 


Need for Study: 


Evaluation of the hydrologic flow components in a watershed requires 
knowledge of the underground flow system. The water status for any area 
is dependent upon the integrated inflow and outflow. 


A&A common water-use practice in mountainous watersheds is the diversion of 
water from streams for irrigation. The increased consumptive use (evapo- 
transpiration) greatly affects the quantity of water yielded from these 
watersheds. A considerable portion of the diverted water in such irrigated 
areas is returned to the stream as both surface and underground flow. 


EinzZ 


A strict accounting of the watcr ontecring and leaving an irrigated areca 
necessary for development of procedures to predict the subsurface flows é 
and the losses to cvapotranspiration. Such prediction schemes when developed 
can be utilized in digital or analog models to evaluate the variations in 

the ground-water flow system whon cvapotranspiration and irrigation practices 
are altered. 





Design of Experiment and Procedure to be Followed: 


The field phase of the study will be conducted on a 62-acre irrigated field, 
Figure 1, in the Reynolds Creek Experimental Watershed. Instrumentation 
for collection of data will include a network of observation wells and 
piczometers at critical locations, a network of soil moisture access:*holes, 

Parshall flumes to measure input ont output of surface water, a recording 
raingage, a hygrothermograph, a wind recorder, a net radiometer, and a 
Iysimeter for independent evaluation of evapotranspiration. 


Pumping and injection tests will be used to detcrmine aquifer constants, * | 
and several recorders will oc installed to monitor continuous variations 
in ground-water levels. a 


Scmiweekly, weekly or biweekly collections of all data will be made during | 
the year, being most intcnsificd during the irrigation scason. 


Experimental Data and Observations: 


Three years of soil moisture data, 1968-1970, were collected for this oth | 
by use of the neutron probe. These data are used in conjunction with changes | 
in the water table for cvaluation of totai water losses. Measurements for —_ 
this 3-year period were made at 22 installations, 16 of which extended well || 
below the water table. When analyzing these data for cach irrigation season, | 
one phenomena consistently prevailcd--soil moisture at specific levels | 
within the saturated a by varied from site to site, as well as between levels. 
of cach site (Figure 1). A systematic time variation was evident in that 
the values were lowest in June and highest in October. This systematic time 
ariation permitted corrections to be made to yield constant soil water 
eae for depths continually below the water table. 


Several phenomena have been considered in an effort to explain the time- 
dependent soil moisture variation below the water table: 


. Changes in composition of the soil water. 

Variation in the volume of gases trapped below the water table. 
Confining sressurc. 

Exsolution and dissolution of dissolved gascs. 

Changes in the bulk density. 

Equipment crror as a result of temperature scnsitivity. : | 


° e 
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cei RST nace 7 Water loss determinations were completed using corrected 
soil water values and piezometer water level data at 13 sites throughout 
the 62 acre irrigated field. These loss values including specific yield @ | 
and any transpiration losses by the alfalfa, are for a soil column extena 
from a depth of 10 inches to the water table. The acrial distribution of 
these calculated values indicates that the total water loss below 10 inches 








1/ Figure follows page ))-7 . 
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is relatively larger in the alluvial soils as compared to those near the 
edges of the field which are shallow soils develoned from bedded clay 
and silts. 


This method of using the corrected soil moisture data and water level data 
is believed to be the most satisfactory method of computing the total water 
joss in situ. 


The water loss information was used to evaluate daily losses due to ¢vapotran- 


Sspiration at each site in the study area adjacent to a soil moisture access 
tube. Losses in the upper portion of the soil profile, determined by 
gravimetric sampling, were included in the daily cvapotranspiration loss. 
This method proved reasonabic for the conditions that prevail. : 


The cocfficient of horizontal permeability was cvaluated using Hovorslev's 
(1952) method which considers the piczometer screen geometry and basic time 
lag. his method does not permit computation of the ratio of horizontal to 
vertical permeability. 


Comments, Interpretations, and Future Plans: 


A considerable amount of effort was put forth to try to explain the time-_ 
dependent soil moisture variation below the water table. These changes, as - 
illustrated in Figure 1, occur at individual levels at each site as well 

as between sites. Several factors were investigated: 


1. Changes in composition of the soil water. The capture of thermal 


neutrons by such clements as cadmium, boron, lithium or chlorine, 
which may be present in the water, could cause yearly variations in 
the apparent soil moisture content. It is not imown how much the 
concentration of these clements must vary in order to produce a 
change of 10 percent in the neutron count. However, the range of 
variation of the concentration must be reasonably large as there was 
no information found on this subject in the literature. From chemical 
analyses run to date, only chlorine of the above-mentioned elements 
was detected, and only in very minor amounts. At present the 
possibility of composition changes in the soil water does not appear 
to be the principal factor causing a time-dependent soil water change 
below the water table, although it is still a possibility. * 


2 


2. Variation in the volume of gases trapped below the water table. 
A change in temperature of the ground water would cause volumetric 


changes of entrapped gases. However, an annual variation of :4°F 

at a depth of 15 feet produces variations of only 1 sercent in the 
volume of a gas. Because of the lowest water table elevation in the 
late fall, soil temperature would most likely be highest at this 
time. This would result in a maximum volume of the gas and a minimum 
soil moisture content through the displacement of ground water 

by a gas. However, a maximum soil moisture contcnt occurs during 

the late fall, the exact opposite of the expected result if 
temperature were affecting the volume of entrap»ed gases below the 
water table. 
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Confining pressure. The volume of entrapped gases below the 

water table would also vary with the confining pressure. The 

major factor affecting the confining pressure is the ycarly variation 
in the water table. A maximum elevation of the water table would 
minimize the volume or any dissolved gases and would consequently 
result in a maximum soil moisture content. However, a minimum soil 


6 


moisture content occurs in June, even at derths which were continually 


below the water table for several years. Another fact which tends 
to discount a significant volume change of dissolved gases because 
of changes in the confining pressure is the lack of an abrupt: change 
in the soil moisture content during irrigation. It was concluded 
that changes in the volume of entrapped gases caused by a change in 
the confining pressure were unable to account for the observed soil 
moisture changes below the water table. 


Exsolution and dissolution of dissolved gascs. Another possible 
phenomenon that must be considered to explain the changes in soil 
moisture below the water table with the passage of time is the 
exsolution and dissolution of dissolved gas from the water. A 
maximum temperature and a minimum confining pressure both occur in 
December. These two factors would favor the exsolution of any gases 
dissolved in the ground water. Thus, a maximum volume of gas 
exsolving from the ground water would have accumulated by December. 
4 maximum volume of gas should coincide with the minimum percent 
o£ soil moisture by volume. With the minimum soil moisture values 
occurring in June, cxsolution of dissolved gases from the ground 
water is not a satisfactory explanation for the increase in the 
percent soil moisture by volume below the water table. 


The evidence gathered to date docs not indicate that dissolution 

is the primary factor causing the cyclical variation in the percent 
soil moisture by volume below the water table with the passage of 
time. The water table has not declined to an elevation sufficiently 
low enough to allow the entire soil column penetrated by the soil 
moisture access tube to become partially drained. Partial drainage 
is necessary if dissolution of entrapped gases is to account for the 
observed changes in the percent soil water by volume below the 

water table with the passage of time. For cxample, the water table 
declined to only 11.0 feet below the soil surface at soil moisture 
access tube site 822 during the period from April to December 1969. 
Soil moisture determinations were made to depths of 17.3 feet below 
the soil surface. Other piezometers indicate that the water table 
rose approximately 2.5 fcet from December 1963 to April 1969. It is 
quite probable that a portion of soil moisture access tube 822 
remained below the water table for 1 or more years. It would not 

be likely that the partially drained soil and the soil below the 
lowest water table attained would have the same percent by volume of 
entrapped gases. Thus, the soil above and below the lowest water 
table attained should behave in a different manner. Vertical soil 
moisture profiles do not reveal any abrupt change in the behavior of 
the percent soil moisture by volume with the passage of time in the 


k 


vertical direction. It was concluded that the dissolution of entrappe 


gases could not account for the observed changes in the percent 
soil moisture by volume with the passage of time for depths that 
were continually below the water table. 
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5. Changes in bulk density. Yearly fluctuations in the water table 
could cause changes in the soil bulk density with a variation in the 
vertical stress. Since the actual compression of the soil was very 
small compared to the apparent soil water variation, changes in the 
bulk density were discounted as being capable of significantly 
affecting the soil moisture content below the water table. 






Equipment crror as a result of temperature sensitivity. The eifect 
of daily ambient air temperature on the scalers and rate mcter has 
also been considered as a possible phenomenon causing the variation 
of the apparent soil moisture content below the water table. This 
study is still progressing, but it docs appear as though the rate 
meter and possibly the two scalers are tempcrature sensitive. One 
problem with the data is that no record was made of the scaler used 
to collect the data. At present, the temperature sensitivity of 
the scalers and rate meter appears to be the most probable factor 
causing a time-dependent soil moisture change below the water table. 





The varigtions in percent volume soil water may seem insignificant for most 
uses of soil water information. However, for this particular study, soil 
water values were used to assess specific yield and ecvapotranspiration to 
compute water losses from the irrigated field. In comparing the corrected 
data with the uncorrected data, the differences are significant. Using the 
uncorrected data, the water loss determinations, including specific yield 
and transoiration, had a logarithmic mean of 10.4 percent. When using the 
corrected soil water data, the water loss values were 23 percent, which is 
reasonably close to the average water loss of 25 sercent for loam soils based 
on estimated porosity and field capacity data (Israclson and Hanson, 1967, 
pp. 411). The water yield values calculated from soil moisture and water 
level data are for a soil column extending from a depth of 10 inches-to the 
water table. 


The soil water losses in the upper 10 inches of the soil were determined by 
repeated gravimetric samples at 4 sites within the fieid. Linear regressions 
on the resulting plots of soil moisture against time indicated that the 
greatest rates of soil water loss from the upper 16 inches of the soil 
decrease from an average of 0.044 in./day in the upper field to 0.017 in./day 
in the lower field, which is finer textured than the upper field. These rates 
are appliable from May 30, two days after significant rainfall ceased until 
the estimated permanent wilting point was attained during the period from 
July 16 to August 2, 1969. Any soil moisture losses from the upper 10 inches 
ot the soil attributable to cvanotranspiration or gravity drainage was 
measured gravimetrically. ; 


An estimate of the long term effects such as ground-water outflow on’the 
water level declines for each shallow piezometer located within the 

irrigated field were based on the assumption that evapotranspiration was 
negligible immediately following the cutting of the alfalfa. First or second 
order regressions were fitted through a plot of the daily water level declines 
beginning immediately following the cutting of the alfalfa and terminating 
about two weeks after the cutting of the alfalfa. These equations were then 
used to calculate the net long term effects on daily water level changes at 
the site when the first crop was cut. This proccdure was repeated for the 
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second crop cutting in August. These two data points and the water level y 
data in the fall following the growing season were used to construct an iC 
equation defining the net long term effects on piezometer water level 

changes at the site throughout the period of study. The net result of long 
term effects such as ground-water inflow or outflow at the piezometer sites 
was measured during periods when evapotranspiration was assumed to be 
negligible. Ic is believed that the long term effects are caused primarily 

by the ground-water flow ates ae the long term effects primarily providing 
an estimate ofthe change in the water level at a sie. 





The daily losses due to evapotranspiration at each site adjacent to a soil 
moisture access tube were cvaluated by subtracting the daily water level 
decline caused by long-term effects from the observed daily ground-water 
decline and multiplying the remainder by the appropriate water loss factor. 
The water loss due to gravimetric soil moisture loss in the upper 10 inches 
was also included in the daily evapotranspiration loss. Preliminary results 
indicate daily evapotranspiration losses up to 0.38 in./day during very 

windy hot conditions with a mature crop. The daily cvapotranspiration appears 
to vary aerially in a manner dependent upon wind direction, andvelocity. This 
method of estimating daily evapotranspiration results in values which are 
reasonable for the field conditions. 


The evapotranspiration water losses for two hourly intcrvais between June 19 
and July 13, 1969, were calculated using the daily cvapotranspiration water 
losses. The two hourly changes in the piezometer water levels were corrected 
for changes in the barometric pressure. The barometric cfficiencies of the y 
piczometers were approximately 5 percent. Piezometer time lags and amp 1itud 
ratios (Hvorslev 1952) are satisfactory for measuring two hourly changes in 
the fluid potential of water adjacent to the piezometer screen. The: largest 
basic time lag for a piezometer equipped with a continuous water level 
recorder was 7.4 min., which would result in an equalization ratio of 90 
percent in 1/7 min. following a change in the fluid potential of the water 
adjacent to the well screen. Since all the piezometer-soil moisture sites 

did not receive their final irrigation at the same time, and particularly the 
lysimeter site which received its last irrigation on May 20 as compared to 
June Il for the surrounding field, adjustments were made to the daily gravi- 
metric and ground-water losses occurring below a depth of 10 inches. Final 
calculations of the two hourly evapotranspiration losses have not been 
completed at this time. The two hourly evapotranspiration water losses appear 
to have much larger relative errors associated with them because piezometer 
time lags and barometric cfficiencies must be taken into account. 


The coefficient of horizontal permeability was evaluated using Hworslev's 
(1952) method which uses piezometer screen geometry and basic time lag. 

The population frequency distributions for permeability are ail log normal. 
It is not possible to determine the ratio of the cocfficient of vertical 

to horizontal permeability using Hvorslev's technique. The mean permeability 
for the alluvium assuming Ky/K, = 1.0, 10, and 100 are 1.57 gal./d./é£t. 

eth de Belief Cp if Bt on nte ON ied 100 scala hel /£t.2 respectively. The mean permeability 
for ort Boston Ranch Unit adler the alluvium assuming K/K, = 1,0, 10, 
and 100 are 0.518 gal./d./ft.? 0.690 gal./d./ft.* and 0.863 gal./d./£t.2 é 
respectively. Even though it was assumed the Kyp/K, = 10, the effect of this' 
assumption on the resulting mean horizontal permeability is rather small 
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Figure 1, Variations in percent volume soil moisture at depths 
below the water table at selected sites. (Sites 142 and 
742 located in clay soils; 241 located in alluvial soils.) é 
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g eft: CRIS Work Unit No. SWC-012-fBo-2 Code No.: Idi.-Bo-103, 4 
Litle: Geochemistry of ground-water flow systems. 
Location: Northwest watershed Research Center, Boise 
Idaho. 
Personnel: iG... Stephenson and John Zuzel; 


Dr STE ot the. Sharp, University of Nevada. 





Date of Initiation: January 1968. Zxpected Duration: January 1971, 
Objectives: 


Il, sevelopmenit of chemically oriented methods for determining 
flow characteristics of ground-water systems of western 
ranzelands., 


2, 4 better understanding of the changes in water chemistry 
whicn occur within such systems, primarily as the result 
of reactions between water and geologic materials. 


* oF Need for Stucy: 


Dissolved constituents of natural waters result from combined reactions 
of water, soil, and rocit,and to a varying extent from man's wastes. 
In general, information is needed on the origin and occurrence of dis- 
solved constituents within a flow system and the full Significance of 
each portion of the system in mobilizing and transporting these con- 
stituents, 


There is a definite need for improvement in the knowledge of the chem- 
istry of ground water in flow systems associated with western range- 
lands, Cherical parameters offer possibilities as effective tools in 
unrzveling some of the more perplexing hydrologic problems heretofore 
not subject to unique solution by Panearional hydrologic methods. Prob- 
lems possibly amena>le to improved solution by geochemical methods 
include (1) age, amount, and residence time of ground weter;: (2) mixing 
of water frora te elise sources; (3) definition of arenas of recharge, 
lateral flow, and discharge; and (4) amount and location of recharge and 
discharge. 


Man's imposition of a variety of management conditions on rangelands 
. , not only upsets long-term hydrologic relationshins, but -lso disturbs 
®) \/ long-standing geochemical relationships, causing changes in water 

} quality which promise to accelerste over the years sahend. 
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Field Phase: It will be necessary to perform a complete hydrologic 


analysis of the flow system(s) under investigation, This will re- 
quire assembling all available hydrologic and geologic date, continued 
collection of ground-water data (water levels), possibly some expan- 
sion of present well networks, and installation of flumes end/or weirs 
at springs, /ater samples for analysis will be collected periodically 
at wells and springs, probably on 2 monthly basis. Saroples of pre- 
cipitation and of soil water for enelysis will be collected, using spe- 
ciel rain gazes and lysimeters. flow-measuring devices will be 
placed at the outlets of springs. Some chemical dye tracer tests and 
artificial infiltration studies may be performed, (Soil and bedrock 
Samples will be collected at various dleces for analysis and K xboratory 
studies, 


x 


Laborotory Fhase: Water samples will be analyzed for common chem- 


ical constituents and, in some instances, trace elements, Some 

samples will se analyzed for isotopic constituents, Some soil and 
rock samples will be analyzed chemicolly and mineralosically. Column 
or batch equilibrations will be mede of water and soils/rocks to observe 
chemical reactions which may be occurring in the syster(s). 
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Experimental Osta and Observations: 
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A sufficient number of water chemistry analyses, over - suiteble time 
span, permits statistical analyses of this data to be performed with a 
high level of significance. The relationships tested were positions in 
the flow systera as indicated by potenticl surface elevation and variation 
in ionic ratios throughout the basclt ground-water flow system. Anal- 
yses for the sedimentary aquifer system, and several local flow sys- 
tems, have not been evaluated so no eitempt will be made to comment 
on their significance, 


A correlation motrix was computed correlating water-level elevation 


with various ionic ratios which occur in the water samples. The toble 
of correlation coefficients and the voricbles are given in Table 1,. 


Table I points up several items of importance concerning the water 
chemistry of the basalt ground-water flow system. Silice (SiQ) cor- 
relates positively. This means thet progressing downstream in the 
flow system,the total percent of silice in the system decrenses, From 
previous year's analyses, silicn was always higher in the beseltic 

aters than in the waters of the sediment aquifers, -nd was considered} 
to be an incicctor to identify one system from another. However, the 
full significance of the presence of silice and its varintion wes not 
understood. 
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TABLE 1. --Correlation coefficients of variables versus 
water-level elevation, 
Varicble Water Level levation 


_ 
e 


Woter Devel Jlevetion—-— * 
Total 2gquivelents Per Million 


Total Anions + Total Cations/ Si@ 


e 


' 
QO 
oO oO 


OD . 


Cot Me/iNi + = if - 0,569 
Total Alsclinity /C1 + SO, — OaAcG 
Total Dissolved Solids - C, o13 

VU,uue 


SiC, 
ASRS oases aeRO a mo le Seattle Came AAA, NII Sein ee 


* No. of cnses = 75 for each Vrriczle 


All the other varicbles considered in Txble 1 show 2 notative 
correlation witn potential surface elevation, indicating - zenereal 
increcse in ionic concentration downstrenm in the flow system. This 
has been previously mentioned but can now be verified with the 1970 


Sampling anc cnalysis. 


Figure 1 shows an avernge hydrograph of the potential surface of four 
representative wells in the basalt flow system from 1967 through 1970, 
The avernge changes in the ionic retio of total anions + total cations 
(TDS) / SiO,, cetermined from 7: enses, is plotted -bove the hydro- 
graph in this figure. The effect cecharge has on the changes in this 
ionic ratio is evidence when reviewing Figure 1. The ing time and 
the variation in the ionic ratio are cependent upon the intensity of the 
recharge event and the time of yenr. The relation of the recharge 
event anc tre related change in TOO/SiO, is indicated by the arrows 
in Figure 1, In 1969 © large rechar-e 4évent was recorded, “A drop 
from 10.75 to $.C0 was recorded in the ionic ratio. ~urins the summer 
months the ionic concentrations incrense and 2 summer or fol] event, 
such ns in september of 1969, influence the ionic ratios to - lesser 
Gegree becouse the dilution effect is less significant. 


Figure 2 shows the relationship betweon position in the brsalt flow 
system, oS expressed by the elevntion of the potentinl surface -t ench 
of four wells, and the ratio of total go scolyed solids to silics:, 

Hach point shown in the figure ropresents an avernge of 1& potential 
surface clevations nnd 18 ratios of totcl dissolved solids to silice at 
each well site. 
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Potential surface hydrograph showing relation of recharge events to changes in ionic 
ratio of ground-water samples, 
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Figure 2, 


Regression showing relation of ratio of TDS/SiO, to position in eroundc-water 


flow system. 


r= —-0.996 
Y= -—37.89X + 4525.98 
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The high degree of negative correlation (-0.996) indicates that as ry 
one proceeas downstream in the basalt flow system, the ratio of 

total dissolved solids increases while silica content decreases re- 
sulting in a predictableincrease in the ratio of total dissolved solids 
to -silica, 





While total dissolved solids versus elevation and silica content versus 
elevation both produce statistically significant correlation coefficients, 
these relotionships are much less clear than when the ratio of total 
dissolved solids to silica is used in the correlation. 


Coraments, Interpretations, and Suture Plans: 


From the present correlation results, SiO, shows considerable in- 
stability. It varies from one sample site to another because of the 
presence or lack of secondary silica deposits which occur mostly in 
the upver portions of the flow system. The correlation matrix. of 
Table 1 shows that SiO, correlates positively. This means that 
progressing downstream in the flow system, the total percent of silica 
in the system decreases. This further explains Pigure 2, which shows 
the ratio of TDS/SiO, changing downstream in the flow system. The 
decrease in the ratio is the result of the decrease in 5iO0,. The 
reason for the decrease in SiO. is probably caused by diminishing é 
source of secondary silica and “increased dilution in that direction, : 
However, the lower values for SiC, in the lower noreas of the basaltic 
aquifer system are still considerably higher than in the flow system 

in the lacustrine sediments. 


When water chemistry information from other wells in the basaltic 
ground-water flow system is analyzed and plotted, as described above, 
and if the results fall within the limits of the regression of Figure 

; ao method has then been found to locate positions within a flow sys- 


tem, This would be very beneficial in water quality mapping. 


One well from which water chemistry information was obtained for a 
short record of time was used to test the above hypothesis. Tata 

from this well, when plotted as water level elevation against the ratio 

of TlDS/SiO,, fell well within the limits of the regression equation of 
Figure 2. “One test is not significant, but tests are now being run to 
try to develop this work further by using water chemistry data ‘already | 
gathdred for wells in the sediments and from local flow systems. | 
Additional samples from wells in the basalt grounc-water flow system 
will also be analyzed, The plan is to develop a family of curves like 
the one on Figure 2, for the flow systems in the lacustrine sediments, 
local flow systems, and the basalt. - 
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CRIS Work Unit No.: SWC-014-£Bo-3 Code No. Ida-Bo-107.1 
Title: Sediment yield from rangeland watersheds. 
Location: Northwest Watershed Research Center, Boise, Idaho. 
Personnel: ; C. W. Johnson, D. L. Schreiber, G. R. Stephenson, 





W. R. Hamon, G. A. Schumaker 





Date of Initiation: September 1, 1969 
Expected Duration: Originally zlanned - December 1974 


Present recommendation - December 1974 


Objectives: 


1. To determine the relationships between sediment yield and variables 
describing hydraulic and hydrologic factors, and site and watershed 
- Characteristics which influence sediment yield. 


N 


To formulate a sediment yield prediction procedure for rangeland water- 
sheds in the Northwest. 


Need for Study: 


Information on sediment yield is almost entirely lacking for millions of 
acres of predominantly sagebrush rangeland under Government land manage- 
ment and private ownership in the Northwest United States. There is a 
growing concern for soil losses from intensively-grazed rangelands, 
sediment damage to reservoirs, and erosion of stream channels. 


Most rangeland watersheds in the intermountain Northwest have large areas 
of relatively steep hillslope topography and these areas need to be 
delineated for treatment to reduce erosion. Also, sediment yield informa- 
tion is needed for evaluating the benefits of watershed management and land 
treatment programs of the Bureau of Land Management and Soil Conservation 
Service. 


Range sites found in the Reynolds Creek Experimental Watershed repre- 
sent a large percentage of the rangeland in the Northwest and studies 
of sediment yield are essential for the development of sound management 
practices and in planning appropriate multiple use of these lands. Good 
land management decisions require information on how vegetative changes, 
fencing and land use alters the sediment yield potential of rangeland 
watersheds. The sources and quantities of these sediments need to be 
determined and measured so that research data can be used to predict 
sediment yield for ungaged areas in terms of available information on 
soils, climate, physiography and use. Research is also needed to adapt 
the universal erosion equation to rangelands. 
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Design of Experiment and Procedures to be Followed: 


Suspended and bedload sediment yield from plots, channels and watersheds 
will be measured by use of pumping sediment samplers, splitting devices, 
catchments and hand sampling. Scour and fill within designated reaches 
of channels will be measured by photogrammetric techniques utilizing 
captive balloons or helicopters for obtaining photographs. 


Runoff hydrographs from plots and watersheds will be measured through 
flumes, weirs, or in tanks to obtain data for computation of hydraulic 
parameters, The effects of precipitation impact and splash on soil 
erosion and overland flow will be studied using 14 ft. by 72 ft. plots 
and simulated rainfall. Critical velocity and depth of flow which 
influence erosion for particular soil and slope conditions will also be 
determined from plots. Plots, microwatersheds and watersheds will be 
located on various soil types in different precipitation zones. A wide 
range of slope length, slope area, aspect and relief ratio will be 
represented. 


Erodibility of watershed soils will be determined from rainfall simulator 
plots to determine the effect of parent material, and soil texture, structure, 
aggregation, organic matter, and pH on soil movement. Large-scale soil 
movement will be traced with colored or tagged particles. 


Data on vegetative, litter, and rock cover will be obtained from a related 
research outline where the procedure will be described in detail. ca 
Physiographic and geomorphic factors will be determined throughout the 
Reynolds Creek Watershed. Detailed determinations of aspect and percent 
Slope will be made by use of a systematic grid and a data reduction pro- 

gram. Hypsometric (area-altitude) information obtained for the study 

sites will be used to designate othcr areas of corresponding sediment 
production potential. Geomorphic factors will be related to sediment 

yield measurements. 


Rainfall intensity and duration data are available from a network of 
precipitation gages and snow data are available from snow courses, snow 
pillow sites and other snow-measuring sites. Research Outlines Nos. 100.1 
and 102.1 describe the instrumentation and procedure. 


Evaluation of management and land use factors will be accomplished by 
selection of watersheds in areas with the same management and then 
fencing off one- to three-acre areas to exclude livestock and establish 
maximum vegetative cover. 


Experimental Data and Observations: 
1. Overland Flow Hydraulic Parameter Study 
To study surface runoff and associated hydraulic parameters (velocity, 


depth, roughness, type of flow) which influence sediment yield, a laboratory 
study was conducted by D. L. Schreiber in coopération with the ARS Erosion 


SN ooo ar tae aoe cae ade eee A oan ie ict ee eee deal 


ik as eae Ee is ee 


; 4 ‘ a a oR TTY STON SST Ca FOC RET LO LO 
PEP Sa Ma ROL ORS Pee eT a AOR MATE A Per aL A 


13-3 


} 
. neat Laboratory at the Palouse Conservation Field Station, Pullman, Washington. 
(CRIS Work Unit No.: SWC-016-fPul-1; Research Outline No.: Wash. Pn-66-15, 
terminated June, 1970.) 


A mathematical model was developed to simulate and predict overland flow 
hydrographs. The kinematic form of the nonlinear partial differential 
equations of unsteady, Spatially varied, open-channel flow were solved 
simultaneously by numerical integration. The mathematical model was 
verified by laboratory data obtained in the Erosion Laboratory Rain Tower, 


Since surface roughness is practically impossible to measure quantitatively, 
parameter optimization was used in conjunction with the mathematical model 
to provide a set of representative, synthetic, laminar-flow resistance 
parameter values. Three different univariate Optimization schemes were 
devised for this study. The first scheme produced “wholly optimized: 
hydrographs for which a Single value of the resistance parameter can be 
used to describe the resistance to flow both during and after the occurrence 
of rainfall. The second scheme produced ‘two-piece optimized’ hydrographs 
that were considered to be composed of two sections - during rainfall and 
after rainfall. An optimal value o£ the resistance parameter was determined 
for the rainfall section before the recession-section resistance parameter 
was evaluated, The third scheme produced ‘sectionally optimized’: hydro- 
graphs that were divided into sections or time intervals of either 10 or 

oN 12 seconds. The optimal resistance parameter value for each section was 

) Kot determined in succeeding order. 


Results of the three optimization schemes were compared to determine the 
sensitivity needed to describe adequately the resistance to Laminar flow. 
Results of the third scheme, sectional optimization, were also used to 
provide synthetic data for use in an attempt to develop predictive equations 
for the laminar flow resistance parameter, 


Detailed descriptions of the laboratory facilities and procedures, the 
experimental and synthesized data, and the development of the mathematical 
model and optimization schemes arc contained in D. L. Schreiber's Ph.D. 
thesis, “Overland Flow Simulation by a Nonlinear Distributed Parameter 
Model,*’ Washington State University, June 1970. Representative data and 
results are indicated below in Table 1 and in Figures 1-4. (All figures 
follow page 13-11). 


2. Sediment Yield Studies 
a. Sediment sampling equipment and facilities 


Sediment facilities and a description of instrumentation Operating in 1970 
are listed in Table 2. 
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TABLE 1.--Optimal values of laminar-flow resistance parameter K 
for two-riece ontimized hydrographs. 
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Test Slope, in Intensity, in Duration, inj} K During K After 

Number } feet per feet | inches per hour minutes } Rainfall | Rainfall 
(1) (2) (3) (4) aso: (5) (S) 
10 0.0075 2392 10 36 LJ 
bz 0.0075 otbk 4.0 29 16 
19 0.0465 2 eye 1.0 39 15 
20 0.0465 2. 84 4.0 thr a 
24, ! 0.0465 1. 06 4,0 20 13 
27 0. 0806 34 13 
22 0.0806 40 16 
33 25 14 
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TABLE 2.--Sedimentation facilities and instrumentation 


Northwest Watershed Research Center, 1970. 
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NO. 


036068 


046017 


046084 


048077 


097098 


116083 


138012 


166076 


015094 


022007 


LOCATION DESIGNATION 


NAME 


Outlet Weir 


Salmon Cr. 
Drop-Box Weir 


Macks Cr. 
Drop-Box Weir 


Summit Drop- 
Box Weir 


Nancy Gulch 
Runoff Plot 


Tollgate 
Drop-Box Weir 


Upper Sheep Cr. 
Drop-Box Weir 


Reynolds Mtn. 
V-Notch Weir 


Ertete”-Rabbit Cr 
Drop-~Box Weir 


Rabbit Cr. 
Flume 


SEDIMENTATION 

INSTRUMENTATION 
Hand Sampling 
Hand Sampling 
P65. D/P Urb ne 
Sampler 


Gravity Sampler 
Bedload Catchment 


Tanks and 
Splitter 


P.S. S/7 Pumping 
Sampler 


Chickasha Pump Samp. 
Bedload Catchment 


Chickasha Pump Samp. 
Bedioad Catchment 


Chickasha Pumping 


Sampler 





DESCRIPTION 


During Flood Events 


During Flood Events 


Float Actuated 


Float Actuated 
Following Events 


Following Events 


Float Actuated 


Float Actuated 
Following Events 


Float Actuated 
Following Events 


Float Actuated 


Float Actuated 


Hand sampling of suspended sediments at the large weirs during flood 
events continued through 1970; however, limited personnel and travel 
difficulties prevented collection of much valuable data at sites where 


pumping samplers were not installed or not in operation. 
the P.S. 67 pumping samplers were common during the past year due to 

inoperative float switches, pumps, solenoids, sample trap mechanisms and 
Plumbing leaks, battery failures and freezing also 
The Interagency Sedimentation 


electronic controls. 
caused frequent maintenance and poor data. 
Laboratory has made numerous changes and improvements on the samplers 

during the past year and two of the four samplers have been altered by 


Failures of 


a E> installing new trap mechanisms, air solenoids, tank float switches and 
NY ( % 
} 


ae electronic controls. 


f instrumentation and loss of much important data. 


» eee Mee re a ETE TA 7 RR I ET IT EI I ET Ee TIT er pee OY RR RA nt TN 


These changes have caused delays in completing 


eg Sr TY I - 


: : ETO NS Ys Hn ON PR BE TR OEE Ter STO POL AI PATIO LY RRR OSH OER EDS STs UY PORT OT ERTIES OS GT 





The Chickasha pumping samplers cannot be effectively used at sites where 
quantities of sand may be pumped through the system because the pump 
impellers soon fail. H-wever, these samplers have proven very reliable 
at sites where only fine grained sediments are in suspension. The 
Chickasha sampler is an excellent device for instrumentation of small 
watersheds and fine grained sediments. Also, screening of the intakes 
to prevent plugging by flood debris is very important since no back- 
flushing system is available with the sampler. 


The purchase, installation, and maintenance costs for automatic sediment 
sampling equipment and facilities has been extremely high because of 

(1) rapidly increasing initial cost of equipment, (2) continuing redesign, 
(3) need for spare parts when repairs and improvements are being made, 

(4) need for heated facilities and (5) frequent servicing, adjustment, and 
repair? 


Two gravity samplers for suspended sediments were designed and fabricated 
at the University of Idaho under the direction of Asst. Prof., Myron Molnau, 
and one of the units was installed at the Summit Drop-Box Weir. However, 
the sampler failed frequently during field testing and improvements are in 
progress. At sites where free overfall exists, the gravity sampler operates 
without use of a pump and a complicated switching mechanism. 


Instrumentation of a 14.ft. by 72-ft. runoff and sediment plot was completed 
ae in 1970 on the Nancy Gulch watershed. Runoff and sediment are measured in 
tanks following each runoff event and the system capacity accommodates one 
inch of runoff. 





Catchments for measurement of bed léad at weir sites have been constructed 
on three watersheds as shown in Table 2. The accumulations of sediment 
are measured following runoff events and then cleaned for the next event. 


b. Sediment yield from January runoff 


January 1970 precipitation amounts exceeded previous records for any month 
at several weather stations in Idaho. Generally, the precipitation occurred 
on 20 consecutive days from January 9 through 28. Similar mid-winter 
storms have occurred on the Reynolds Creek Experimental Watershed in six 

of eight years since records began and caused some degree of flooding. 
Comparative data on precipitation and peak flow are shown in Table 3 for 
the major mid-winter storm events. Usually, the rain is accompanied by 
warm temperatures and wind which causes rapid melting of the snow on frozen 
ground. However, the 1970 flood peak occurred on January 27th after 19 © 
consecutive days of warm temperatures and alternating rain and snow which 
thoroughly thawed and saturated the soil. 
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TABLE 3.--Precipitation and neal: flow for mid-winter storm events 
Reynolds Creak Experimental Watershed. 


Year Storm Peale! Precipitation at Selected Stat ions2/ 


Period ~« F Low ey RE 4/ 5/ 
R.G. 076459 Beis. LODGOs. R.G,. 176407. 





(OSES) (INCHES) (INCHES) (INCHES) 
1963| Jan. 28-Feb. 3 2330 ZUG 2.60 SPSS 
1964} December 3800 eos, En ye, 10.536 
1965} January 1100 Daye bed Apes 
1967] January 266 Leni Saf 2 9.49 
1968; February : a | 1.00 3.92 3.00 
1969} January 900 Ja, 0o 3G. 14 2203 
1970} January 730 eee 9.14 ae 


i/ Outlet Weir, drainage area of 90 sq. miles 
2/ Unshielded gages 

3/ Raingage elevation 3915 ft. 

4/ Raingage clevation 5410 ft. 


5/ Raingage elevation 6800 ft. 


Sediment yield per unit area was greatest at the lower elevation watersheds 
as a result of January storm runoff because che deeper snow at higher 
elevations absorbed most of the rain. Sediment data from three representa- 
tive stations at various elevations are summarized in Table 4. Storm 
characteristics, runoff, and sediment yield for the January 1970 event were 
very similar to the event of January 1969, as shown in Figure 5. On the 
basis of nearly ten years of record, mid-winter runoff events of near-flood 
magnitude occur about three years out of four. The usual combination of 
warm temperatures, wind, rain, snowmelt, and frozen soil have caused the 
major floods of record with associated high sediment losses. The January 
1970 event appears somewhat typical of mid-winter events common to inland 
areas of Northwestern U.S. aithough of greatcr-than-normal magnitude. 
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TABLE 4.--Summary of runoff and suspended sediment from selected water- ® 


sheds, 1970. 











Watershed Drainage Runoff Suspended Sediment Yield 
and No. Area January March-June January March-June 
(Acres) (Inches) (Inches) (Tons/acre) (Tons/acre) 


Reynolds Mt. 


East, 166076 1G0 0.38 en Tas 0 0.216 
Reynolds Cr. 

LOL Usivs 

116063 DSe353 ae 2 Pi Osab Ly, 0.409 
Macks Cr. at 
Weir, 046084 7,846 1.01 Dad5 0.465 0.054 


c. Sediment yield from spring snowmelt runoff 


The water content of the snowpack at higher elevations of the Reynolds Creek 
Experimental Watershed was greater than normal at the beginning of the spring 
snowmelt season; however, January storms caused melting of much of the snow- 
pack below 5,500 feet elevation. igure 6 shows the graphs of temperature, @ 
runofi, and suspended sediment for the snowmelt seasons of 1969 and 1970. 
Typically, the major portion of the seasonal sediment yield is associated 

with peak flow during periods of highest temperature. Table 4 above contains 
a summary of data on runoff and sediment yield for 1970 at selected stations. 
Sediment concentrations are extremely Low during periods of low flow and do 
not contribute significantly to the total annual sediment yield amounts. 


d. Runoff and erosion plots 


A 14-ft. by 72-ft. piot in the Nancy Gulch ‘Jatershed was instrumented for 
measurement of runoff and sediment from natural events. An S-inch high metal 
border, collector, tanks, splitter, and FW-1 recorder have been installed 

for measurement of runoff and sediment from this typical range site. Also, 
about 20 other plots have been staked at five other sites which represent a 
wide range ox elevation, slope, cover, and soil conditions. The completed 
plot instrumentation is designed to catch one inch of runoff from the plot 

ard to be pumped and cleaned following each runoff event. 


e. Cooperative sedimentation researcn 


Cooperative sedimentation and hydrology research between the Northwest 
Hydrology ,Research Center and the Water Resources Research Institute, 
University of Idaho, Moscow, Idaho, ted to the development and fabrication 
of two gravity-flow suspended sediment samplers for use on a watershed near 
Moscow, Idaho (Thompson Watershed), and on the Reynolds Creek Watershed. 
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However, no runoff events have occurred at these stations since installa- 
tion of the samplers. The 1970 annual report to the Agricultural Research 
Service by Myron Molnau and D. J. Davis contains data on precipitation, 
runoff, sediment yield, temperature, radiation, humidity, wind, and evapora- 
tion for the Thompson Watershed or nearby stations. During 1970 precipita- 
tion was about 30 inches, runoff about 5 inches and suspended sediment 
approximately 0.044 tons per acre. The report also shows a schematic 

of the gravity sediment sampler and explains the sampler mechanism and 
controls. 


Comments, Intcrpretations, and Future Plans: 


1. Overland Flow Hydraulic Parameter Study 


An overland flow plane is a very important component of a total watershed 
system. Furthermore, much of the sediment yield from rangeland watersheds 

may originate from surfaces that could be classified as overland flow planes. 
Therefore, one emphasis of this sedimentation study has been and will continue 
to be an analysis of overland flow. 


The synthetic laminar-flow resistance savameter data that was obtained from 
the sectional optimization was used in conjunction with a dimensional 

analysis and a least-squares analysis to yicld predictive resistance parameter 
relationships. As a result, the following power equations relate the re" 
Sistance parameter to the precipitation number P (precipitation rate divided 
by the product of unit discharge and downstream depth) during rainfall and 

to the Weber number W (the product of unit discharge, downstream velocity, 

and fluid mass density divided by surface tension) after rainfall: 


R= 6.95 po- 1431 


and 


1666 


0. 
K = 15.64 W 


The predictive power equations were tested by including them in the numerical 
solution of the mathematical model. Resulting predicted hydrographs com- : 
pared closely to the observed hydrographs, as shown in Figures. The 
exponents of the power equations are small in magnitude (reflecting a small 
slope for the regression lines), thus Supporting a hypothesis that two care- 
fully selected constant values of the resistance parameter are adequate to 
reproduce the laboratory hydrographs. This hypothesis is supported by the 
results shown in Figures 2 and 4. The value of the resistance parameter | 
that is chosen for the rainfall period must be larger in magnitude than the 
value chosen for the recession period, as indicated in Table 1 and Figures 

2 and 4. The fact that two constant valucs may be used satisfactorily 
suggests that the hydrographs are fairly insensitive to changes in the 
precipitation number and the Weber number. 


In another study a stochastic streamflow model was used to sequentially 
generate synthetic records for the Cascade River et Marblemount, Washington. 
Statistics of generated flow sequences compared closely to the respective 
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of the 100-year maximum and minimum flows also compared reasonably to thos 
obtained by sequential generation. 


Statistics of the 36-year historic floy sequence. Gumbel method tothe 


Future plans for further overland floy analysis under this research out Line 
include a field evaluation of the mathematical model, Field plots, 14 fe; 
by 72 ft., have been selected at various locations on the Reynolds Creek 
Experimental Watershed. Borders and water and sediment measuring devices 
will be installed during 1971 and 1972. Artificial rainfall will be applied 
to the plots with a portable Sprinkler system. Runoff hydrographs will 

be measured in conjunction with sediment yields in an attempt to rclate the 
overland flow parameters to sediment yields. The runoff hydrographs will 
also be used to test further the mathematical model of overland flow. 


Other future plans include developing a new research outline during 1971 
for further study of surface runozs., Further analysis of the overland flow 
model will be conducted using previously published data by other authors, — 
This analysis will provide an unbiased check on Schreiber's model, since 

it was developed from a limited scot of laboratory data. This analysis is 

a necessary prerequisite to the plot-runoff field testing program. The 
overland flow model will then be expanded to include a channel network and 
several overland flow planes. The nay Surface runoff mathematical model — 
will be tested by runoff data obtained from the laboratory (Palouse Rain 
Tower and other published data) and the field (Reynolds Creek Experimental 
Watershed and other published data). A comparison of results obtained by 
the proposed surface runoff model will also be made to results obtained 

by the instantaneous unit hydrograph approach. & 





Publications completed and in preparation since last year's report are the 
following: 


Schreiber, “D. L.-1970, “Overland flow simulation by a nonlinear distributed 
parameter model. Ph.D. thesis. College of Engineering, Washington State 
University, Pullman, Washington, June. 


Schreiber, D. L., and Bender, D. L. Sequential generation of annual stream- 
flow. SWC approval granted September 1970 for publication in the Trans- 
actions of the ASAE. Submitted to ASAE November 1970, 


Schreiber, D. L.; and Bender. D. 1. Obtaining overland flow resistance by 
optimization. SWC approval requestcd January 1971 for publication in 
the Journal of the Hydraulics Division, ASCE. 


2. Sediment yield studies 


There is a definite need for a greater number of sediment measuring stations 

in Idaho, Eastern Oregon and Eastern Washington if erosion, stream sediments 
and reservoir sedimentation are to be measured, understood and controlled. 

Only a few of the numerous stream gaging stations are instrumented for sedi- 
ment sampling although serious erosion and sedimentation problems are evident 
in many places. Therefore, a networ!: of sediment measuring stations are i 
needed to supplement long term data from stream-gaging stations. 


The accurate measurement of total sediment from plots, watersheds, and 
large drainage basins requires the use of the best available sampling devices 
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and techniques. Therefore, efforts must be continued to develop and im- 
prove automatic suspended sediment samplers for field use. 


Experience at the Northwest Watershed Research Center has proven (1) that 
procurement, installation and maintenance cf automatic sediment samplers 
are very expensive at remote sites when operated to monitor winter floods 
during freezing weather, (2) that less expensive Chickasha and Interagency 
P.S. 67 pumping samplers are useful only at certain sites and have limited 
reliability, and (3) that improved sampling devices and techniques are 
needed under the difficult field conditions encountered. 


Above average winter precipitation and snow accumulation in 1970 caused 
runoff and sediment yield similar to that of several other years of record 
and is, probably, representative of conditions to be expected in about three 
out of four years. Soil loss from major Reynolds Creek Watersheds approaches 
one inch in 100 years on the basis of limited data. 


With continued support and cooperation of the Bureau of Land Management, 
measurement and analysis of sedimentation on plots, from source areas and 
from streams will continue, improve and increase. Also, all available 
data is being processed and analyzed for publication. 


The gravity suspended sediment samplers at Moscow, Idaho (Thompson Watershed) 
and at the Reynolds Creek Watershed will be thoroughly field tested during 
natural runoff events to determine the reliability of the equipment. The 
usefulness of these samplers for runofi-sediment plots and small watershed 
studies will be investigated. A raper and slides describing the gravity 
samplertere presented at the 25th Annual Meeting of the Pacific Northwest 
Region of the American Socicty of Agricultural Engineers, October 7-9, i970, 
as noted: 


Johnson, C. W., and Molnau, Myron. 


A gravity sediment sampler for drop-box weirs. 25th Annual Meeting 
of Am. Soc. of Agr. Engineers, Oct. 7-9, 1970, Bozeman, Montana. 
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Fipure 3, Observed and predicted hydrographs (based on power equations for K) for 
Test 27, 
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Figure 4, Observed and predicted hydrographs (based on K = 24 during rainfall and 
K = 14 after rainfall) for test 27. 


ca feed WOES al ead 


t 


| 
| 
a i 
{ 


Pe 
t 
¢ 
+ 


ate pt oe 


+ 4 

4 

ete , 
+ 





HE eSB baae: 


' : 
F 
H , ' ' 
| ; 
{ : 
t 
ee z ae reer a Oe ear ae . 
‘ { a7 i 
i t . ! 
: i i , Her 
1 male aes jee a e—6-- _ uss - = 
4 A ‘ ! . 
; h ! i 
i he os ) 
{ ' , i 4 
| j i 
ae 5 : u 
i i i 
nt : \ ! 
+-- - = = ~ os = 4 —— 
i H oi ; t : 
Z ae 5 ws +—4 : eee 
t Doty ’ f ; 
\ ; A. I 
' | 1 } i ' 
Ef Ta, . Tow seks | 
ek oe 4 ont |B 4 os eae ee 
' ' ' j ‘ { . ' 
{ aaa 
' i ' 
, i oD 
+ : ey oe 
t ; 
: i Ye ' x A Da 
: H } i 1 
? } i] ' i 
ae te 
} : , = 
f i ; [4 
: 4 ie t 
‘ { ‘ , 1 | { 
’ ! ‘ ' 





Day of Month 


Ha 

2S I eee. 2 eh eees 
1 t u ‘ A ! ; Hi 
: i 1 





at anata 





os (aid SEY Ok 
Be eee, 
HS 


i 

1 
; 

i | 


(seysut) 


*dyderg poqernunssy (Aeq/suo],) pl@IA Jusautpag *sng AT eq (*s°y°o) a8aeyosiq ATyeq e8er9ay 


Precipitation, runoff, and sediment yield (suspended) for Macks Creek 
1970, 


Watershed, January, 


Figure 5. 


Sediment (Tons/day) Maximum daily temperature (°F) 


Daily discharge (c.f.s.) 


Figure 6. 
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Temperature, runoff and suspended sediment yield, Reynolds Creek 
at Tollgate, 1969 and 1970. (Temperature from Reynolds Weather 
Station) 











SUMUATION OF Ii:PORTANT FINDINGS 


SiJC-O11-fB0-1 Influence of climatic, DLiologic, and 
SS 


physical factors on rangelend water~ 
Shed hydrology 


Data collected during a rain-on-snow event demonstrated that 
the snowpack can store a considerable anount of Water -0oOf i the 
1.52 inches of precipitation that fell as rain (measured by 

an unshielded gage), only 1.44 inches of water cane out the 
bottom of the snovpack. Therefore, at least 0,38 inch of rain- 
water Was stored within the ll-inch snowpack, with an initial 
Water equivalent of 1.28 inches, (ida-D5o0-100.1) 

Computed precipitation, using data fron shielded and unshielded 
gages, for the water year 1669-70 in the Reynolds Creek Experi- 
mentval Watershed, exceeded the precipitation caught by un- 
shielded gages by about 1 inch where the total precipitation 
caught by the unshielded gace was 10 inches, but exceeded the 
catch by nearly 14 inches at higher clevations where the total : 
catch was about 31 inches. The unshielded gage catch was 

only about 69 percent o£ the computed precipitation at the 
higher elevations where the greater perecntage of precipitation 
occurred as snow and where wind speeds are greatest, 
(Ida-S0-102.1) 


f&. developed mathematical model for computin actual precivita~ 
tion £xom shielded and unshielded gage data, independent ox 
wind and type of precipitation, utilizes wind functions for 
the ratios of catch by shielded and unshielded gages to actual 
precipitation, Data fox the direct computation of actual 
precipitation are not available, except in the Reynolds Creek 
watershed, since precipitation is niversally observed by use 
of only unshielded or shielded gages. Coefficients have been 
derived by use of the mathematical model, and data on catch by 
shielded and unshielded gages, temperature, and wind to convert 
catches by either unshielded or shicided gages to estimates 

of actual precipitation as a function o* tind anc temoerature, 
The ratic for unshielded gage catch to actual precipitation 
for temperatures of 23° to 32°F. is 0.71 for a wind of 10 
mepeh. and is 0,36 for a tind of 30 mepeh. The coevfficients 
derived makes it possible to convert precipitation méasured 

by ordinary gages to estimates of actual precipitation by use 
of wind and temperature data, (ida-Do0-102.6) 


Data collected during a long period oz high continuscus snormelt 
Pay 


further demonstrates the differential response, of the 12 .it,. 
butyl pillow and the universal gage to decreases in water 


equivalent of the snowpack,. Under high melt conditions the 
water equivalent measurement as made -y the universal gage 
Clearly exhibits distinct diurnal changes as compared to only 


continuous linear decreases for the 12 £t, ~utyl pillow 
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measurement. This difference in response appears to be associ- 
ated with dixferential settling of the snowpack on the flexible 
surface of the 12 ft. butyl pillow and the temporery mounding 
of water on its surface. (Ida-50-102.7) 





Studies o£ channel roughness in laboratory flumes and natural 
channels have shovm the complex interaction of: (Ll) shape; 
spacing, and size of bed roughness elements; (2) geometry and 
slope of the channel; and (3) £low velocity and turbulence.s 
Parameters to describe various bed element spacing, size 
distribution, intensity and effective height have been useful 
in comparing randon and regular flume roughness patterns to show 
that maximum flow resistance for a specified intensity occurs 
with random spacing of elements. <.1S0, narrow channels have 
greater flow resistance than vide channels with equal bed 
roughness. (ida-Bo-104.1) 


Removal of a dense sagebrush growth by grubpbing and by spray- 
ing with 2-4-5-T resulted in a threefold increase in herbage 
yield compared with an untreated control area after two growing 
seasons with no grazinge. Measured soil moisture depletion was 
greater and snov accumulation less where brush had been removed.e 
Studies at other sites with different brush cover, annual pre- 
cipitation, soils and grazing use are in progress to fully 
verify the favorable forage yield increases Imdlcaced on the . 
first study site and to verify any differences in soil moisture 
depletion. (Ida-S50~-105,.%) = 





A mathematical model of a steady-state, two-dimensional flow 
system resulting from snormelt infiltration on a viatexrshed 

slope has been partially yverizied by 1970 field data obtained 
from the north slope of Upper Sheep Creek ‘Jatershed (WJ-17, a sub-= 
basin of the Reynolds Creek Experimental “Yatershed ‘ele 

Analysis of the field data indicate that steady-state con= 
ditions were not fully achieved curing the snovmelt period of 
1970, The mathematical model does, however, indicate what 

types of additional observations should be obtained. Limited 
field data verify to a degree, the indications from mathe- 
matical solutions that the saturated hydraulic conductivity ot 
the soils (1) decreases with depth below the ground surface, 
and (2) decreases upslope from the stream channel Sor about one- 
third the distance, then increases sharply at this, point, and = 
decreases again upslope until near the top of the slope where 

it once more increases, (Ida-50-105.5) 

Mathematical models o£ transient, partially saturated, one= 
dimensional vertical and three-dimensional axisymnetric, flow 
through soils have been developed and partially verified by 
laboratory data obtained from a soil from the :eynolds Creek 
Experimental “Jatershed. A modified Burdine Theory provides 

a functional relationship for obtaining, £rom saturation= 
capillary pressure data, the change of relative hydraulic be 
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conductivity with capillary pressure. Reasonaply good agree- 
ment is attained between observed conductivity-pressure data 
and the results predicted by the modified Burdine Theorye 
Solution results from the two mathematical models indicate that 
boundary effects on circular infiltrometers Significantly alter 
the flow pattern, even reducing the saturation at the surface 
centerline appreciably over that which would exist for the 

Same application rate over an infinite area, (Ida-Bo-105.6) 


A rainfall simulator-gamma probe-infiltrometer with a combina- 
tion of capillary needle sizes, air and water pressure, is 
capable of duplicating the median drop sizes of natural rainfall 
up to intensities of 4 inches/hour. The control of drop sizes 
and intensities also make it vossible to produce a desired 
rainfall energy even though terminal velocities are not reached, 
(Ida-Bo-105.6) 


Laboratory tests of a new, commercially available two-probe 
(gamma) density gauge utilizing a tracking, differential pulse 
height discriminator for temperatures ranging from 37° to 80°F, 
indicates that density readings (both means and variances). are 
not significantly different. The probable error in the measure- 
ment of soil moisture on a large laboratory sample over a 
Saturation range of 55 to 85 percent was 0.63 percent. 

( Ida-B0-105 e 6 ) 


Estimates of evapotranspiration by the Bowen ratio method and 
by energy balance-combination equations are in agreement for 
irrigated alfalfa, but the Bovien ratio estimates are consider-~ 
ably smaller for semiarid sagebrush rangelands. The apparent 
difficulty in using the Sovwen ratio method in semiarid environ-— 
ments lies in the assumption of equality in the transfer co-~ 
efficient of heat and water vapor and in obtaining reliable 
measurements of very small vapor pressure gradients. The 
popular energy balance-combination equations for computing 
evapotranspiration that require humidity data in addition to 
Surface temperature and radiation data that are required in 

the basic energy balance equation are of linited use since 
their results are redundant. ‘their only purpose is in the 
formulation of an equation for estimating potential evaporation, 
where for a saturated surface, the inclusion of air humidity 
eliminates the need fox surface temperature data. (Ida-5o0-106.1) 


Water balance computation for watersheds in southwestern Idaho . 
during the water year 1969-70 indicate that nearly all precipi-~: 
tation was lost to evapotranspiration in areas with less than | 
15 inches of precipitation while 67 percent tiaas lost in the 20-: 
to 40-inch precipitation zone. The water yield was 10 inches 
from the 20- to 40-inch precipitation zone, (Ida-5o0-106.1) 


For a semiarid low sagebrush site, the surxace albedo was 
observed to be relatively constant at 13 percent in early Junee 
The net radiation was estimated with reasonable accuracy as a 
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linear function of incoming solar radiation. A close 
correlation was also found between the soil heat flux and net 
radiation. A regression of evapotransniration on net radiatio§ 
yielded a simple correlation of 0.77. (Ida-So0-106.1) 





Siv¥C-012-fBo=-2 Ground water in relation to management of 
rangeland watershed in the Northwest, 


A systematic change was noted in percent volume soil water in 
the saturated zone at constant depths of measurement under an 
irrigated field, The cause may be due to instrument error as 

a result of temperature sensitivity. Soil water information, 
including corrected soil water values for the saturated zone, 
and changes in the piezometric water levels, were used to assess 
specific yield and evapotranspiration to compute total water 
loss. The loss values computed to date average 23 percent, 
which is reasonable for loam soils such as these as mentioned 
in the literature, (Ida-Zo0-103.3) 


From a correlation analysis of water level elevation versus 
ionic ratio from ground-water samples in a basalt aquifer flow 
System, it was found that silica (SiO0o0) decreases dovmstream 

in the flow system, A negative correlation (-0.996) for total 
dissolved solids to silica (TDS/Si05) versus water level 
elevation, shows that TDS increases dovmstream in the flow 
System. Analyses to date indicate that these characteristics 
are predictable and can be used with confidence for water ) 
quality mapping purposes in the area of Study. (Ida-Bo-103.4) 





SiJC-O11-£50-3 iifect of runoff, precipitation, climate, 
soil, vegetation, land use, and land form 
con sediment yield, 


A mathematical model developed to simulate and predict over~ 
land flow hydrographs was verified by laboratory data. Param~ 
eter optimization provided a set of representative, synthetic 
resistance parameter values. The synthetic data and a dimen- 
Sional analysis yielded predictive resistance parameter rela-~ 
tionships. Resulting predicted hydrographs compared closely 
to the observed hydrographs, Satisfactory results are also 
obtained with a constant resistance parameter value during 
rainfall and a constant of lesser Magnitude during recession. 
AS indicated by the results of this study, sediment yield from 
overland flow planes is probably initiated during rainfall 
periods, since turbulence effects are greater. (Ida-bo-107.1) 


Two years of experience with three interagency P.S. 67 pumping 
samplers have resulted in numerous failures involving almost 
every mechanical and electronic part of the sampler systemse 
Therefore, continued development and improvement of automatic » 
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suspended sediment sampling equipment is very important. Also, 
a network of suspended sediment sampling stations is urgently 
needed in Idaho, Nevada and Eastern Oregon to compliment long 
term runoff records throughout the area, (Ida-30-107.1) 


Sediment yield from Southwest Idaho watersheds above 6000 feet 
elevation is mainly contributed during the spring snovmelit 
season. In contrast, watersheds from 4000 to 6000 fect eleva- 
tion contribute the greatest sediment during mid-vrinter runoff 
events which occur when soils are frozen or nearly saturateds. 
However, extremely dry or wet years may cause a change from the 
somewhat normal pattern, Delineation of major sediment source 
areas under the extremely variable topography, soil and cover 
conditions found in the Reynolds Creek ijatershed is very 
difficult without runoff-sediment plots to define the factors 
which greatly influence erosion and sediment production from 
specific sites. (Ida-So-107.1) 
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